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Abstract 


A discrete  collisional  theory  of  dissociation  under 
equilibrium  and  non-equilibrium  conditions  using  a quantum 
kinetic  approach  to  molecular  energy  transfer  was  investi- 
gated. Semi-classical  probabilities  of  vibrational 
exchange  based  on  a distorted  wave  treatment  were  calcu- 
lated for  the  molecular  species  of  interest.  Included  in 
this  analysis  was  the  specific  evaluation  of  the  contribu- 
tion of  rotation  and  the  source  of  rate  enhancement  in  pure 
gases . 

Using  an  experimentally  consistent  set  of  rates,  the 
kinetic  equations  for  molecular  dissociation  in  a pure 
diatomic  gas  and  diatomic  gases  in  the  presence  of  an  inert 
diluent  were  formulated  and  solved.  This  solution  estab- 
lished that  in  the  temperature  range  typical  to  the  dissoci- 
ation process,  a pseudo-steady  state  was  rapidly  attained. 

The  delay  or  incubation  time  prior  to  the  establishment  of 
steady-state  dissociation  was  calculated  and  exhibited  good 
agreement  with  experimental  measurements.  An  analytic 
steady  state  form  for  the  vibrational  population  and  dissoci- 
ation rate  was  developed  which  agreed  with  the  time- 
dependent  calculations.  The  influence  of  rotational  effects 
was  considered.  It  was  determined  that  rotational  states 
play  a significant  role  in  the  dissociation  process  and  that 
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the  standard  ladder  model  is  inadequate.  When  rotational 
effects  are  properly  included  in  the  analysis,  good  agree- 
ment with  experimental  data  is  obtained  for  diluted  mixes. 

An  analysis  of  dissociation  in  pure  gases  established 
that  the  observed  rate  enhancement  in  pure  gases  can  not 
be  attributed  to  changes  in  collisional  efficiency  due  to 
mass  variations,  however  V-V  exchanges  do  yield  a rate 
enhancement  consistent  with  that  observed  experimentally. 

The  kinetics  developed  to  address  thermal  dissociation 
under  equilibrium  conditions  were  extended  to  treat  the 
highly  non-equilibrium  conditions  achieved  through  intense 
optical  pumping  or  electron  impact  excitation.  In  addition 
to  those  physical  aspects  addressed  in  the  equilibrium 
treatment,  the  non-equilibrium  treatment  included  a 
detailed  investigation  of  the  electron  pumping  and  electron- 
collisionally  induced  dissociation.  This  aspect  involved 
the  self-consistent  solution  of  the  collisional  Boltzmann 
equation  including  excited  states  and  superelastic  colli- 
sions. A parametric  study  of  heavy  particle  and  electron 
impact  dissociation  rates  revealed  that  under  conditions 
typically  encountered  in  electric  discharges,  electron 
impact  dissociation  dominates  the  heavy  particle  rate. 
However,  under  highly  non-equilibrium  conditions,  the  heavy 
particle  rate  is  significantly  enhanced  over  its  equilibrium 
value. 
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Chapter  I_.  Introduction 


Purpose  and  Objective 
of  This  Study 

Dissociation  and  energy  transfer  in  diatomic  gases  is  a 
topic  embracing  a broad  class  of  physical  phenomena  and  has 
direct  application  to  a wide  variety  of  technical  fields. 

Gas  temperatures  of  a few  thousand  degrees  Kelvin  are 
encountered  in  the  operating  environment  of  re-entry  vehic- 
les, rocket  engines,  and  hypersonic  aircraft.  The  need  for 
a detailed  understanding  of  the  processes  and  phenomena 
associated  with  this  environment  has  spurred  interest  in 
the  chemical  kinetics  of  the  dissociation  and  the  recombina- 
tion of  molecules.  The  development  of  chemical,  electric 
discharge  and  gas-dynamic  lasers  has  accelerated  interest  in 
the  allied  field  of  intermolecular  energy  transfer  and  iden- 
tified the  dissociation  or  recombination  processes  as  funda- 
mental to  the  development  of  new  and  reliable  laser  systems. 
Our  interest  in  the  dissociation  process  grew  out  of  studies 
(ref  1)  of  the  characteristics  of  closed  cycle,  electric 
discharge  lasers.  Within  these  studies,  we  established  that 
plasma  chemistry  has  a critical  influence  on  discharge 
stability  and  laser  performance.  The  dissociation  reaction 
initiates  the  chemical  evolution  of  the  gas  discharge,  and 
is  therefore  fundamental  to  an  understanding  of  the  kinetics 
of  closed  cycle  lasers. 
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While  the  modern  literature  and  texts  contain  numerous 
studies  of  energy  transfer  and  dissociation  in  diatomic 
molecules,  recent  measurements  of  dissociation  (refs  2,  3, 
4)  and  relaxation  rates  (refs  5,  6,  7,  8)  indicate  substan- 
tial deviations  from  previous  theories.  Most  previous 
theories  of  dissociation  (refs  9,  10)  are  based  upon  a 
steady  state,  "ladder  climbing"  analysis  with  the  vibra- 
tional energy  of  the  dissociating  molecule  making  the  domi- 
nant contribution  to  the  reaction  energy.  Recent  detailed 
calculations  utilizing  this  theory  (refs  11,  12)  have  shown 
that  this  model  fails  to  characterize  the  process  consis- 
tently, and  yields  dissociation  rates  an  order  of  magnitude 
lower  than  observed  experimentally.  On  the  other  hand, 
several  recent  publications  (refs  11,  12,  13,  14,  15)  based 
on  a modified  "ladder"  analysis  have  met  with  varying 
degrees  of  success.  These  analyses  have  been  inconsistent 
with  regard  to  relaxation  rates,  and  yielded  contradictory 
results.  These  recent  works  have  included  a variety  of 
additional  processes  viz;  anharmonic  effects,  multiple 
quantum  transitions,  complex  exchanges  and  rotational 
effects.  However,  these  processes  have  been  treated  indi- 
vidually and  for  specific  conditions  so  that  detailed  com- 
parisons are  difficult. 

The  apparent  anomalies  and  inconsistencies  in  the 
theories  of  thermal  dissociation  hctve  not  tempered  the 
very  recent  extension  of  the  standard  "ladder"  model  to 
highly  non-equilibrium  conditions  (refs  16,  17,  18).  These 
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non-equilibrium  studies  have  been  confined  to  an  analytical 
steady  state  analysis  and  have  examined  the  dissociation 
reaction  rate  under  conditions  of  intense  optical  or  col- 
lisional  pumping.  Rotational  effects  and  the  influence  of 
electron  impact  dissociation  have  not  been  included  in  the 
analysis . 

This  study  is  addressed  to  establishing  a consistent, 
discrete  collisional  theory  of  dissociation  under  equili- 
brium and  non-equilibrium  conditions  and  to  five  correspond- 
ing objectives. 

The  first  objective  is  to  examine  the  standard  ladder 
model  of  dissociation  and  evaluate  its  predictive  utility 
when  realistic  transfer  rates  are  incorporated  into  the 
analysis.  A variety  of  gases  are  analyzed  over  a broad 
range  of  temperatures.  Anticipating  the  failure  of  the 
standard  "ladder”  model,  the  second  objective  is  to  evalu- 
ate the  role  of  rotational  effects  in  modifying  the  col- 
lisional activation  process  and  hence  to  determine  to  what 
extent  does  rotation  alter  the  calculated  dissociation 
rates.  While  rotational  effects  have  long  been  invoked  as 
leading  to  enhanced  dissociation,  few  serious  and  extensive 
analyses  have  been  performed.  Within  the  evaluation  of  the 
standard  and  rotational  models,  the  third  objective  will 
be  to  examine  the  time  dependent  characteristics  of  the 
dissociation  process  and  assess  the  influence  of  the  models 
on  the  observed  incubation  times.  The  fourth  objective  is 
to  analyze  the  effect  of  V-V  exchanges  and  of  changes  in 
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collisional  efficiencies  due  to  mass  variations.  This 
objective  will  address  the  resolution  of  the  inconsisten- 
cies that  have  arisen  in  the  literature  regarding  rate 
enhancement  in  pure  gas.  The  fifth  objective  is  then  to 
extend  the  analysis  to  the  highly  non-equilibrium  situation 
of  intense  optical  or  electron  impact  pumping  and  to  param- 
eterize the  dissociation  process  in  this  regime. 


In  this  chapter  we  begin  with  a review  of  collisional 
energy  exchange  mechanisms  in  molecules,  an  examination  of 
the  theories  of  vibrational  energy  transfer,  and  an  intro- 
duction to  the  dissociation  process. 

In  Chapter  II,  a historical  review  of  the  theories  of 
dissociation  is  presented  which  provides  a guide  to  a num- 
ber of  more  complete  review  papers.  The  physical  concepts 
and  critical  assumptions  of  previous  theories  are  then  dis- 
cussed. Here  the  dissociation  terminology  of  vibrational 
depletion,  incubation  time  and  vibrational  bias  is  intro- 
duced. Previous  calculations  are  reviewed  and  results 
cited.  The  contradictions  and  inconsistencies  present  in 
previous  works  are  then  highlighted  and  attention  is  pointed 
towards  the  need  for  the  present  analysis. 

Chapter  III  is  devoted  to  the  development  of  the  Master 
Equation  formalism  for  diatomic  gases  diluted  in  inert 
gases  and  pure  diatomic  gases.  The  numerical  technique 


used  for  the  time-dependent  solution  is  discussed  and  a 
steady  state  solution  derived. 
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Building  on  the  formalism  developed  in  Chapter  III, 
three  simplified  models  of  dissociation  are  presented  and 
examined  in  Chapter  IV.  Results  of  detailed  numerical  cal- 
culations of  dissociation  rates  based  on  these  models  are 
presented  in  Chapter  V and  compared  with  experimental  data. 
The  time-dependent  behavior  of  the  dissociation  process  is 
also  examined.  In  addition,  the  extent  of  rate  enhancement 
due  to  changes  in  collisional  efficiencies,  complex  energy 
exchanges  and  coupling  to  excited  electronic  states  are 
also  discussed. 

Chapter  VI  deals  with  the  highly  non-equilibrium  con- 
ditions achieved  through  electron  impact  excitation  or 
intense  optical  pumping.  In  order  to  address  electron 
impact  dissociation,  a solution  of  the  collisional  Boltzmann 
transport  equation  is  developed  and  the  electron  impact 
rates  are  examined  for  parameters  characteristic  of  elec- 
trical discharges.  Next,  the  characteristics  of  the  vibra- 
tional energy  distribution  are  discussed  within  the  frame- 
work of  an  analytical  steady-state  solution  of  the  Master 
Equation  for  strongly  pumped  conditions.  Time  dependent 
and  steady  state  results  are  then  presented  in  a parametric 
study  of  the  heavy  particle  and  electron  impact  dissocia- 
tion processes.  A comparison  is  made  of  the  efficiencies 
of  electron  impact  dissociation  and  heavy  particle  dis- 
sociation. Finally,  Chapter  VII  summarizes  the  new  aspects 
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of  the  results  of  this  study  and  presents  some  considera- 
tions and  recommendations  for  additional  theoretical  and 
experimental  study. 

Energy  Transfer  Mechanisms 

In  considering  the  dissociation  of  diatomic  molecules 
in  pure  gases  or  in  highly  diluted  mixtures  with  an  inert 
gas,  five  distinct  energy  transfer  mechanisms  resulting 
from  the  collisional  interaction  of  heavy  particles  can  be 
identified.  These  are: 

1.  Vibration-Translation  exchange  (V-T) 

2.  Rotation-Translation  exchange  (R-T) 

3.  Electronic-Translation  exchange  (e-T) 

4.  Vibration-Vibration  exchange  (V-V) 

5.  Electronic-Vibration  exchange  (e-V) 

A V-T  exchange  is  characterized  by  a collision  of  a 
molecule  with  another  atom/molecule  in  which  the  molecule 
gains  or  loses  vibrational  energy  from/to  the  translation 
mode.  The  R-T  and  e-T  processes  represent  a similar  exchange 
except  they  involve  the  transfer  of  rotational  and  elec- 
tronic energy  respectively.  V-V  and  e-V  exchanges  are 
classified  as  complex  exchanges.  AV-V  exchange  may  occur 
when  two  molecules  collide.  One  molecule  gains  vibrational 
energy,  the  other  molecule  loses  vibrational  energy,  and 
any  energy  mismatch  is  compensated  by  the  translational 
mode.  An  e-V  exchange  parallels  a V-V  exchange  except  that 
electronic  and  vibrational  energies  are  exchanged. 
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In  the  present  work,  it  will  be  assumed  that  the  rota- 
tional mode  is  in  equilibrium  with  the  translational  mode; 
and  we  will  therefore  omit  an  explicit  treatment  of  rota- 
tional energy  exchanges  during  collisions.  The  assumption 
of  rotational  equilibrium  is  easily  justified  for  low  lying 
rotational  states,  being  consistent  with  the  experimental 
observation  of  rapid  rotational  energy  relaxation  during 
collisions  (ref  24) . However,  when  very  large  values  of 
the  rotational  quantum  number,  J,  are  considered,  the  assump- 
tion is  questionable.  We  retain  this  assumption;  however, 
based  on  the  simplification  it  introduces  and  the  lack  of 
any  experimental  or  theoretical  exchange  rate  data  for  very 
high  J levels.  e-T  exchanges  will  also  be  neglected.  Their 
omission  is  justified  by  the  adiabatic  nature  of  such  an 
exchange  in  the  temperature  range  considered.  The  explicit 
evaluation  of  e-V exchanges  is  also  omitted  in  the  present 
analysis.  However,  the  possible  role  played  by  e-V  exchange 
in  the  dissociation  process  will  be  discussed  when  we 
examine  dissociation  under  non-equilibrium  conditions  in 
’ Chapter  VI.  Thus  in  this  work,  we  will  concentrate  on  the 
V-T  and  V-V  mechanisms  as  we  proceed  with  the  development 
of  a collisional  theory  of  dissociation. 

Theory  of  Vibrational 
Energy  Exchange 


Having  identified  the  critical  energy  transfer  mech- 
anisms as  V-V  and  V-T  exchange,  we  proceed  by  outlining  the 
development  of  the  theory  of  vibrational  energy  exchange. 


Landau  and  Teller  (ref  19)  first  examined  the  theory  of 
vibrational  relaxation  in  1936  using  a classical  mechanical 
approach.  Schwartz,  Slawsky,  and  Herzfeld  (ref  20)  published 
an  analogous  quantum  mechanical  theory  in  1952.  This  "SSH" 
theory  has  received  widespread  use  and  represents  the  stan- 
dard with  which  most  relaxation  data  are  compared.  The  SSH 
analysis  is  based  upon  an  adiabatic  "distorted  wave" 
quantum  mechanical  theory  of  scattering  developed  by  Jackson 
and  Mott  (ref  21) . Since  the  appearance  of  the  "SSH" 
theory,  an  extensive  amount  of  theoretical  work  on  colli- 
sional  energy  exchange  has  been  accomplished.  This  includes 
both  rotational  and  vibrational  energy  exchange.  This  work 
has  been  reviewed  in  the  landmark  work  of  Herzfeld  and 
Litovitz  (ref  22) , and  by  Takayanagi  (ref  23) , Cottrell  and 
McCoubrey  (ref  24) , Rapp  (ref  25) , and  Nikitin  (ref  26) . 

Subsequent  improvements  in  the  basic  SSH  theory 
centered  around  the  anharmonic  features  of  the  molecular 
potential.  Rich  (ref  14)  improved  the  SSH  approach  by 
using  an  "exact"  form  of  the  Morse  transition  probabilities. 
Mies  (ref  27)  identified  the  critical  dependence  of  the 
transition  probability  upon  the  ratio  of  the  diagonal 
matrix  elements.  In  1965  Keck  and  Carrier  (ref  28)  were 
able  to  extend  the  SSH  formalism  through  a numerical  calcu- 
lation and  curve  fitting  procedure  to  provide  a smooth 
transition  from  the  adiabatic  to  non-adiabatic  limits  of 
the  V-T  exchange.  In  1970,  Rich  and  Treanor  (ref  29)  pre- 
sented a comprehensive  review  of  vibrational  relaxation  in 
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gasdynamic  flows.  The  early  seventies  brought  in  wide- 
spread efforts  to  deal  with  new  laser  requirements  for  the 
analysis  of  vibrational  energy  transfer  in  CO  and  CO^  laser 
systems.  At  this  time  detailed  numerical  calculations  of 
V-V  and  V-T  rates  started  to  surface.  These  calculations 
have  been  described  by  Secrest  (ref  30)  who  also  reviewed 
their  application  to  rotational  and  vibrational  energy 
transfer.  Numerical  calculations  have  been  a useful  gauge 
of  the  validity  of  the  various  approximations  in  the  ana- 
lytic theories.  These  assumptions  include  those  of  an 
exponential  interaction  potential  and  the  applicability  of 
first  order  perturbation  theory.  The  recent  fully  quantum 
calculations  of  Rabitz  (ref  31)  and  the  semi-classical  cal- 
culations of  McKenzie  (ref  32)  represent  the  state  of  the 
art.  These  calculations  themselves,  however,  have  been 
quite  limited  in  number  of  gases  and  temperature  ranges 
considered.  This  is  easily  understood  since  these  analyses 
are  quite  lengthy  from  a computational  standpoint  and  are 
limited  by  the  accuracy  of  the  interaction  potential  surface 
used . 

When  applying  collisional  energy  exchange  theory  to 
the  problem  of  molecular  dissociation,  three  aspects  of  the 
exchange  rates  must  be  considered. 

1.  Temperature  dependence 

2.  Magnitude 

3.  Variation  wi  vibrational  state 
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The  vibrational  level  scaling  of  the  rates  reflects 
the  influence  of  the  anharmonic  nature  of  the  vibrational 
energy  levels  on  the  calculated  transition  probabilities. 
When  we  consider  highly  excited  states  the  consequent 
limitations  on  the  validity  of  perturbation  and  simple 
harmonic  oscillator  theories  must  be  realized  and  can  be 
accounted  for  by  using  Morse  transition  probabilities  and 
the  Mies  (ref  27)  assymetrical  correction  factor.  For  a 
more  detailed  discussion  see  Appendix  A.  In  addition  to  the 
influence  of  molecular  anharmonicity  on  the  transition  ele- 
ments, the  variation  in  the  spacing  of  the  vibrational 
energy  levels  also  changes  the  nature  of  the  collisional 
interaction.  Due  to  this  variation  in  energy  level  spacing, 
collisions  involving  low  lying  vibrational  levels  are 
adiabatic  in  character  while  collisions  involving  vibra- 
tional states  near  the  dissociation  limit  are  non-adiabatic 
or  impulsive.  Experimental  data  on  excited  state  V-T  mea- 
surements are  limited,  but  comparison  with  the  available 
data  of  Hancock  and  Smith  (ref  33)  show  favorable  agreement 
with  the  SSH  theory  as  modified  by  Keck  and  Carrier.  Based 
on  this  comparison  with  experiment  and  numerical  calcula- 
tions of  McKenzie  (ref  32),  the  SSH  theory  as  modified  by 
Keck  and  Carrier  (ref  28)  will  be  utilized  throughout  the 
present  analysis  and  fitted  to  experimental  data  where  avail- 
able. 
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The  Dissociation  Process 


The  previous  sections  have  dealt  with  the  physics  of 
collisional  energy  transfer  and  vibrational  energy 
exchange.  The  energy  transfer  mechanisms  of  V-V  and  V-T 
exchange  have  been  defined  and  discussed.  As  the  vibra- 
tional excitation,  relaxation,  and  excitation  sequence  pro- 
ceeds, excitation  to  higher  and  higher  vibrational  levels 
occurs.  If  this  excitation  proceeds  sufficiently  far, 
molecular  bonding  is  overcome  and  two  free  atoms  are  formed 
the  molecule  dissociates.  Thus  if  the  energy  of  a diatomic 
molecule  exceeds  the  dissociation  energy,  the  molecule  will 
decompose  at  a frequency  on  the  order  of  vibrational  fre- 
quency. This  frequency  is  much  greater  than  the  gas 
kinetic  collision  frequency  and  therefore  this  phase  of  the 
reaction  may  be  considered  as  instantaneous.  As  a result, 
the  effective  reaction  rate  is  determined  by  the  rate  of 
excitation  to  the  dissociating  state.  Thus  in  order  to 
investigate  the  dissociation  process  and  thermal  activation 
it  is  essential  to  know  the  behavior  of  the  collision 
induced  transition  probabilities. 

Experimental  evidence  from  resonance  emission  spectro- 
scopy (refs  34,  35,  36)  and  theoretical  calculations  (refs 
31,  32)  indicate  that  energy  transfers  involving  the 
exchange  of  more  than  one  or  two  vibrational  quanta  are 
extremely  improbable.  Whereas,  the  probability  of  popu- 
lating the  vibrational  states  from  those  immediately  below 
has  a relatively  high  value.  This  restriction  on  the 


number  of  quanta  exchanged  in  a collision  leads  one  to  view 
the  excitation  of  the  vibrational  states  during  relaxation 
as  a stepwise  process.  The  activation  of  excited  vibra- 
tional states  therefore  occurs  primarily  through  single 
quantum  exchanges  and  the  overall  energization  of  the  vibra 
tional  manifold  can  be  characterized  as  a stepwise  process 
in  which  excitation  proceeds  up  the  "ladder"  of  vibra- 
tional energy  levels  to  high  lying  states  from  which  dis- 
sociation occurs.  This  concept  of  stepwise  excitation  to  a 
dissociating  state  is  the  essence  of  the  "ladder  model"  of 
dissociation  and  relaxation. 

The  collisionally  induced  dissociation  of  diatomic 
molecule,  A2,  can  be  represented  as  a simple  chemical  pro- 
cess 

A —*■  A+A (i-D 

The  rate  equation  for  the  reaction  (1)  takes  the  form 

(1-2) 

where  the  number  densities  of  the  particles  have  been 
denoted  by  their  symbols.  The  reaction  rates  depend  on 
the  type  of  molecule  or  atom,  and  gas  temperature.  In  pre- 
dicting the  dissociation  reaction  races  from  the  properties 
of  the  reacting  molecules  and/or  at  >ms,  one  is  concerned 
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with  the  evaluation  of  the  dissociation  rate  coefficient 

k. . The  recombination  rate  coefficient  k is  related  to 
d r 

kd  through  the  law  of  mass  action  or  rate  quotient  law: 


where  K (T)  is  the  equilibrium  constant.  Due  to  the  develop- 

G 

ment  of  the  shock  tube  as  a reliable  tool  for  high  tempera- 
ture kinetic  studies,  extensive  data  exist  on  the  dissoci- 
ation coefficients  of  H 2,  0 2,  N2,  Br2 , *2  and  C0  ^ases  diluted 
in  an  argon  mixture.  Data  are  also  available  for  pure 
diatomic  gases,  although  these  data  are  not  quite  as  extensive. 
The  dissociation  rate  in  pure  gases  is  generally  enhanced 
over  that  value  obtained  when  the  diatomic  gas  is  diluted 
in  an  inert  gas . The  magnitude  and  temperature  dependence 
of  the  enhancement  varies  from  gas  to  gas.  We  will  con- 
sider the  possible  sources  of  this  enhancement  in  Chapter  V. 

The  results  of  experimental  measurements  of  the  dis- 
sociation coefficient  are  often  fitted  to  the  form 


^ = C/fmexp  ) 


where  8=1/  (kT),  T is  the  translational  temperature,  Dq  the 
dissociation  energy,  and  C a constant.  The  Arrenhius 
activation  energy,  Eact  , is  given  by 


13 


~j5  c^1  n 


(1-5) 


Thus  for  the  form  given  in  (1-4) , the  activation  energy  is 


£act.  = Q * m/f 
= ZL  + 7~ 


(1-6) 


Table  I,  extracted  from  Johnston  and  Birks  (ref  9),  pre- 
sents experimental  data  on  the  observed  activation  ener- 
gies, Eobs,  and  the  power  of  the  temperature,  m,  for  a 
given  reactant  (A),  gas  diluent  (X),  and  temperature  range. 
Note  that  the  observed  activation  energies  are  less  than 
the  dissociation  energies  so  that  the  empirical  parameter 
m is  always  negative.  This  fact  serves  as  a guide  in 
evaluating  the  various  models  of  dissociation  that  we  will 
introduce  in  Chapter  IV. 
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TABLE  1 


EXPERIMENTAL  DATA  ON  ACTIVATION  ENERGIES  FOR  THE 
DISSOCIATION  OF  HOMONUCLEAR  DIATOMIC  MOLECULES 


A 

X 

T,K 

■.  ~ -r:  -.i:..  — — - . g= 

Do, kcal 

Eobs.,koal 

m 

F2 

Ar 

1300-1600 

37.1 

30  + 4 

-2.5 

Ne 

1650-2700 

24  + 5 

-3.0 

Ar 

1300-1600 

27.3  + 2.5 

-3.4 

Xe 

1300-1600 

31.1 

-2.1 

C12 

Ar 

1700-2500 

57.0 

48.3 

-2.1 

1600-2600 

50 

-1.7 

1700-2600 

41  + 5 

-4.0 

1600-2600 

45  + 2 

-2.9 

1700-2600 

48.3 

-2.0 

Br2 

Ar 

1300-1900 

45.5 

41.4 

-1.3 

1400-2700 

41.4 

-1.0 

1200-2200 

32.4 

-3.8 

X2 

Ar 

1000-1600 

35.6 

29.7 

-2.3 

850-1650 

30.4 

-2.1 

°2 

Ar 

5000-18000 

118.0 

110 

-0.4 

3800-5000 

106  + 5 

-1.4 

4000-6000 

108 

-1.0 

3400-7500 

108 

-1.2 

3000-5000 

114 

-0.5 

°2 

3000-5000 

98 

-2.5 

4000-7000 

91 

-2.5 

mm 
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TABLE  1 — Continued 


A 

X 

T,K 

DQ,kcal 

E , ,kcal 

obs. 

m 

2500-4000 

108 

-1.5 

3000-6000 

100  + 5 

-1.0 

2600-7000 

85 

-3.7 

11 

H2 

Ar 

2800-4500 

103.3 

97 

-0.9 

i 

3000-5300 

97 

-0.8 

2800-5000 

97 

-0.8 

: 

| 

2900-3790 

97 

-1.0 

Xe 

3000-4500 

100 

-0.4 

H2 

3000-4500 

*> 

92 

-1.5 

H 

3000-4500 

100 

-0.4 

°2 

Ar 

3000-4800 

105 

97 

-1.0 

M J 

Ar 

3000-4900 

97 

-1.0 

N2 

Ar 

6000-9000 

225 

218 

-0.5 

1 

6000-10000 

204 

-1.3 

N2 

6000-9000 

218 

-0.5 

i 

. 

6000-10000 

198 

-1.7 

1 

N 

6000-9000 

203 

-1.5 

NOTES : 

Dq — Dissociation  Energy 


E , — Observed  activation  energy 


Chapter  II.  Dissociation  Theories; 


Approaches  to  the  Problem 

Having  discussed  vibrational  energy  transfer  and  intro- 
duced the  dissociation  process,  let  us  turn  to  previous  work 
on  dissociation  theories  and  review  the  approaches,  assump- 
tions, and  results.  It  is  convenient  to  divide  these 
approaches  into  three  groups:  the  classical  approach,  which 
treats  the  molecular  energy  levels  as  a continuum  and  gives 
no  reference  to  the  quantum  nature  of  the  activation  or 
reaction  processes;  the  statistical  approach,  where 
reactions  are  considered  as  systems  crossing  boundaries  in 
a multi-dimensional  phase  space;  and  finally,  the  colli- 
sional  approach,  involving  a discrete  kinetic  analysis  of 
molecular  activation  and  reaction.  We  will  consider  these 
approaches  in  turn. 

Classical  Theory 

The  early  theories  of  the  dissociation  of  diatomic 
molecules  used  classical  mechanics  and  kinetic  theory. 

Within  this  framework,  it  was  assumed  that  in  order  for  a 
molecule  to  dissociate  as  a result  of  a binary  encounter, 
the  relative  kinetic  energy  of  the  collision  partners  must 
only  be  in  excess  of  some  activation  energy.  Under  this 
"Available  Energy  Theory"  the  specific  rate  of  the  reaction 


Az  + X-~  4+  X'X 


(2-1) 


at  a temperature  T,  is  given  by  (ref  37) 

<^  = TZ  exp  (-/$£),)  (2-2) 

where  Z is  the  collision  rate,  P a steric  factor,  B=l/kT 
and  Dq  is  the  dissociation  energy.  This  approach  has  been 
extended  to  include  the  influence  of  internal  energy  in  all 
participating  modes.  At  dissociation  temperatures,  mole- 
cules will  have  significant  rotational  and  vibrational 
energy.  Assuming  equipartition  of  this  energy  over  2s 
classical  dynamical  variables  and  classical  excitation  of 
these  modes,  the  rate  constant  is  given  by  (ref  37): 

-i  ' '*/>(-/&•)  <2'3' 


This  form  is  similar  to  that  discussed  previously  in 
Chapter  I and  is  frequently  used  to  correlate  experimental 
data  as  presented  in  Table  I.  Rice  (refs  38-39)  attempted 
to  refine  this  theory  by  assuming  that  on  the  average  a 
molecule  that  is  about  to  dissociate  upon  a collision  will 
have  an  energy  within  kT  of  the  dissociation  energy.  This 
assumption  coupled  with  a correction  for  electronic  state 


E I 

multiplicity  resulted  in  a significant  increase  in  the 
effective  Z in  (2.3).  Careri  (refs  40-41),  however,  cited 
Rice's  omission  of  rotation  and  upon  including  rotational 
effects  he  obtained 

Jd  - PZ  O'  exp(/#/»')]exp  C-pA.) 

[fr~PDo  (jgjt?  ' G/^)]  (2-4) 
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where  v is  the  vibrational  frequency  of  the  ground  state. 
Rice  (ref  39)  criticized  this  result.  Although  agreeing 
that  rotation  should  be  included,  he  pointed  out  that  add- 
ing the  rotational  term  to  the  potential  energy  curve  would 
produce  a weak  maximum  at  some  inter-nuclear  separation 
greater  than  the  equilibrium  separation.  This  maximum  is 
known  as  the  rotational  barrier.  Rice  stated  that  the  bar- 
rier would  result  in  a reduction  in  the  dissociation  rate. 

No  quantitative  evaluation  of  this  effect  was  made  in  the 
paper. 

In  addition  to  raising  the  question  of  rotational 
effects,  the  theories  of  Rice  and  Careri  have  resulted  in 
the  formation  of  two  different  theoretical  "camps"  regarding 
vibrational  bias  or  how  dissociation  is  distributed  over  the 
vibrational  levels.  Rice  calculated  the  dissociation  rate 
by  assuming  a strong  bias  of  the  dissociation  reaction 
toward  high  vibrational  states.  He  assumed  that  dissociation 
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proceeded  dominantly  from  vibrational  states  within  kT  of 
the  continuum  and  that  it  was  necessary  for  molecules  to 
climb  the  "ladder"  of  vibrational  states  from  the  lower 
states  to  these  upper  levels  prior  to  undergoing  a colli- 
sion resulting  in  dissociation.  Conversely,  Careri  felt 
that  the  dissociation  process  was  more  evenly  distributed 
over  the  various  vibrational  states  (that  is,  had  a weaker 
vibrational  bias) . In  his  view,  dissociation  was  likely  to 
occur  by  direct  transitions  from  low  vibrational  levels  as 
well.  We'll  return  to  this  point  when  we  review  the  colli- 
sional  approach,  but  at  this  time  state  that  this  question 
of  bias  has  yet  to  be  resolved  and  the  two  "camps"  remain 

i 

separated. 

It  is  easy  to  cite  the  limitations  of  the  classical 
approach  especially  if  we  are  considering  highly  non- 
equilibrium  circumstances.  The  assumptions  of  equilibrium 
and  of  every  collision  possessing  sufficient  energy  result- 

V 

ing  in  reaction  are  questionable.  The  rate  of  collision- 
ally  induced  dissociation  should  be  amenable  to  a calcula- 
tion based  on  the  detailed  mechanics  of  the  collision  with- 
out invoking  such  strong  and  limiting  assumptions.  The 
task  of  establishing  a more  detailed  analytic  theory  is 
certainly  worthwhile  and  has  been  pursued  utilizing  both 
the  statistical  and  the  collisional  approaches. 

Phase  Space  Theories 

The  statistical  theories  approach  the  problem  of  dis- 
sociation and  recombination  by  reasoning  that  near  the 
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dissociation  limit,  transitions  between  states  are  very 
rapid  and  a classical  approximation  is  valid.  Keck  (refs 
42,  43,  44)  states  that  because  of  the  high  density  of 
levels  near  the  dissociation  limit,  the  classical  Master 
Equation  is  appropriate  for  treating  dissociation  and 
recombination.  In  the  limit  of  small  energy  transfers  this 
integro-dif ferential  equation  may  be  reduced  to  an  equiva- 
lent diffusion  equation.  Reactions  are  then  considered  as 
systems  crossing  boundaries  and  described  as  the  motion  of 
a representative  point  in  an  18  dimensional  phase  space. 
Wigner  (refs  45,  46,  47)  pioneered  research  in  this  area  by 
considering  the  phase  space  volume  accessible  to  the  three 
recombining  atoms.  The  phase  space  is  divided  into  a 
recombined  region,  where  two  atoms  have  recombined  in  the 
presence  of  a third  body,  and  a free  region.  The  division 
into  two  regions  is  based  on  the  fact' that  in  the  recombined 
region  the  relative  kinetic  energy  is  negative  and  in  the 
free  region  it  is  positive.  The  trajectories  of  the  atoms 
represent  lines  in  phase  space  and  the  reaction  rate  is 
represented  by  the  number  of  such  trajectories  entering  the 
recombined  region.  This  analysis  only  yields  an  upper  bound 
to  the  recombination  rate  since  it  neglects  recrossing  of 
the  phase  boundary  surface.  Some  ambiguity  may  exist  in 
definition  of  the  recombined  region.  In  order  to  remove  the 
ambiguity,  the  recombined  region  of  phase  space  is  taken  as 
a "trial  surface,"  which  is  then  varied  subject  to  the  con- 
dition that  the  rate  constant  obtained  is  a minimum  and  a 
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least  upper  bound  is  obtained.  Refinements  of  the  theory 
over  the  past  10  years  have  included  the  evaluation  of  a 
recrossing  parameter  and  a correction  for  departures  from 
equilibrium.  Suitable  choice  of  adjustable  parameters  in 
the  theory  has  yielded  rate  constants  of  the  correct  order 
of  magnitude  and  temperature  dependence  for  a variety  of 
diatomic  molecules  (ref  109) . 

Collisional  Approach 

A collisional  approach  to  dissociation  is  undertaken  in 
the  present  work.  Adoption  of  this  approach  is  encouraged 
by  the  rapid  advances  in  vibrational  kinetic  studies 
relating  to  laser  technology  and  the  new  availability  of 
experimental  data  on  non-equilibrium  vibrational  distribu- 
tions and  kinetic  rates.  Moreover,  the  collisional  approach 
offers  the  unique  prospect  of  treating  the  highly  non- 
equilibrium situations  of  intense  optical  pumping  and  ana- 
lyzing the  new  generation  of  excimer  laser  systems. 

Several  detailed  collisional  theories  of  molecular  dis- 
sociation have  surfaced  over  the  past  twenty  years.  As 
brought  out  in  Chapter  I,  since  the  decomposition  of  the 
molecule  can  be  assumed  to  be  instantaneous,  the  observed 
dissociation  rate  is  governed  by  the  rates  of  thermal  acti- 
vation and  energy  exchange  of  the  molecule.  We  can  investi- 
gate this  thermal  activation  providing  we  know  the  transition 
rates,  V-V  and  V-T,  of  the  various  vibrational  states  of 
the  molecule.  The  calculation  of  the  dissociation  rate  is 
then  based  on  the  solution  of  the  Master  Equation,  a system 


of  coupled  kinetic  equations  describing  transitions  between 
all  the  bound  and  free  states  of  the  diatomic  molecule.  The 
Master  Equation  formalism  will  be  developed  in  the  next 
chapter.  Here  we  review  previous  collisional  approaches  to 
the  dissociation  problem.  The  main  differences  in  the  vari- 
ous calculations  relate  to  the  internal  modes  considered, 
the  mechanisms  of  energy  exchange,  and  the  vibrational  bias 
of  the  dissociation  reaction.  These  differences  lead  to 
variations  in  the  magnitude  and  temperature  dependence  of 
the  calculated  rates. 

First  let  us  consider  those  works  concerned  with  a 
single  quantum  stepwise  excitation.  The  early  work  of 
Nikitin  (refs  48,  49,  50)  treated  the  dissociation  of 
diatomic  molecules  in  the  presence  of  an  inert  diluent.  The 
molecule  was  represented  as  non-rotating  Morse  oscillator 
and  each  molecular  vibrational  state  was  treated  as  a dis- 
tinct species.  Within  this  analysis  he  established  the 
importance  of  changes  in  adiabaticity  with  respect  to 
vibrational  transitions  and  demonstrated  that  most  of  the 
steps  in  the  vibrational  excitation  process  are  near  adia- 
batic, occurring  with  small  probability  for  large  exchanges 
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in  vibrational  energy.  In  the  temperature  range  of  interest 
for  most  studies  of  dissociation,  it  was  determined  that 
only  those  levels  within  kT  of  the  dissociation  energy  were 
definitely  non-adiabatic.  In  this  region,  collisions  are 
impulsive  and  transitions  occur  with  a probability 
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approaching  unity.  The  following  expression  for  the  dis- 
sociation rate  was  obtained 


j 

I 

1 


= Z ^ e*p(fE*)/Px,K-l  J (2-5) 

« f 

where  denotes  the  vibrational  energy  of  levels  v=0  to  d, 
and  P„  v . the  V-T  probability.  Upon  replacing  the  summa- 
tions  with  an  integration  he  obtained 

^ = Z ) O'  exP  (:^LU)J  (2-6) 

revealing  an  activation  energy,  Eactf  less  than  the  dissocia- 
tion energy,  Dq.  The  dissociation  rate  was  then  averaged 
over  all  rotational  states  assuming  rotational  equilibrium 
and  making  allowance  for  the  fact  that  the  depth  of  the 
radial  potential  will  decrease  with  increasing  angular 
momentum.  Application  of  this  analysis  was  limited  to  Br^ 
at  one  temperature. 

In  a subsequent  paper  by  Stupochenko  (ref  51) , the 
Master  Equation  was  solved  under  the  assumption  of  the 
existence  of  a steady  state  in  all  vibrational  states.  He 
concluded,  as  Nikitin  had,  that  a strong  vibrational  bias 
existed  for  the  dissociation  reaction.  In  addition,  the 
existence  of  vibrational  depletion  was  discussed.  Vibra- 
tional depletion  results  from  the  coupling  of  vibrational 
and  dissociative  kinetics  and  is  manifested  in  a reduction 
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in  the  population  of  the  vibrational  states  from  their 
equilibrium  value.  This  depletion  is  most  significant  for 
highly  excited  vibrational  states.  Stupochenko  provided  no 
detailed  calculation  of  transfer  rates  or  quantitative  cal- 
culations of  the  dissociation  rate. 

These  previous  works  favor  the  "ladder"  model  concept 
and  conclude  that  the  rate  determining  step  in  the  dissocia- 
tion reaction  is  the  activation  of  lower  states.  Concurrent 
with  these  publications  were  the  works  of  Montroll  and 
Shuler  (refs  52,  53,  54,  55,  56)  in  which  they  considered 
the  vibrational  relaxation  of  molecular  systems.  They  con- 
sidered a system  of  simple  harmonic  and  later  anharmonic 
oscillators.  In  the  latter  investigation  (ref  54) , the 
solution  was  achieved  through  a numerical  solution  of  the 
Master  Equation.  In  all  circumstances,  they  concluded  that 
single  quantum  transitions  are  most  likely  and  that  relaxa- 
tion and  dissociation  processes  are  "ladder-like."  Cal- 
culated dissociation  rates  (ref  55) , however,  were  several 
orders  of  magnitude  lower  than  experimental  values. 

Pritchard  (ref  57)  performed  a quantitative  investigation 
of  the  kinetics  of  dissociation  of  a diatomic  gas  in  which 
he  numerically  calculated  Morse  matrix  elements  for  use  in 
the  Jackson-Mott  (ref  21)  transition  probabilities  for 
hydrogen.  Bound-bound  and  bound-free  transitions  were 
evaluated.  He  showed  that  single  step  transitions  dominate 
and  that  a depopulation  or  depletion  of  high  vibrational 
levels  would  occur. 
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Benson  and  Fueno  (ref  58)  modeled  the  recombination 
process  assuming  a cascaded,  stepwise  process.  Landau-Teller 
(ref  19)  transition  probabilities  based  on  simple  harmonic 
oscillator  matrix  elements  were  used.  Rotational  effects 
were  accounted  for  by  introducing  a long-range,  polarization 
term,  to  the  potential.  A steady-state  solution  was 
obtained  and  the  recombination  coefficient  calculated  for  a 
variety  of  gases.  Order  of  magnitude  agreement  was  obtained. 

Treanor  and  Marrone  (refs  59,  60)  analyzed  the  effect 
of  vibrational-dissociation  coupling  using  a modified  form 
of  the  Landau-Teller  relaxation  equation.  They  exercised 
great  care  in  incorporating  experimentally  determined 
relaxation  times  and  equilibrium  rate  constants.  Various 
empirical  biases  for  the  dissociation  reaction  were  investi- 
gated ranging  from  an  Available  Energy  approach  to  a more 
strongly  biased  dissociation  probability.  The  results  of 
this  analysis  supported  the  ladder  model  to  a degree,  but 
suggested  a range  of  participating  vibrational  levels  rather 
than  a single  transition  occurring  from  the  last  bound 
vibrational  level. 

At  this  point  in  the  historical  development  of  colli- 
sional  rate  theories,  the  "ladder"  model  had  obtained  a 
strong  vote.  Most  analyses,  however,  were  quite  limited 
and  qualitative.  Vibrational  structure  of  the  ground  elec- 
tronic state  was  the  by-word  with  little  critical  attention 
being  devoted  to  incorporating  a detailed  experimentally  con- 
sistent set  of  V-T  ra.tes,  examining  the  role  of  rotation 


26 


and  electronically  excited  states,  or  examining  the  effect 
of  complex  collisions.  Time  dependent  solutions  of  the 
Master  Equation  were  in  the  wings  awaiting  better  numerical 
methods  and  faster  computers.  We  will  now  discuss  the  "New 
Era"  of  collisional  theories  during  which  significant 
strides  were  taken  to  address  the  areas  cited  above. 

Although  published  in  1959,  Herzf eld's  (ref  61)  work 
on  molecular  dissociation  revealed  an  excellent  perspective. 
The  paper,  although  marred  by  several  errors  which  obviate 
his  quantitative  results,  attempted  to  utilize  a realistic 
set  of  kinetic  rates  based  on  the  SSH  theory  using  Morse 
oscillator  wave  functions  and  an  exponential  interaction 
potential.  In  addition,  he  specifically  considered  the 
effect  of  V-V  transfer  on  the  dissociation  of  a pure  gas. 

The  influence  of  rotational  states  on  the  dissociation 
process  was  discussed  by  Bauer  and  Tsang  (ref  62)  in  a 
thought-provoking  article.  Here  they  called  attention  to 
several  mechanisms  which  effectively  couple  translational 
and  vibrational  states  through  the  intermediary  of  chemical 
reactions.  They  proposed  that  the  transfer  of  energy 
involving  both  rotation  and  vibration  at  approximately  con- 


stant total  energy  can  be  induced  by  distant  collisions. 
Consideration  of  such  collisions  gave  new  importance  to  the 
rotational  structure  and  provided  alternative  paths  for 


rapid  vibrational  relaxation.  Bauer  and  Tsang  hypothesized 
that  this  may  account  for  the  fact  that  the  observed  dis- 
sociation rates  are  considerably  higher  than  those  calculated 


under  the  assumptions  of  the  standard  ladder  model.  The 
fact  that  detailed  calculations  of  dissociation  rate  for  a 
variety  of  gases  using  the  "ladder"  were  not  completed  at 
this  time  left  this  hypothesis  unproved. 

Snider  (ref  63)  considered  the  effect  of  angular  momen- 
tum on  the  dissociation  rate  in  diatomic  molecules  assuming 
rotational  equilibrium.  Using  a classical  approach,  he 
analyzed  the  influence  of  rotation  on  the  mechanics  of  the 
collision.  His  results  did  not  differ  substantially  from 
those  of  Nikitin  (ref  50) . 

A comprehensive  study  of  the  "ladder"  model,  vibra- 
tional depletion,  and  V-V  exchange  was  presented  by  Rich 
(ref  64).  Following  the  work  of  Herzfeld  (ref  61),  he  used 
the  Jackson  and  Mott  (ref  21)  near-adiabatic  distorted  wave 
treatment  with  an  exponentially  repulsive  interaction  poten- 
tial in  an  approximate  form  that  enabled  him  to  obtain  a 
simple  analytic  form  for  the  matrix  elements.  Under  the 
approximation  used,  a single  quantum  selection  rule  emerged, 
while  still  retaining  the  effects  of  anharmonicity  in  the 
matrix  elements.  This  form  was  used  in  an  SSH  calculation 
of  the  transition  probabilities.  A time  dependent  solution 
of  the  coupled  rate  equations  was  obtained  in  addition  to  an 
analytical  expression  for  the  quasi-steady  state  dissocia- 
tion rate.  Good  agreement  was  obtained  with  experimental 
data,  however  the  Mies  assymmetrical  correction  factor 
(Appendix  A) , was  not  properly  accounted  for  resulting  in 
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a serious  overestimation  of  the  transition  probabilities.* 
Rich  proposed  that  V-V  exchange  accounted  for  the  observed 
rate  enhancement  in  pure  gases  and  calculated  the  cross- 
section  for  resonant  V-V  exchange.  Based  on  the  magnitude 
of  the  calculated  resonant  exchange  cross-section,  he  con- 
cluded that  it  was  sufficient  to  account  for  the  rate 
enhancement. 

Dove  (ref  65)  investigated  the  vibrational  relaxation 
and  dissociation  of  hydrogen  diluted  in  an  argon  heat  bath 
by  performing  a numerical  integration  of  the  coupled  rate 
equations.  Collisional  transition  probabilities,  for  V-T 
exchange,  were  calculated  by  the  SSH  method  modified  for  a 
Morse  oscillator.  The  time-dependent  Master  Equation  was 
numerically  integrated  and  vibrational  depletion  discussed. 
After  correcting  his  transition  probabilities  to  agree  with 
experimentally  determined  relaxation  data,  the  calculated 
dissociation  rates  were  a factor  of  30  too  low.  Multiple 
quantum  transitions  were  later  included  in  the  analysis  and 
resulted  in  an  increase  in  the  dissociation  rate  of  less 
than  15%.  Moreover,  the  computed  activation  energy  was  too 
large.  In  a subsequent  paper  by  Dove  (ref  66) , the  ladder 
approach  was  modified  to  take  explicit  account  of  rotational 
effects  using  the  matrix  elements  of  Roberts  (ref  67).  This 
distributed  the  dissociation  among  the  vibrational  levels 
and  resulted  in  an  increase  in  the  dissociation  rate  and  a 

♦The  correction  factor  was  assumed  to  be  independent  of 
vibrational  state. 
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reduction  in  the  activation  energy.  This  calculation 
yielded  good  agreement  with  experimental  data  both  in  rate 
magnitude  and  temperature  dependence. 
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The  relationship  of  the  observed  activation  energy  to 
vibrational  bias  was  examined  by  Johnston  (ref  68) . Noting 
that  the  experimentally  observed  activation  energy  is  always 
less  than  the  dissociation  energy,  and  decreases  with 
increasing  temperature,  he  analyzed  the  consistency  of  vari- 
ous models  of  dissociation  with  respect  to  these  criteria  in 
a steady  state  analysis.  He  concluded  that  a ladder- 
climbing model  for  the  reaction  process  with  dissociation 
occurring  only  from  the  top  vibrational  level  was  unaccept- 
able and  that  diatomic  molecules  dissociate  from  all  levels. 

The  influence  of  V-V  transfer  on  the  dissociation  rate 
was  investigated  by  Kiefer  (ref  69) . He  used  a steady-state 
approximation  for  a single-quantum  step  Master  Equation,  and 
concluded  that  the  rate  was  only  slightly  enhanced  when  V-V 
transitions  were  included.  Using  the  "ladder"  model  he 
achieved  "quite  satisfactory  agreement"  for  dissociation 
and  recombination  rates  in  Oj , H ^ and  HC1  and  provided 
"qualitative  explanations  for  nearly  all  the  unusual  fea- 
tures" of  the  measured  rate.  No  rotational  effects  were 
included. 

(As  you  can  now  see,  contradictions  and  inconsistencies 
have  arisen  in  the  literature  cited  regarding  the  appli- 
cability of  the  ladder  model,  the  role  of  rotation,  and  the 
influence  of  V-V  effects.  We  will  continue  with  this  review 
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of  the  collisional  theories,  however,  and  later  evaluate 
the  status  of  collisional  dissociation  theories.) 

Pritchard  examined  the  rotational  contribution  to  dis- 
sociation and  recombination  in  a series  of  papers  (refs 
70,  71,  72)  and  found  it  to  be  significant.  In  addition, 
he  addressed  the  relation  of  rotational  transitions  to 
observed  activation  energies  and  their  temperature  depen- 
dence. He  concluded  that  inclusion  of  rotation  adequately 
accounted  for  these  variations. 

Recently  Kiefer  (ref  73)  has  taken  a new  and  refresh- 
ing look  at  the  problem  of  dissociation.  He  concludes,  in 
light  of  recent  measurements  of  dissociation  and  relaxation 
rates,  that  there  can  be  at  most  a very  weak  vibrational 
bias  for  dissociation  from  vibrational  states.  The  ladder 
model  is  out.  His  analytic  approach  demanded  consistency 
in  the  magnitude  and  temperature  dependence  of  the  dissoci- 
ation rate  and  induction  or  incubation  time  (the  time  to 
establish  steady  state  dissociation) . He  also  investigated 
the  effect  of  vibrational  depletion.  He  points  out  that 
depletion  is  so  severe  that  any  theory  which  assumes  that 
dissociation  proceeds  only  from  highly  excited  vibrational 
states  is  unacceptable.  This  is  an  extremely  interesting 
comment,  since  only  three  years  ago  (ref  69)  the  same  author 
used  a ladder  model  and  achieved  "quite  satisfactory  agree- 
ment." He  concludes  at  this  time,  however,  that  dissocia- 
tion must  take  place  to  an  appreciable  extent  from  low 
lying  vibrational  states  in  accordance  with  the  early  ideas 
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of  Careri  (ref  41)  and  the  recent  theory  of  Johnston  (ref 
68) . 

Time-dependent  solutions  of  the  vibrational  Master 
Equation  for  hydrogen,  including  dissociation,  have 
recently  been  obtained  by  Kewley  (ref  74)  using  the  stan- 
dard ladder  model  without  rotation.  The  rates  used  were 
of  the  SSH  form  (ref  20)  as  modified  by  Keck  and  Carrier 
(ref  28)  and  fitted  to  the  experimental  relaxation  data. 

The  calculated  dissociation  rate  showed  fair  agreement  in 
magnitude  and  temperature  dependence  with  experimental  data. 
When  V-V  exchanges  changes  were  included  in  the  analysis, 
it  was  noted  that  they  resulted  in  only  a small  enhancement 
of  the  dissociation  rate. 


Contradictions  and  Inconsistencies 

Through  this  review  of  the  collisional  theories  of 
dissociation,  it  has  become  apparent  that  despite  the  intro- 
duction of  new  computational  techniques  and  the  existence 
of  reasonable  transfer  rate  data,  inconsistencies  and  con- 
tradictions exist.  In  order  to  highlight  this  fact  and 
capsulize  the  historical  review  of  collisional  rate  theo- 
ries, a status  table  of  dissociation  theories  is  presented 
in  Table  II.  In  this  table  we  list  the  author  and  year  of 
the  work  cited,  evaluate  the  success  of  the  theory,  classify 
the  approach  and  vibrational  bias,  and  list  the  gases  con- 
sidered. In  addition,  we  point  out  in  those  works  consider- 
ing the  effects  of  V-V  transfer,  whether  or  not  this 
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exchange  mechanism  led  to  an  enhancement  of  the  dissocia- 
tion rate.  Observe  that  calculations  utilizing  the  ladder 
model  have  extended  over  the  past  twenty  years  and  have 
yielded  agreement  ranging  from  good  to  poor.  Fully  time 
dependent  calculations  have  been  limited  in  scope  tnd 
yielded  contradictory  results.  Steady  state  analyses  based 
on  the  ladder  model,  while  having  treated  a wider  variety 
of  gases,  also  have  yielded  contradictory  results.  Investi- 
gations of  rotational  effects  and  empirical  biases  have 
been  limited  to  steady  state  and  equilibrium  analyses,  with 
time  dependent  calculations  concentrating  on  H ^ • 

In  the  present  study,  we  shall  address  many  of  the 
inconsistencies  which  have  arisen.  We  will  attempt  to 
develop  a collisional  model  of  dissociation  that  provides 
reasonably  accurate  values  for  dissociation  and  recombina- 
tion coefficients,  exhibits  an  acceptable  temperature 
dependence,  and  permits  an  evaluation  of  rate  enhancement 
effects  in  pure  gases.  We  will  utilize  both  a detailed  time 
dependent  and  steady  state  analysis  in  a thorough  and  con- 
sistent  evaluation  of  the  kinetics  of  the  dissociation  pro- 
cess and  study  the  following  questions  which  have  been 


2.  What  role  is  played  by  rotational  structure  and 


must  it  be  included  for  accurate  rate  determination? 

3.  Is  it  possible  to  select  a particular  model  based 
on  the  observed  incubation  times? 

4.  Can  the  observed  rate  enhancement  in  pure  gases 
be  attributed  to  complex  exchange  phenomena? 

These  questions  will  be  addressed  by  considering  a variety 
of  gases , not  limiting  ourselves  to  hydrogen  as  have  others , 
and  utilizing  a fully  coupled  time-dependent  solution  of 
the  Master  Equation,  to  validate  any  steady-state  assump- 
tion or  analysis. 


If 
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Chapter  III.  Development  of  the  Master  Equation 

Introduction 

In  this  chapter  we  will  develop  the  Master  Equation 
formalism  which  will  be  applied  to  a collisional  investiga- 
tion of  the  dissociation  reaction.  The  Master  Equation 
approach  to  the  dissociation  problem  is  attractive  because 
it  maintains  the  discrete  nature  of  the  individual  states 
of  the  molecule,  is  exact  in  principle,  and  provides  a 
better  understanding  of  the  course  of  internal  relaxation 
and  its  coupling  to  the  dissociation  reaction.  In  addi- 
tion, this  approach  may  be  easily  transferred  to  the 
related  problem  areas  of  kinetic  modeling  of  lasers. 
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System  and  Assumptions 

The  physical  system  under  consideration  is  a pure  gas 
consisting  of  diatomic  molecules  or  a gas  mixture  of  a 
diatomic  gas  diluted  in  an  inert  gas,  such  as  argon.  We 
will  model  the  dissociation  of  diatomic  gases  as  observed 
under  experimental  conditions  typical  in  shock  tube  studies. 
Under  these  conditions  the  shocked  gas  is  suddenly  raised 
to  a high  temperature  and  maintained  at  constant  volume. 

The  experimental  dissociation  rates  are  inferred  by  using 
a light  absorption  (ref  102)  or  interferometric  technique 
(ref  97)  to  measure  the  density  profiles  behind  the  shock. 

In  our  analysis,  chemical  reactions,  other  than  dissociation 
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and  recombination,  will  be  excluded.  The  energy  of  the 
system  is  contained  in  the  translational,  rotational,  vibra 
tional,  and  electronic  modes  of  the  molecules.  Energy  is 
exchanged  or  transferred  through  the  various  modes  by  means 
of  collisions.  Translational  and  rotational  energy  are 
exchanged  readily,  and  therefore  these  modes  rapidly  equili 
brate  only  requiring  of  the  order  of  100  collisions  (ref 
24).  In  this  analysis,  we  will  be  concerned  primarily  with 
the  exchange  of  vibrational  energy  resulting  from  V-V  and 
V-T  exchanges  and  with  the  development  of  a formalism  for 
treating  the  vibrational  and  dissociative  relaxation  of 
diatomic  gas. 

The  overall  dissociation  process,  although  represent- 
ing one  of  the  simplest  of  all  collisionally  induced  gas 
phase  reactions,  is  extremely  complex.  The  process  is 
governed  by  a set  of  coupled  non-linear  differential  equa- 
tions for  each  of  the  internal  molecular  states  as  well  as 
the  continuum.  Formally  each  state  population  will  be 
governed  by  an  equation 


where  Nj  is  the  concentration  of  molecules  in  state  l,  and 

a, 

the  R's  represent  generalized  rate  coefficients  for  the 
population  or  depopulation  of  state  l . This  set  of  equa- 
tions is  commonly  referred  to  as  the  Master  Equation  and 


is  generally  applicable  (refs  75,  76,  77)  for  determining 
the  rate  law  for  the  dissociation  reaction. 

In  principle,  the  rate  of  molecular  dissociation  could 
be  calculated  by  incorporating  the  detailed  mechanics  of 
the  molecular  collisions  into  the  generalized  rate  coeffi- 
cients and  solving  the  resulting  set  of  equations  for  each 
and  every  state.  Speaking  realistically  though,  this  is 
not  possible  and  probably  not  desirable.  Hydrogen,  simple 
Hj , would  require  close  to  350  states  to  be  tracked  while 
Iodine  would  require  over  50,000  states  to  treat  the 
ground  electronic  state  alone.  The  attempt  in  this  analy- 
sis will  be  to  obtain  a simplified  Master  Equation  while 
retaining  only  those  physically  important  mechanisms 
required  to  describe  and  understand  the  process  of  dissoci- 
ation. To  that  end,  the  following  assumptions  are  intro- 
duced . 

i.  The  translational  and  rotational  degrees  of  free- 
dom are  in  complete  Boltzmann  equilibrium. 

ii.  The  molecules  in  each  vibrational  state  are 
treated  as  separate  species. 

iii.  Only  single  quantum  vibrational  exchanges  occur, 

iv.  The  molecule  is  represented  as  a Morse  oscillator 
rigid  rotator  system. 

~iven  the  experience  of  previous  works,  the  results  of 
2xperimental  relaxation  studies,  and  the  state-of-the-art 
ji  collision  theory,  the  above  assumptions  are  reasonable 


for  the  scope  of  this  research.  As  justification  for  (i) , 
one  may  offer  the  standard  response  of  the  extremely  rapid 
relaxation  of  the  translational  and  rotational  degrees  of 
freedom.  Since  rotation  is  involved  in  the  dissociation 
process,  it  is  possible  that  this  equilibrium  may  be  per- 
turbed. However,  to  effectively  analyze  the  extent  of  any 
departure  from  equilibrium  would  require  rate  data  that  is 
simply  not  available.  Assumption  (iii)  may  be  viewed  as 
quite  severe;  however  the  restriction  to  single  quantum 
transitions  has  been  examined  previously  and  has  been 
shown  to  have  little  effect  on  the  dissociation  rate  (refs 
65,  90) . This  fact  is  consistent  with  the  adiabatic  nature 
of  the  multiple  quantum,  V-T  exchanges.  In  addition,  it 
points  out  that  the  interplay  of  the  dissociation  and 
vibrational  relaxation  processes  results  in  an  increased 
vibrational  depletion  and  a resultant  saturation  of  the 
dissociation  rate.  A Morse  potential  was  selected  because 
of  its  simplicity  and  ease  of  application  to  theoretical 
heavy  particle  collision  probabilities.  The  Morse  poten- 
tial is  able  to  match  very  well  the  spectroscopically 
observed  vibrational  energies  and  retains  the  important 
anaharmonic  features  of  the  molecular  potential.  Here, 
vibrational  energies  a he  ‘measured  from  the  dissociation 
limit  and  within  the  Morse  representation  are  given  by 


(ref  77) 


» 
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for  the  Morse  potential 


V(r)=  Q.  ft-  exp  l«(r-r*)3j 


Dg  is  the  well  depth,  M the  reduced  mass,  rg  the  equil- 
ibrium separation,  u)  the  oscillator  frequency.  Values  of 
the  required  constants  were  obtained  from  Herzberg  (ref 
78)  . 

Transition  Bates 

Let's  now  discuss  the  rates  that  are  to  be  incorpor- 
porated  into  the  Master  Equation.  The  pumping  and  relaxa- 
tion process  controlling  the  state  populations  can  be 
separated  into  two  collisional  classes:  collisions  between 
heavy  particles  and  electron  impact  collisions. 

The  heavy  particle  transfer  rates,  V-V  and  V-T  were 
discussed  in  Chapter  I.  In  summary,  these  processes  can 
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be  represented  by  the  following  reactions  under  the  single 
quantum  assumption: 


N +N  N . . + N Intramolecular  V-V 

v w v+1  w+1 

N +N  N , . + N Intramolecular  V-T 

v w v+1  w 

N +X  — -*•  N , . + X Intermolecular  V-T 

v v+1 

It  is  conventional  to  express  the  heavy  particle 
kinetic  rates  as  the  product  of  the  collisions  frequency, 

Z,  the  applicable  two  populations,  and  the  transition 
probability.  For  the  V-T  processes  we  utilize  the  SSH  form 
(ref  22)  as  modified  by  Keck  and  Carrier  (ref  28)  and  fit 
this  form  to  the  data  of  Millikan  and  White  (ref  78)  and 
Miller  and  Millikan  (ref  80).  The  transition  probability 


for  v+1  + v is  then  written 


and  0 is  the  SSH  group  of  molecular  constants  (ref  22), 
while  PQ  is  obtained  from  the  experimental  vibrational 
relaxation  times. 

We  have  similarly  utilized  the  Keck  and  Carrier  form 
(ref  28)  of  the  SSH  theory  for  the  V-V  probabilities  result 
ing  from  short  range  forces  and  augmented  this  with  a term 
resulting  from  the  long  range  force  contribution  following 
Sharma  and  Brau  (ref  81)  where  appropriate.  The  resulting 
probability  expression  for  v+1 -*■  v and  w ■+  w+1  is 


where  the  short  and  long  range  contribution  are  designated 

y 

by  S and  L.  The  short  range  contribution  is  functionally 
similar  to  the  SSH  V-T  model  while  the  long  range  terms  are 
complex  and  are  matched  by  a simple  algebraic  form.  The 
details  of  both  terms  are  given  by  Herzfeld  and  Litovitz 
(ref  22),  Sharma  (ref  81),  and  Fisher  (ref  82). 

The  electron  impact  collision  rates  are  utilized  only 
in  the  non-equilibrium  portion  of  the  present  work  and  are 
obtained  from  a steady  state  solution  of  the  collisional 
Boltzmann  equation  for  the  electron  energy  distribution. 

We  will  discuss  this  treatment  in  Chapter  VI  and  for  the 
moment  merely  include  the  effective  electron  impact  rates 
in  the  development  of  the  Master  Equation  for  completeness. 


< 
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Master  Equation  Formalism 


Having  discussed  the  physical  system  and  the  rates, 
let  us  proceed  to  write  the  rate  equations  for  the  instan- 
taneous population  of  the  vth  vibrational  state,  Nv« 

v 

* ■ * £ faWA)*  ",  fo,,,  w" 


C3-5) 


-VvcJy  -+  r^A*"  v*c±/r..y* 


Here  the  first  term  accounts  for  electron  impact  collisions 
and  the  last  five  terms  account  for  heavy  particle  colli- 
sions with 


IP**)  = electron  number  density 

= electron  excitation  rate  from  vibrational 
V*  state  v to  k 


V+LV  V,‘ 


V 


k+! 

V 


K Y*I 


= v-v 

and 

transfer 
k ->-k+l 

rate 

for 

v+1  -*■  v 

= V-T 

transfer 

rate 

for 

v+1  V 

= V-T 

transfer 

rate 

for 

v -*•  v+1 
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dyC/sec) 

n/cm3) 
ry  (c*fcec) 
A(c»3) 


= dissociation  rate  out  of  level  v 


= vibrational  population  of  level  v 


= recombination  rate  into  level  v 


= atom  number  density 


For  rare  gas  mixtures,  the  VT  rates  are  summed  over  all 
species  M . In  addition,  it  should  be  pointed  out  that 

S 

microscopic  reversibility  or  detailed  balance  requires 
that 


<:=  VCZ  **/°C/Sfc.rEsr„-£w.J 


where 


i/c. 


'expC-fCP'-Fy)) 

l Qrot.(Ai)  l 


is  the  specific  vibrational  equilibrium  constant  and  with 
the  g's  representing  the  electronic  degeneracy  and  the  Q's 
the  single  particle  partition  functions.  It  should  be 


noted  that 


£ CK] 
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For  further  computations  and  manipulations  it  is  help- 
ful to  introduce  the  following  "lumped"  form  of  the  Master 
Equation 


/iv  ~ ft- I ,V  ~ (fyv+i  +/°Y,V-l) 

- </v  - rv  A* 


where  <5~y  = electronic  pumping  term  for  level  v 


Analytical  solutions  of  these  (v*+l)  coupled  non- 
linear rate  equations  for  an  anharmonic  oscillator  have 
not  been  obtained.  However,  if  we  consider  a highly 
diluted  gas  mixture,  and  consequently  obtain  a linear  set, 
the  formal  solution  is  well  known,  in  terms  of  a series  of 
exponentials  in  the  eigenvalues  (ref  83) . Numerical  pro- 
cedures could  be  utilized  to  determine  the  roots  of  the 
(v*+l)th  order  polynomial.  Rather  than  taking  this 
approach,  we  will  attach  the  general  problem  with  a numeri- 
cal integration  of  the  rate  equations. 
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Numerical  Method 

Given  an  initial  state  population  distribution  and  a 
set  of  suitably  defined  rates,  we  wish  to  perform  a time- 
dependent  integration  of  the  rate  equations.  Since  the 
relaxation  rate  is  governed  by  the  transition  probabilities, 
V-V  and  V-T,  which  range  over  many  orders  of  magnitude,  the 
resulting  time  constants  for  the  coupled  set  differ  sub- 
stantially in  magnitude.  Simply  stated  the  system  of  equa- 
tions is  stiff. 

Gear  (ref  84)  has  summarized  the  general  philosophy  of 
stiffness.  Let  us  here  discuss  stiffness  in  the  context  of 
vibrational  relaxation.  Vibrational  relaxation  kinetics 
include  fast  and  slow  reactions.  The  decay  times  t can  be 
related  to  the  eigenvalues  y of  the  Master  Equation  Jacobian. 
The  fast  reactions  control  numerical  stability  and  step 
size,  while  the  slow  reactions  determine  the  truncation 
error.  In  this  study  the  "multistep  implicit  stiffly 
stable  scheme"  proposed  by  Gear  (ref  85)  was  utilized. 

The  Gear  scheme  utilizes  a predictor-corrector  technique  to 
perform  the  stiff  integration.  The  predicting  technique, 
that  insures  stiff  stability,  is  coupled  to  a correcting 
scheme  that  utilizes  the  Jacobian  to  accomplish  time  step 
adjustment.  The  "Gear  package"  is  an  extremely  user 
oriented,  flexible  code,  and  efficiently  handles  the  inte- 
gration of  the  Master  Equation. 
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Steady  State  Solutions 
of  the  Master  Equation 

A complete  time  dependent  solution  of  the  Master 
Equation  is  complex  and  computationally,  fairly  time  con- 
suming. In  addition,  a direct  integration  while  providing 
results,  often  fails  to  provide  the  insight  that  can  be 
obtained  from  an  analytic  approach.  For  this  reason,  we 
will  also  investigate  steady  state  solutions  to  the  coupled 
rate  equations. 

Any  steady  state  analysis  is  based  on  differences  in 
characteristic  times.  The  relevant  times  here  are  those  of 
vibrational  relaxation  and  dissociation.  If  the  processes 
of  vibrational  relaxation  and  dissociation  are  separable 
with  the  rate  of  vibrational  relaxation  being  greater  than 
the  rate  of  dissociation,  a quasi-steady  state  may  be 
assumed  to  exist  in  the  vibrational  populations.  A more 
general  existence  of  this  steady  state  has  been  documented 
(refs  64,  86,  87)  and  numerous  analyses  based  on  the  assump- 
tion have  evolved  (refs  48,  49,  63,  64,  69). 

We  shall  assume,  and  later  numerically  validate,  the 
existence  of  a steady  state  in  the  following  development  of 
the  vibrational  distribution  and  dissociation  rate.  The 
steady  state  form  to  be  derived  will  then  be  applied  to 
the  situation  of  thermal  dissociation  and  the  highly  non- 


equilibrium situations  obtained  under  electric  discharge 
or  optical  pumping  conditions.  Let  us  proceed  by  writing 
the  generalized  rate  equations,  where  we  will  omit 
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recombination  terms  since  we  are  far  from  equilibrium.  It 
will  be  assumed  that  the  steady  state  achieved  is  charac- 
terized by  the  vibrational  distribution  assuming  a self- 
similar form,  that  is 


/'v  consta/it" 


From  (3-6) 


(3-7) 


^ ■ V*  + Vv  • - H,fly  Hd,  + e, 


0 

m -ZL  ^ 


Defining 


G-ZLfti  +Hidi~  ec} 

S L*o  ^ 


(3-8) 


(3-9) 
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and  substituting  in  (3-8)  yields 


G;  ~ N. 


1 i*' 


(3-10) 


M - 


_ '/-V'  N/-'  + &r‘ 


(3-11) 


r . 


/*/-/ 


in  view  of  (3-7)  ^-/ 

N,  = ±ii  Ht'.  + 

upon  noting  G.=0,  we  obtain 

3 y» 


(3-12) 


N = fyr  % . T- 


(3-13) 


Explicitly  relating  the  V-V  and  V-T  rates  to  the  p's 
after  invoking  detailed  balance  we  obtain  the  simplified 
identification 


//V  * •'  l,tL  Wf?> x *w + ^7 

u-t,w 


(3-14) 


where  «„  = 6 (E,‘E,.,I. 
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Equation  (3-13)  is  an  extremely  simple  steady  state  solu- 
tion and  will  provide  insight  as  to  limiting  forms  of  the 
vibrational  distribution  and  the  sensitivity  and  extent  of 
vibrational  depletion.  For  a specified  set  of  V-V  and  V-T 
rates,  and  any  distribution  of  dissociation  rates  we  can 
use  the  steady  state  form  of  (3-13)  to  calculate  both  the 
vibrational  distribution  and  the  steady  state  dissociation 
rate.  In  addition,  it  is  a stable  representation  of  the 
distributions  amenable  to  rapid  solution  using  iterative 
techniques.  These  results  will  be  presented  in  Chapter  V. 

Here,  let  us  first  establish  the  various  limiting 
forms  for  the  distribution.  With  no  dissociation  or  pump- 
ing,  Gj=0 . Then  Nj=Nj_^e  j;  the  distribution  is  Boltzmann 
as  anticipated.  The  Treanor  distribution  (ref  29) 


Hj  - & exp  ( 4 j ) 


or  = exp  (<f>- 

t/f-t 

where  <J>  is  the  vibrational  temperature  parameter,  must 
also  be  obtained  when  only  V-V  terms  are  retained.  Sub- 
stituting into  (3-13)  we  obtain 


(3-15) 
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which  reduces  to 


Thus  again  a consistent  form  is  yielded. 

Next  let's  look  at  a situation  typical  of  many  shock 
tube  experiments  where  a diatomic  gas  is  highly  diluted 
in  an  inert  gas.  Under  these  circumstances  we  may  neglect 
the  V-V  and  e terms  and  obtain  from  (3-13) 


(3-16) 


From  (3-16)  we  observe  that  the  dissociation  process  causes 
the  vibrational  populations  to  deviate  from  the  equilibrium 
distribution.  This  deviation  results  in  a reduction  in  the 
dissociation  rate  below  the  value  that  would  be  obtained  if 
no  deviation  occurred.  The  extent  and  character  of  this 
departure  from  equilibrium  is  dictated  by  the  behavior  of 
the  summation  term  in  (4-8) . This  term  is  controlled  by 
the  transfer  rate,  V-T,  and  the  magnitude  and  bias  of  the 
dissociation  rates,  dv-  Note  that  for  the  special  case, 
where  dy  is  a constant,  there  is  no  depletion.  For  the 
V-T  rates  and  all  biases  considered  in  this  study, 
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dissociation  monotonically  decreases  the  population  of  the 
vibrational  states  below  their  equilibrium  value. 

Let  us  now  develop  an  alternate  steady  state  solution, 
based  on  the  analysis  of  Nikitin  (ref  50) . This  solution 
yields  an  analytic  form  for  the  dissociation  rate  and  does 
not  require  iteration.  However,  no  information  regarding 
the  vibrational  distribution  is  obtained.  Consider  a 
diatomic  gas  diluted  in  an  inert  gas  such  as  argon,  with 
dissociation  proceeding  only  from  the  last  bound  level  to 
an  unbound  pseudo- level.  From  (3-5)  , upon  neglecting  V-V 
terms , we  obtain 


(3-17) 


with 


= I « 

Writing  out  this  first  few  rate  equations  explicitly  we 


obtain 


fa  /v,  %,-z- .z) 

fa  rt, ^z-fi(zn,,*z ‘ 

fa  z + w} 


for  k“0,l,...v*  and  with 


ft,m° 


(3-18) 
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As  mentioned  earlier,  the  general  solution  can  be  repre- 
sented* in  terms  of  the  eigenvalues,  ym,  of  the  resulting 
matrix  B 


resulting  in  a general  solution  of  the  form 


(3-19) 


where  the  eigenvalues  are  ordered  in  increasing  order  of 

V • 
m 

Within  the  framework  of  the  standard  ladder  model,  dis- 
sociation will  proceed  from  the  last  bound  level  to  an 
unbound  pseudo-level,  with  a rate  given  by  ?v*+1  v*e  v*+l 

— /C 

= P e v*+l.  If  P ~0  which  implies  no  dissociation,  the 
matrix  B yields  the  eigenvalues,  y°'  corresponding  to  the 
thermal  relaxation  of  the  system  with  y°=0.  If  P^O , all 
values  of  ym  will  be  displayed  in  the  same  direction  rela- 
tive to  y°.  Therefore  y >y°  when  dissociation  is  present. 

If  y^,  the  distribution  will  reach  a pseudo-steady-state 

and  exponentially  decay,  reflecting  the  thermal  dissociation 

....  *■  ' »■ 

♦Provided  B is  non-defective,  ’’hich  we  assume. 
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or  decomposition  of  the  molecule.  Assuming  that  this 
inequality  holds,  the  smallest  eigenvalue,  yQ,  correspond- 
ing to  the  dissociation  constant,  can  be  found  by  expanding 
the  determinant  in  the  power  series  of  y while  retaining 
only  first  order  terms.  For  a two  level  system  the  rate 
equations  are : 


- Z [k £ 

(3-20) 

4 • z{ 


where 


The  determinant  is  given  by. 


~(?o 


- o 


(3-21) 


where  * ex/o  [En'En.i}) 

Setting  the  determinant  equal  to  zero,  we  obtain 

/?o 4 /0(%o  +%,/  •*)"  */  Ko  ” O 


Thus  A = 

x %o04+)+% 


±Zi!L  = U-+-L.  * ' 7 


22) 
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Similarly  for  a three  level  system. 


//-  + Cl  +-L  V—  X 

'jz  £ ?a  1 •*/•£*£  V •*/  _/ 


(3-23) 


Since^Vr^Vk  = exp  (En“En-k-l*  1 this  form  Can 


be  simplified  to  yield, 


(3-24) 


/(.  -X  Te 

-/■  j— 

rv  [f  +e  J p,,° 

• £ & It  ’*]} 


We  can  now  generalize  to  a system  with  (v*+l)  levels  to 
obtain. 


(3-25) 


This  differs  from  the  form  obtained  by  Nikitin  (ref  50) 
He  obtains 


* 4‘o  ( L *=1  J 


(3-26) 


Equation  (3-25)  is  a general  form  which  will  be 
utilized  later  in  a steady-state  analysis.  As  pointed  out 
earlier,  in  contrast  to  (3-13)  , it  yields  no  information 
on  the  vibrational  distribution  and  is  limited  to  the  case 
where  dissociation  occurs  only  from  the  last  bound  state. 
However,  for  a given  set  of  V-T  rates  and  a specified  dis- 
sociation rate,  it  does  provide  the  dissociation  constant 
and  does  not  require  iteration.  For  the  moment,  in  order 
to  attempt  to  establish  an  upper  bound  on  the  dissociation 
constant  and  examine  depletion  effects,  we  will  assume  all 
P.  . , =1.  This  can  later  be  refined  by  replacing  the 
P.  -i_i  by  effective  average,  <P.  . ,>.  Let  us  now  assume 
that  terms  wi th  large  i dominate  the  summation  due  to  the 


term  exp  (BE^) . Then  the  summation  over  s rapidly 


approaches  the  standard  vibrational  partition  function  and 
can  be  removed  from  the  summation.  Thus 


x'-.'v 


(3-27) 


Recalling  that  for  a Morse  potential  with 


(3-28) 


E^  can  be  represented  as  (Ref  61) , 


5-- 


with  the  highest  bound  level,  v*,  given  by 


y*  — c/~  A ~ & 


(3-29) 


with  6 a positive  number  less  than  1.  Defining  the  qth 
vibrational  level  below  v*  and  relating  it  to  i,  we  obtain 


or 


then 


i+;/z  = c/-cf~  J 


(3-30) 


Substituting  in  (3-27)  we  obtain 


Q,ti « 

j=o 


(3-31) 


Where  upon  approximating  the  summation  by  an  integral  we 
obtain 
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- Z tuJ  e~^D° 
Qv/6.  [b  Q4tJ//x- 

So  that 


iu/e 


■A D. 


J r-ZA ^ 


(3-32) 


At  this  point  we  should  note  that,  under  the  assumption  of 
P-  • ,=1;  the  rate  of  dissociation  from  the  last  bound  level 

1 f 1“  J. 

assuming  a Boltzmann  distribution  for  the  vibrational  levels 
is  given  by 


Z e 


Qv/6. 


(3-33) 


Thus  we  may  identify  the  term,  [ttD  kTK  / ("hw)  , as  a depletion 
factor  resulting  from  the  coupling  of  vibrational  transfer 
and  dissociation.  This  term  reduces  the  dissociation  rate 
below  its  equilibrium  value  and  introduces  a stronger 
negative  temperature  dependence  to  the  rate  which  results 
in  a lowering  of  the  activation  energy. 
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Conclusion 

In  this  chapter,  we  have  developed  the  formalism, 
theory,  and  numerical  approaches  required  for  time-dependent 
and  steady  state  solutions  of  the  Master  Equation.  Recog- 
nizing the  stiff  nature  of  the  set  of  simultaneous  differ- 
ential equations,  we  have  selected  the  Gear  (ref  85)  inte- 
gration technique  for  their  solution.  This  time-dependent 
solution  will  be  used  to  examine  the  temporal  behavior  of 
the  dissociation  process,  assess  the  degree  of  vibrational- 
dissociative  coupling,  and  validate  the  subsequent  use  of 
both  steady  state  analyses. 


Chapter  IV.  Dissociation  Models 


\ 


i 


General 

Having  established  the  Master  Equation  formalism,  and 
developed  the  computational  approaches  to  be  used  in  this 
study,  we  are  now  prepared  to  address  the  dissociation  pro- 
cess itself.  Within  the  framework  of  the  four  assumptions 
given  in  Chapter  III;  Morse  oscillator,  single  quantum 
transitions,  rotational  equilibrium,  and  distinct  vibra- 
tional states,  we  will  now  present  three  models  of  the  dis- 
sociation reaction.  All  three  are  based  upon  the  ladder 
climbing  model,  the  differences  lying  in  how  rotational 
effects  are  treated.  Through  a consistent  examination  of 
these  models,  we  will  attempt  to  resolve  the  inconsisten- 
cies and  contradictions  that  have  arisen  regarding  the 
vibrational  bias  of  the  dissociation  reaction,  the  role  of 
rotation,  and  the  source  of  rate  enhancement  in  pure  gases. 


Standard  Ladder  Model 

The  first  model  is  designated  as  the  standard  ladder 
model.  It  has  been  used  extensively  and  recently  (refs  65, 
88)  in  time  dependant  and  steady  state  analyses  of  dissoci- 
ation. Application  of  this  model  has  yielded  inconsistent 
results,  as  was  brought  out  in  the  review  of  :ollisional 
theories  of  dissociation.  The  model  can  be  understood  by 
referring  to  Fig.  4-1  which  displays  a typical  potential 


j 
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INTERNUCLEAR  SEPARATION 

Fig  4-1  Standard  Ladder  Model 
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energy  diagram  using  the  Morse  representation.  Super- 
imposed is  the  vibrational  energy  structure  with  the  states 
being  identified  by  the  vibrational  quantum  number  v.  Dg 
is  the  potential  well  depth  and  Dq  the  dissociation  energy. 
During  the  relaxation  process,  the  molecules  are  excited  to 
higher  and  higher  vibrational  states  and  dissociate  by 
making  a transition  to  the  pseudo-level,  designated  by  the 
dashed  line.  Rotational  structure  is  neglected.  This 
model  was  used  by  Montroll  and  Shuler  (ref  56),  Dove  and 
Jones  (ref  65) , and  Kewley  (ref  74)  for  studies  of  the  dis- 
sociation of  hydrogen  diluted  in  argon  and  by  Wengle  (ref 
12)  for  studies  of  hydrogen  and  oxygen  in  an  argon  heat 
bath. 

Rotational  Structure 

As  an  introduction  to  the  models  which  include  the 
effects  of  rotation,  we  discuss  in  this  section  rotational 
energy  and  structure,  and  the  influence  of  rotation  on  the 
dissociation  process.  Fig.  4-2  shows  how  the  potential 
energy  curves  of  a diatomic  molecule  depend  on  the  rota- 
tional quantum  number  J.  The  Morse  potential,  V , is  now 
augmented  by  the  centrifugal  potential,  V'.  This  term 
results  in  a reduction  in  the  potential  well  depth  and 
shift  in  the  potential  minimum  to  a larger  internuclear 
separation.  A rotational  or  centrifugal  barrier  arises  and 
is  usually  of  sufficient  height  to  result  in  the  formation 
of  one  or  more  quasi-bound  states  which  are  metastable  with 
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respect  to  the  normal  dissociation  limit  for  J=0.  These 
states  can  predissociate  via  a penetration  of  the  barrier, 
but  in  most  cases  this  is  of  only  secondary  importance 
(refs  70,  71),  due  to  the  width  of  the  barrier.  Because  of 
this  barrier,  the  effective  dissociation  energy  as  mea- 
sured relative  to  the  v=0,  J=0  state  is  increased  as  a 
result  of  rotation.  The  molecular  energy  levels  are  also 
modified  from  those  given  previously  in  (3-2)  and  are  now 
given  by  (ref  77) 


While  examining  rotational  effects,  some  interesting 
consequences  previously  cited  (refs  62,  70,  90,  91)  should 
be  mentioned.  Rotational  effects  will  alter  the  molecular 
partition  functions  and  provide  new  pathways  to  the  dis- 
sociating pseudo-state.  The  graphical  representation  in 
Fig.  4-3,  proposed  by  Bauer  and  Widom  (ref  62),  gives  a 
schematic  representation  of  these  new  possible  vibrational- 
rotational  energy  transfer  processes.  At  each  point  in  the 
plot  the  total  energy  is  given  by  the  sum  of  the  abscissa, 
(rotation) , and  ordinate,  (vibration)  values.  The  heavy 
curve  denotes  the  dissociation  limit,  and  the  dashed  lines 


represent  constant  total  energy.  Note  again  that  a mole- 
cule must  accumulate  more  total  energy  if  it  is  to  dissoci- 
ate from  a high  J state,  which  would  imply  a reduction  in 
the  dissociation  rate  relative  to  the  J=0  state.  This 
reduction  is  tempered,  on  the  other  hand,  by  the  fact  that 
the  statistical  weight  of  this  state  is  greater,  increasing 
linearly  with  J.  Energy  exchanges  are  easily  represented 
on  this  diagram  as  vectors.  Dissociative  transitions 
involving  changes  only  in  vibrational  energy.  A,  will  be 
vertical.  Pure  rotational  transitions,  B,  will  be  hori- 
zontal, while  vibrational-rotational  transitions  occurring 
at  constant  energy,  C,  will  be  inclined  at  45°.  Transitions 
of  the  type  A and  C will  be  considered  in  this  study,  as 
they  appear  to  be  the  most  efficient  from  a required  energy 
standpoint  and  can  be  handled  within  reasonable  computer 
storage  requirements.  An  A type  transition  requires  that 
the  molecule  surmount  the  rotational  barrier,  while  a C 
type  transition  involves  the  simultaneous  reduction  in  J 
and  increase  in  v.  The  C process  is  not  the  same  as  tunnel- 
ing through  the  rotational  barrier.  It  does  have  a similar 
effect,  in  that  the  barrier  is  lowered  below  the  final 
vibrational  state  and  dissociation  occurs.  The  importance 
of  such  collisions,  especially  for  heavy  molecules,  has 
been  intuitively-argued  (refs  62,  90)  and  recent  trajectory 
calculations  (ref  92)  have  confirmed  their  importance.  We 
will  discuss  equal  energy  transitions  in  greater  deal  in  a 
following  section. 
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Let's  recap  the  implications  of  rotational  effects 


before  proceeding  to  introduce  the  two  rotational- 
vibrational  dissociation  models.  As  a result  of  the  vari- 
ation of  the  rotational  barrier  with  J and  the  variation  of 
the  rotational  population  with  gas  temperature,  the  effec- 
tive molecular  potential  is  temperature  dependent.  In 
addition,  including  rotational  structure  introduces  new 
dissociative  states  and  transitions  which  will  be  found  to 
enhance  the  dissociation  process  and  distribute  the  reaction 
among  the  vibrational  states. 

Rotational  Models 

In  this  section,  two  rotational-vibrational  models  of 
the  dissociation  process  will  be  introduced.  Both  models 
will  include  the  rotational  effects  discussed  in  the  pre- 
vious section,  but  will  differ  with  respect  to  the  form  of 
the  assumed  vibrational-rotational  coupling.  The  form  of 
the  assumed  coupling  is  directly  related  to  the  question 
of  vibrational  bias  and  consequently  vibrational  depletion. 
By  considering  these  models,  in  addition  to  the  standard 
ladder  model,  we  can  evaluate  the  sensitivity  of  the  dis- 
sociation process  to  the  distribution  of  the  dissociation 
rates  among  the  vibrational  levels.  As  in  the  standard 
ladder  model,  dissociation  will  be  assumed  to  occur  upon 
making  a transition  to  an  unbound  pseudo- level. 

Referring  to  Fig.  4-2,  we  introduce  the  first  rota- 
tional model  which  will  be  designated  the  uncoupled  model. 
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In  this  model,  each  J value  wilj  represent  a distinct 
molecular  species  and  vibrational  ladder.  This  model  will 
be  analyzed  under  the  standard  ladder  model  procedures 
using  a time-dependent,  steady-state  or  analytic  approach. 

In  contrast  to  the  uncoupled  model,  the  second  rota- 
tional model  assumes  that  rotational  states  are  closely 
linked  or  strongly  coupled  to  a given  vibrational  state. 

This  model  is  called  the  coupled  model  and  is  represented 
in  Fig.  4-4.  Here  we  plot  the  J=0  vibrational  potential 
and  superimpose  the  J=0  vibrational  energy  level  structure. 
To  the  right,  the  rotational  energy  level  structure,  E(v,J) 
for  v=0  and  v=6  is  also  presented.  Under  the  strongly 
coupled  assumption,  dissociation  now  proceeds  from  the 
last  bound  rotational  state  of  every  vibrational  level. 

Dissociating  States 

The  first  step  in  the  analysis  of  collisionally 
induced  dissociation  is  the  identification  of  the  dissoci- 
ating states.  These  are  the  last  bound  or  quasi-bound 
vibrational  states,  v* ( J) , from  which  a dissociative  tran- 
sition (to  the  unbound  pseudo-level) , takes  place.  The 
states  were  identified,  within  the  adopted  Morse  formalism, 

by  a numerical  calculation  of  the  bound  vibrational  states 
* 

using  (4-1) . The  standard  ladder  model  omits  rotation  and 

•There  exists  some  uncertainty  associated  with  the 
exact  specification  of  the  dissociating  states.  This  uncer- 
♦ainty  does  not  critically  o feet  the  present  analysis,  but 
i«  liscussed  further  in  Appendix  B. 
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is  analyzed  as  the  J=0  case.  The  rotational  models,  how- 
ever, consider  all  values  of  J for  which  there  exists  a 
bound  or  quasi-bound  vibrational  state. 

Figs.  4-5  to  4-10  show  the  calculated  potential  curves 
for  H2 » ^ , O2 » CO,  Br2 / and  parameterized  in  J.  The 
plots  are  linear  with  the  ordinate  given  by  the  normalized 

energy,  E/Dq.  The  abscissa  is  the  normalized  internuclear 
r-r 

separation  - . The  family  of  curves  represent  values  of 
e 

the  rotational  quantum  number,  J.  The  value  of  J and  the 
last  bound  vibrational  level  v*  are  presented  in  the  plot 
legend.  The  influence  of  the  centrifugal  potential  is 
striking,  resulting  in  a reduction  in  the  depth  of  the 
potential  well  and  a shift  in  the  potential  minimum  to 
higher  energies.  Of  more  significance  is  development  of 
the  rotational  barrier,  which  is  of  sufficient  height  to 
form  one  or  more  quasi-bound  states.  These  states  have 
energy  greater  than  the  J=0  dissociation  energy,  D^. 

Having  identified  the  dissociating  states,  it  is  now 
necessary  to  determine  their  relative  importance  based  on 
their  population.  In  equilibrium  the  population  of  a spe- 
cific v,  J state  is  given  by 

/j&J'J  _ C2  J+/)exp  (rtfCCvyT)) 

Ntot  Qtot 

c 

where  N.  . is  the  total  number  of  molecules  and  Q.  . the 
tot  tot 

total  partition  function.  The  population  of  the  last  bound 
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SEPARATION 

Fig  4-5  Effective  Potential  Energy  Curves  for  the  Hydrogen  Molecule 
(Curves  correspond  to  values  of  J»0  to  Jmax.  a*-  AJ  increments.  Last  bound 
vibrational  state  for  a given  J is  v*,) 


0.50 
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SEPARATION 

Fig  4-10  Effective  Potential  Energy  Curves  for  the  Iodine  Molecule 
(Curves  correspond  to  values  of  J=0  to  J at  AJ  increments.  Last  bound 
vibrational  state  for  a given  J is  v*.) 


vibrational  state,  v* (J) , for  a given  value  of  J will  be 


given  by 


/V Cv*(J‘),'7)  - (2J+1)  exp *>J>) 
Mfot  Qttf 


(4-3) 


Assuming  that 


£(/&>,*)  = 0*S&) ) 


where  6 represents  the  variation  in  the  rotational  barrier 
height , then 


^tot  Qtot 


(4-4) 


For  all  gases  examined 


fa  <j,«7 


where  a1  is  a constant  depending  on  the  gas  considered. 
Incorporating  this  dependence  into  (4-4)  yields 


/y/H _ ( zJ+Oexp  &A>)ex£> 

A\ tot  Qtot 


(4-5) 
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which  gives  the  population  of  the  last  bound  vibrational 
state  as  a function  of  J and  temperature.  It  should  be 
noted  at  this  point  that  N(v*(J),  J)  is  distributed  in  a 
manner  similar  to  the  normal  rotational  distribution 
(ref  78) 


N(j)  frJ>/)exp  C-/0tc BJfr+i))  (4_6) 

^tot  Qrot 

The  peak  value  of  the  N(v*,J)  distribution  varies  approxi- 

ic 

mately  as  (T)  as  does  the  peak  of  N(J)of  (4-6).  Typical 
distributions  of  N (v* , J) /N (v* , J)  are  given  in  Fig. 

4-lla,b  and  are  designated  by  the  circled  points  with  the 
scale  shown  on  the  right.  Superimposed  on  these  curves  are 
plots  of  v* (J)  as  a function  of  J.  One  can  see  from  curve 
(11a),  that  the  maximum  population  occurs  for  j=260,  with 
a corresponding  v*=58,  while  from  (11b),  the  maximum  in 
N(v*,J)  has  shifted  down  to  J=200  and  v*=68.  For  the 
coupled  rotational  model,  it  is  necessary  to  obtain  the 
ratio  of  the  population  in  the  lowest  bound  rotational 
state,  J,  for  a given  vibrational  state  to  the  total  popula- 
tion in  that  vibrational  state.  This  ratio  is  obtained 
from  4-3  and  is  given  by 

21  = Hfv(sh  £1 (4_7) 

'%  C*Stf)  exp  f/SSfaj)) 
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This  factor  is  then  incorporated  into  the  coupled  rota- 
tional model  to  determine  the  weighting  of  the  dissociation 


rate  dv,  given  by 


NCv,*)- 


(4-8) 


with  R , representing  the  sum  of  the  V-T  and  V-V  rates 

V / V T _L 

leading  to  a transition  from  state  v to  v+1.  It  should  be 
noted  that  when  rotational  structure  is  neglected  this 
representation  degenerates  to  the  standard  ladder  model 
form  with 


By  introducing  rotational  structure  we  have  increased 
the  number  of  states  in  molecular  structure  being  considered 
and  also  increased  the  number  of  dissociating  channels. 

Let  us  now  examine  if  this  rotational  structure  enhances 
the  possibility  of  dissociation;  and  if  so,  how  does  the 
enhancement  vary  with  temperature  and  gas  considered.  For 
this  purpose,  we  will  assume  here  that  all  populations  have 
their  equilibrium  distribution  and  compare  the  ratio,  t\,  of 
the  total  dissociating  state  population  to  the  total  popula- 
tion. For  the  standard  ladder  model  the  ratio  is  given 
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7?  = exo  (-/SECvI) 

QVj6 


While  upon  including  rotation,  we  obtain 


(4-9) 


v*o 


(2?+l)**p(-/SE(v,  J*)) 


v?  JL 

/ / (z^O  cxpfaE’CVjJ’)) 


(4-10) 


V«o  iT*o 


In  Table  IV,  the  values  of  hladder,  nrot  and 

W.  bladder  are  given  for  H2 ' °2'  N2'  C0'  Br2  and  I2* 

(The  meaning  and  implications  of  ti  , also  listed 

in  the  table,  will  be  considered  later.)  Observe  that 
rotation  not  only  increases  the  number  of  dissociation 
channels,  but  increases  the  critical  state  populations 
(i.e.,  hrot  /hladder  is  always  greater  than  1.).  This 
suggests  a resulting  rate  enhancement  for  dissociation. 

The  ratio,  hrot  /^ladder'  typicaHy  a factor  of  five, 
except  for  where  the  enhancement  is  approximately  20. 

Of  equal  importance  with  the  magnitude  of  the  enhancement, 
is  the  temperature  dependence.  The  ratio  nrot  ladder 
decreases  with  increasing  temperature  with  the  power  of  the 
temperature  dependence  ranging  from  ^ = 1.2  in  to 
<J>  = 0.33  in  l£,  Br2 , where  <f>  is  defined  by 
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TABLE  IV 


ft  1 

I | 

R I 


COMPARISON  OF 

DISSOCIATING 

STATES 

Gas 

Temperature 

^ladder 

Vt.  n: 

rot.  ^ladder 

'"’c/^ladder 

»2 

5000 

1.45-5 

2.1-4 

14.6 

11.7 

4000 

1.09-6 

2.09-5 

19.1 

18.5 

3500 

1.68-7 

3.85-6 

22.6 

22.8 

3000 

1.38-8 

3.74-7 

27.5 

27.6 

2500 

3.95-10 

1.37-8 

35.1 

36.2 

°2 

10000 

5.01-4 

2.32-3 

4.64 

6.50 

5000 

2.67-6 

1.75-5 

6.56 

14.10 

3000 

1.54-9 

1.23-8 

8.01 

25.76 

N2 

14000 

5.54-5 

3.27-4 

5.91 

9.2 

10000 

2.87-6 

2.02-5 

7.03 

13.5 

8000 

1.98-7 

1.53-6 

7.77 

17.4 

CO 

15000 

3.29-5 

2.19-4 

6.67 

9.7 

10000 

6.61-7 

5.55-6 

8.40 

15.4 

8000 

3.22-8 

3.00-7 

9.33 

19.9 

Br2 

2000 

2.15-6 

9.94-6 

4.62 

12.6 

1300 

6.67-7 

3.19-6 

4.78 

14.1 

1200 

1.67-9 

9.02-9 

5.38 

22.1 

X2 

2000 

1.76-5 

6.86-5 

3.9 

9.4 

1500 

1.18-6 

5.04-6 

4.26 

12.8 

1250 

1.29-7 

5.82-7 

4.49 

15.7 
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These  results  are  in  good  agreement  with  the  analysis  of 
Pritchard  (ref  72)  for  all  gases  except  H2-  Pritchard  used 
an  exact  potential  curve  for  H2 , while  using  a Rydberg- 
Klein-Rees  form  for  the  other  gases.  The  difference  for  H2 
may  reflect  an  inadequacy  of  the  Morse  representation. 

Using  Pritchard's  values  or  the  Morse  values  for  6 and 
v* (J)  resulted  in  minor  differences  in  the  ultimate  dissoci- 
ation rate  calculated. 

The  negative  temperature  dependence  of  the  ratio, 
nrot  /rl^a(j<jer'  suggests  that  including  rotational  effects 
will  not  only  increase  the  dissociation  rate,  but  also  lower 
the  activation  energy.  Recall  from  Chapter  I,  that  the 
observed  activation  energy  is  less  than  the  dissociation 
energy.  The  dissociation  coefficient,  assuming  an  equili- 
brium vibrational  distribution  and  gas  kinetic  losses  out 
of  the  critical  dissociating  states,  CDS,  is  given  by 


t 


Z(#m  CDs)/rfTot 


For  the  standard  ladder  model. 


J>  - Z ex/?  v*) 


(4-11) 
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While  for  rotation. 


^ *>  T\Mr 


(4-12) 


with  being  a constant.  Thus  rotation  enhances  dissocia- 
tion and  would  appear  to  yield  "m's",  as  given  in  (1-4) 
and  shown  in  Table  I,  less  than  those  obtained  in  the 

standard  ladder  model,  where  m . =m.  ..  -4>.  While  this 

rot.  ladder 

is  promising,  we  must  postpone  final  judgement,  since  up 
to  this  point  we  have  not  allowed  for  the  non-equilibrium 
effect  of  vibrational  depletion.  The  influence  of  vibra- 
tional bias  and  vibrational  depletion  will  be  addressed  in 
the  next  chapter. 


An  Alternate  Path  to  Dissociation 

As  initially  suggested  by  Bauer  and  Tsang  (ref  62) , 
revived  by  Pritchard  (refs  70,  71,  72),  and  confirmed  by 
recent  trajectory  calculations  (ref  92) , there  is  an 
alternate  path  for  collisional  dissociation.  This  path  is 
an  equal  energy  transition  involving  a quasi-bound  vibra- 
tional rotational,  V-R,  level  and  was  identified  earlier  as 
process  C in  Fig.  4-3.  While  addressing  this  process 
and  its  possible  influence  on  the  dissociation  rate,  let's 
classify  the  various  V-R  levels  and  identify  the  pathways 
to  a dissociating  state. 
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Fig.  4-12  portrays  a typical  set  of  effective  poten- 
tials for  a diatomic  molecule.  Fig.  4-12a  represents  the 
ground  rotational  state  potential  and  shows  a transition  of 
the  A type.  Fig  4-3,  from  a bound  vibrational  state  to  an 
unbound  pseudo-level.  This  is  the  typical  transition  util- 
ized in  the  majority  of  ladder  climbing  models.  Fig.  4-12b 
illustrates  the  change  in  the  effective  potential  due  to 
rotational  effects  for  various  rotational  quantum  numbers. 

In  this  picture,  the  A type  transition  can  be  identified 
again  as  a collisional  excitation  to  an  unbound  pseudo- 
level where  the  molecule  "climbs  over"  the  rotational  bar- 
rier. In  this  case,  however,  the  excitation  proceeds  from 
a quasi-bound  state.  In  addition,  the  molecule  can  now 
penetrate  the  rotational  barrier,  process  P,  through  tunnel- 
ing (predissociation) , but  this  has  been  established  (refs 
70-71)  to  have  a low  probability  due  to  the  extent  of  the 
barrier.  There  is  yet  another  path  to  the  dissociating 
states,  identified  by  process  C (see  Figs.  4-12b  and  4-3). 

In  this  process  the  molecule  does  not  surmount  the  rota- 
tional barrier  nor  does  it  tunnel  through  it.  Process  C 
is  an  equal  energy  transition  from  a quasi-bound  (v' , J') 
state  to  an  unbound  (v,  J)  pseudo-level  of  near  equal 
energy,  where  v'  < v and  J'  > J.  It  illustrates  the  mech- 
anism whereby  a molecule  trapped  in  the  quasi-bound  state, 
undergoes  a collision  which  removes  some  angular  momentum 
and  simultaneously  increases  its  state  of  vibrational 
excitation.  If  this  process  is  very  efficient,  it  will 
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proceed  in  a cascaded  fashion  to  successively  reduce  the 
rotational  barrier  height  ultimately  leading  to  the  effec- 
tive disappearance  of  the  barrier  as  depicted  in  Fig. 

4-12c.  An  upper  limit  to  the  dissociation  rate  when  this 
mechanism  is  effective  can  be  obtained  by  now  redefining 
the  critical  state  as  the  last  bound  vibrational  rotational 
state.  This  is  in  contrast  with  our  previous  designation, 
when  discussing  the  A type  mechanism,  of  the  critical  state 
being  the  last  quasi-bound  vibrational  rotational  state. 
Thus  equal  energy  transitions  essentially  remove  the  rota- 
tional barrier  and  lead  to  a new  distribution  of  critical 
dissociating  states,  CDS. 

Let's  examine  and  contrast  the  distribution  of  CDS 
under  conditions  when  mechanism  A or  mechanism  C are  in 
effect.  Fig.  4-13a  presents  the  A type,  "over  the  top," 
picture  for  a quasi-bound  state.  The  distribution  of  CDS's 
is  given  by  equation  (4-3,  4-4)  and  displayed  in  Fig.  4-13b 
The  distribution  is  depressed  at  large  values  of  J due  to 
the  dominance  of  the  rotational  barrier  and  departs  from 
its  initial  linear  behavior  in  J.  Should  equal  energy 
transitions  be  highly  effective,  we  obtain  the  C type 
picture  presented  in  Fig.  4-14.  The  effective  potential 
for  a given  value  of  J is  presented  in  Fig.  4-14a,b  with 
the  unbound  pseudo- level  given  by  a wavy  line,  the  quasi- 
bound states  by  the  dashed  lines,  and  the  bound  states  by 
solid  lines.  Fig.  4-14a  gives  the  effective  potential  and 
vibrational  states  for  an  intermediate  value  of  J,  and 
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identifies  level  v*  as  the  critical  dissociation  state  from 


which  collisional  energization  must  proceed.  As  J increases 
v*  will  range  from  v*  , given  by  J=0  down  to  v*=0  for 

IUclX 

J=J  . J represents  the  largest  value  of  J for  which  the 
v=0  state  is  bound.  For  very  large  values  of  J,  J>Jc, 

Fig.  lib,  no  bound  states  exist.  In  the  equal  energy  trans- 
ition picture  of  type  C,  the  distribution  of  CDS's  is  given 
approximately  by  equation  (4-4)  with  6=0.  The  calculated 
distribution  is  displayed  in  Fig.  4-14c.  Note  that  the  C 
distribution  is  modified  from  that  of  the  A picture,  and 
exhibits  an  enhanced  population  in  CDS's  especially  at  large 
J values.  This  reflects  the  effective  lowering  of  the 
potential  barrier  achieved  through  the  equal  energy  transi- 

l 

tion. 


Table  IV  presents  the  ratio  of  hc/nlad(jer  where 


and  J is  the  maximum  value  of  J for  a given  v that  results 
in  a bound  vibrational  state.  Observe  that  nc/n2adder  > 
^rot  ^ladder  or  that  t^le  critical  state  population  is 
enhanced  in  the  C picture  over  that  obtained  in  the  A 
picture  and  hence  over  the  standard  ladder  model.  The 
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previously  cited  for  the  A process,  and  indeed  even 
exhibits  a stronger  temperature  dependence. 

The  effects  of  changes  in  the  CDS  distribution  and  its 
temperature  dependence  will  be  evaluated  in  the  next 
chapter  where  detailed  results  including  the  non-equilibrium 
effects  of  vibrational  depletion  will  be  presented  and  con- 
trasted for  the  three  dissociation  models 

1.  standard  ladder 

2.  coupled  rotation 

3.  uncoupled  rotation 
which  we  have  developed  here. 
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Chapter  V.  Calculated  Dissociation  and 
Recombination  Rates 

General 

Having  established  the  Master  Equation  formalism  and 
developed  the  kinetic  models  of  dissociation,  we  now  pre- 
sent the  results  obtained  for  the  dissociation  and  recom- 
bination rate  coefficients.  We  begin  with  an  examination 
of  the  temporal  evolution  of  the  dissociation  process  and 
the  behavior  of  the  incubation  times.  Following  the  time- 
dependent  results,  dissociation  rates  in  diluted  mixtures 
are  presented.  Next,  rate  enhancement  in  pure  gases  will 
be  addressed  by  considering  the  aspects  of  changes  in  col- 
lisional  efficiency  due  to  mass  variation  and  V-V  enhance- 
ment. Finally,  the  calculated  recombination  coefficients 
will  be  presented  and  compared  with  experiment. 

Matrix  of  Computational  Capabilities 

As  pointed  out  in  Chapter  III,  the  dissociation  pro- 
cess, one  of  the  simplest  of  all  collisionally  induced  gas 
phase  reactions,  is  nevertheless  very  complex.  While  a 
detailed  time  dependent  solution  of  the  full  set  of  kinetic 
equations  for  each  electronic,  and  vibrational-rotational 
state  is  possible  in  principle,  it  represents  an  undesir- 
able approach  due  to  the  large  number  of  equations  involved 
and  the  lack  of  detailed  exchange  rate  data  for  each  of  the 
states.  In  this  study,  therefore,  we  will  parameterize  the 
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the  dissociation  process  over  a wide  range  of  temperatures 
and  gases  using  the  three  simplified  models  of  dissociation 
that  were  introduced  in  Chapter  IV: 

1.  Standard  Ladder  Model 

2.  Coupled  Model 

3.  Uncoupled  Model 

These  models  are  established  to  permit  an  evaluation  of 
the  influence  of  rotational  effects  and  complex  collision 
phenomena  on  the  dissociation  process  and  represent  differ- 
ent limiting  cases,  or  approximations  to  the  fully  coupled 
problem.  By  examining  the  rates  calculated  under  the  vari- 
ous models,  we  may  establish  a model  of  dissociation  that 
includes  all  essential  elements.  This  model  should  provide 
reasonably  accurate  values  for  dissociation  and  recombina- 
tion coefficients,  exhibit  an  acceptable  temperature  depen- 
dence, and  permit  an  evaluation  of  rate  enhancement  effects 
in  pure  gases. 

Within  the  framework  of  the  three  dissociation  models, 
the  following  capabilities  were  established  and  are  sche- 
matically presented  in  Table  V-l.  A time-dependent  capabil- 
ity represents  the  time-dependent  integration  of  the  coupled 
rate  equations  utilizing  the  Gear  (ref  85)  integration 
scheme  for  stiff  systems  of  equations.  The  steady-state 
solutions  utilize  equation  (3-3)  for  the  ladder  and  coupled 
nodels,  and  a rotational  average  of  equation  (3-25)  for  the 
ladder*  and  uncoupled  model.  The  "Dilute"  column  refers  to 

*No  rotational  average  is  needed  for  the  ladder  model. 
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calculations  of  the  dissociation  rate  for  cases  in  which 


the  diatomic  gas  is  diluted  in  a monatomic  gas  such  as 
argon.  This  is  an  interesting  case  to  consider,  not  only 
because  extensive  experimental  data  exist  for  these  cir- 
cumstances, but  because  the  collisional  processes  are 
especially  simplified.  Under  these  conditions  the  energy 
transfer  process  and  collisional  activation  are  dominated 
by  V-T  exchanges  with  the  buffer  gas,  and  complex  exchanges 
can  be  neglected.  In  contrast,  the  "Pure"  column  refers 
to  calculations  in  pure  diatomic,  (undiluted) , gases. 

Under  these  circumstances,  in  addition  to  the  change  from 
the  dilute  case  in  the  colliding  masses  for  V-T  exchanges, 
the  influence  of  complex  collisions,  V-V,  must  be  included 
and  analyzed.  The  designation  "A"  or  "C"  refers  to  the 
critical  dissociation  states  considered,  "A"  referring  to 
the  "over  the  top,"  ott,  model,  while  "C"  designating  the 
pseudo-barrier  penetration  resulting  from  equal  energy 
vibrational- rotational-translation,  V-R-T,  transitions.  It 
should  be  noted  that  no  detailed  time-dependent  solutions 
for  the  uncoupled  model  were  performed  because  it  would 
involve  the  solution  for  the  full  set  of  vibrational- 
rotational  states,  (over  300  for  Hj  and  50,000  for  I2) . 

The  uncoupled  model  was  not  analyzed  for  the  pure  case, 
due  to  the  linear  form  assumed  in  the  development  of  (4-11) 
Complex  exchanges  could  have  been  included  in  this  develop- 
ment and  the  resulting  equations  linearized,  but  this 
was  not  considered  necessary  due  to  the  agreement  obtained 
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between  the  coupled  and  uncoupled  models  for  the  "dilute" 
cases. 

Utilizing  this  matrix  of  calculational  capabilities, 
we  have  analyzed  the  dissociation  process  in  H2 , 02,  N2, 

CO,  Br2,  I2  considering  dilute  and  pure  gases.  Let's  dis- 
cuss these  results  by  reviewing  the  time  dependent  evolu- 
tion of  the  vibrational  distribution  and  the  dissociation 
rate.  Within  this  discussion,  we  will  demonstrate  the 
existence  of  a quasi-steady-state,  justifying  the  applica- 
tion of  a steady-state  analysis.  Then  the  incubation  times 

of  the  dissociation  reaction  will  be  examined. 

I 

Time  Dependent  Results 

We  wish  to  examine  the  shock  heating  of  a diatomic  gas 

t 

and  the  resulting  thermal  relaxation  and  dissociation. 

This  heating  is  modeled  as  an  instantaneous  rise  in  the 
translational  and  rotational  temperatures  occurring  as  the 
shock  passes,  followed  by  a period  of  vibrational  excita- 
tion and  relaxation.  During  this  relaxation  period,  the 
vibrational  manifold  is  populated  by  heavy  particle  colli- 
sions and  rapidly  attains  a quasi-steady-state.  The 
vibrational  "relaxation"  is  calculated  through  a numerical 
integration  of  the  coupled  vibrational  rate  equations 
using  the  SSH  (ref  20)  exchange  rates  as  modified  by  K„ck 
and  Carrier  (ref  43)  and  fitted  to  the  data  of  Millikan 
and  White  (ref  5)  and  of  Hancock  and  Smith  (ref  33) . The 
purpose  of  the  time-dependent  solution  is  to  establish  the 
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existence  of  a quasi-steady-state  vibrational  distribution, 
to  calculate  the  associated  incubation  time,  and  to  be  able 
to  verify  any  subsequent  steady-state  calculations . 


Before  considering  the  details  of  the  relaxation  pro- 
cess in  terms  of  the  vibrational  populations,  let  us  examine 
for  illustrative  purposes  the  typical  time  profiles  of  the 
molecular  and  atomic  densities  as  given  in  Fig.  5-1,  where 
we  have  adopted  this  data  from  the  calculations  of  Fisher 
(ref  86)  and  Bauer  and  Hilden  (ref  87).  The  solid  lines 
represent  the  molecular  and  atomic  densities.  Superimposed 
on  the  plot  is  the  unidirectional  flux  coefficient,  k^, 
where: 

\ 

4= 

v=o 

Note  that  the  dissociation  process  can  be  divided  into  three 
characteristic  time  regimes:  (1)  the  vibrational  relaxa- 
tion regime  in  which  the  molecule  is  "heated"  through  col- 
lisional  energy  exchanges  and  where  insignificant  dissocia- 
tion occurs,  (2)  the  quasi- steady  state  regime  (Q.S.S.) 
where  the  state  populations  have  obtained  a quasi-steady 
balance  between  collisional  activation  and  dissociation, 
and  where  recombination  is  negligible,  and  (3)  the  chemical 
equilibrium  regime  where  atomic  and  molecular  concentrations 
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Fig  5-1 


Time  Profiles  and  Uni-Directional  Flux 


have  achieved  their  equilibrium  values  and  state  popula- 
tions are  Boltzmann.  The  characteristic  times  associated 
with  the  first  two  regimes  are  the  vibrational  relaxation 
time  and  the  incubation  time.  The  vibrational  relaxation 
time  is  defined  by 


c/Etot  _ CE of  Etot(€l) 
d t 


in  which  Etot  (t)  is  total  vibrational  energy  and  Ewn  is 


EQ 


the  equilibrium  value  of  this  energy.  The  incubation  time 
x^nc,  is  defined  as  the  onset  time  for  Q.S.S.  dissociation 
Observe  that  the  phenomenalogical  rate  equation 


-4 

applies  in  the  region  from  t=Tinc  to  t=10  seconds  and 
that  a steady-state  analysis  is  justified  under  these  cir- 
cumstances. Also  note  that  the  chemical  equilibrium  value 
of  the  dissociation  flux  is  significantly  larger  than  the 
Q.S.S.  value,  reflecting  the  influence  of  depletion  on  the 
Q.S.S  and  the  increased  atomic  dissociation  efficiency. 

Let  us  consider  Regimes  i and  II  and  examine  the 
relaxation  process  and  ensuing  quasi-steady  state.  At 
time=0,  the  translational  and  rotational  modes  assume  a 
high  temperature,  (T=4000°K,  in  this  example). 


99 


Subsequently,  the  vibrational  states  are  populated  via 


V-T  collisions  and  relax  toward  a distribution  compatible 
with  this  temperature  and  consistent  with  the  dissociation 
losses.  A quasi-steady-state  is  then  established.*  Dur- 
ing the  period  of  time,  indicated  as  starting  Tincubation' 
quasi-steady-state  dissociation  takes  place.  The  vibra- 
tional distribution  maintains  a self-similar  form  as  the 
vibrational  populations  decay  exponentially  with  the  same 
time  constant.  This  is  the  time  interval  over  which  most 
experimental  observations  are  made  and  from  which  dissocia- 
tion rates  are  inferred.  During  this  time  period,  the 
molecular  and  atomic  densities  are  far  from  their  equili- 
brium concentration  values  and  the  recombination  process 
does  not  have  a significant  effect  on  the  vibrational  dis- 
tribution or  atom  density.  As  dissociation  proceeds,  the 
molecular  density  decreases  and  the  atom  density  increases 
to  the  point  where  it  is  comparable  with  the  molecular 
density.  Near  this  point  a departure  from  the  quasi- 
steady-state, Q.S.S.,  occurs.  The  reasons  for  this  depar- 
ture are  twofold.  First,  as  discussed  earlier,  during  the 
Q.S.S.,  energy  losses  due  to  dissociation  modify  the  vibra- 
tional distribution  and  vibrational  depletion  results. 

This  modification  of  the  vibrational  distribution  is 

*At  very  high  temperatures,  th^  relaxation  and  dissoci- 
ation coupling  may  be  so  strong  that  the  separation  of  the 
relaxation  and  dissociation  regimes  may  be  invalid. 
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manifested  in  the  reduction  in  the  occupation  of  the  high 
lying  vibrational  states  below  their  equilibrium  values 
and  a consequent  reduction  in  the  dissociation  rate.  As 
recombination  begins  to  become  important,  the  net  loss  from 
the  vibrational  manifold  decreases,  ultimately  reaching 
zero  when  true  chemical  equilibrium  is  established.  There- 
fore, recombination  reduces  depletion  and  the  unidirec- 
tional dissociation  rate  is  enhanced  as  recombination 
becomes  important.  Hence  the  approach  to  equilibrium  is 
accelerated  by  the  recombination  process. 

The  second  reason  for  the  departure  from  the  Q.S.S. 
is  related  to  the  dissociation  efficiency  of  atoms.  The 
product  atoms  are  very  efficient  collision  partners  with 
respect  to  the  parent  molecule,  with  0 atoms  being  almost 
an  order  of  magnitude  more  efficient  than  or  Ar  (ref 
93).  When  significant  atom  concentrations  are  achieved, 
they  will  accelerate  the  approach  to  chemical  equilibrium 
by  increasing  the  dissociation  flux  and  consequently  the 
atomic  density  in  a typical  "chain  reaction"  type  process. 
This  aspect  of  the  reaction  has  not  been  modeled  in  our 
calculations . 

With  the  background  provided  by  reviewing  the  molecu- 
lar and  atomic  profile  data  of  Fig.  5-1,  we  can  now  present 
selected  results  of  our  time  dependent  calculations  in 
Figs.  5-2  and  5-3.  These  results  reveal  the  details  of 
the  vibrational  relaxation  occurring  in  regime  I,  and 
quantify  the  concepts  of  Q.S.S.  and  vibrational  depletion 
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encountered  in  regime  II.  In  Fig.  5-2,  a semi- log  plot 
of  normalized  vibrational  population  as  a function  of 
normalized  vibrational  energy  is  presented  for  selected 
times  in  the  relaxation  process.  In  this  representation 
a Boltzmann  distribution  would  appear  as  a straight  line. 

The  vibrational  distribution  is  given  by  a Boltzmann  dis- 
tribution at  300 °K  at  time  t=0 . The  distribution  evolves 
through  a series  of  pseudo-Boltzmann  distributions  and 
achieves  a Q.S.S.  distribution  for  t>0 . 1 ysec.  These 
results  are  consistent  with  those  of  Montroll  and  Schuler 
(ref  55)  and  Kewley  (ref  74) . It  should  be  noted  that  the 
relaxation  accelerates  upon  achieving  a significant  popula- 
tion in  the  excited  vibrational  states  reflecting  the  rate 
scaling  of  the  V-T  rates  and  the  effect  of  molecular 
anharmonicity . The  Q.S.S  distribution  is  approximated 
well  by  a Boltzmann  distribution  at  the  gas  temperature 
for  E < 0.7  Dq,  however,  a significant  departure  is 
observed  for  E > 0.7  Dq.  This  departure  from  a Boltzmann 
distribution  is  caused  by  the  vibrational-dissociative 
coupling.  A vibrationally  depleted  Q.S.S.  distribution  is 
thus  established.  As  a result  of  vibrational  depletion, 
the  dissociation  rate  is  depressed  below  that  value  that 
would  be  obtained  if  the  distribution  were  undepleted  with 
the  rate  typically  reduced  by  a factor  of  10-20.  Fig.  5-3 
shows  the  results  of  similar  calculation  for  Hj . 

Thus  from  Figs.  5-2  and  5-3  we  have  established  that 
the  vibrational  distribution  rapidly  attains  a quasi-steady- 


state  and  that  a steady-state  analysis  is  reasonable  and 
justified.  In  addition,  we  have  illustrated  that  the 
vibrational-dissociative  coupling  results  in  a vibration 
depletion  and  a significant  rate  reduction. 

Let  us  now  direct  our  attention  to  the  calculation 
of  incubation  times.  Recall  from  Fig.  5-1,  that  the  incu- 
bation time  is  the  lag  time  associated  with  the  onset  of 
quasi-steady-state  dissociation.  Fig.  5-4  presents  the 
unidirectional  dissociation  flux  results  of  the  time  depen- 
dent calculations  of  (^“Ar  mixtures  at  densities  of 

[O2]  - Z.  67 X /o'6 partic /es /c/ri3 
£ArJ =Z.67x  / o/ 9 portictts/c/rj5 


and  at  the  temperatures  noted.  Observe  that  the  incubation 
time  decreases  monotonically  with  temperature  in  contrast 
to  the  dissociation  flux  which  is  increasing  with  tempera- 
ture. The  dashed  curve,  D',  illustrates  the  incubation 
profile  for  pure  02»  and  reveals  a reduction  in  the  incu- 
bation time  resulting  from  V-V  exchanges  in  addition  to 
the  increased  V-T  exchange  rates  for  °2_02  collis*-ons • 
Recombination  has  not  been  included  in  our  calculation, 
so  that  the  transition  to  Regime  III  will  not  be  observed. 

A comparison  of  the  calculated  incubation  times  with  the 
experimental  data  of  Watt  (ref  94)  and  Wray  (ref  95)  is 
presented  in  Fig.  5-5  where  the  results  have  been  normalized 
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results  obtained  in  this  study,  where  identical  V-V  and 
V-T  rates  were  used  in  both  model  calculations,  and  is 
also  supported  by  Wengle ' s (ref  12)  results. 

Dissociation  of  Diatomics  Diluted 
in  an  Argon  Heat  Bath 

Most  of  the  dissociation  data  from  shock  tube  experi- 
ments has  been  obtained  for  various  diatomic  species  highly 
diluted  in  argon.  This  is  an  extremely  simple  physical 
situation  from  a collisional  analysis  viewpoint,  in  that 
collisional  exchanges  are  limited  to  V-T  exchanges  with 
the  argon  buffer.  It  therefore  offers  a test  of  any  colli- 
sional theory  of  dissociation.  In  this  section,  we  present 
the  predicted  dissociation  rate  coefficients  for  the  dilute 
case  based  on  the  calculational  matrix  in  Table  V-l* 

Time  dependent  solutions  for  the  ladder  model  (A) , and 
coupled  model  (A/C)  showed  agreement  within  0.1%  with  the 
steady-state  analysis,  based  on  equation  (3-13),  for  all 
gases  and  for  the  entire  range  of  temperatures  presented 
in  this  study.  Since  this  agreement  was  within  the  conver- 
gence criteria  established  for  the  iterative  solution,  it 
serves  as  a validation  of  these  calculations.  A comparison 
of  time  dependent  results  with  the  steady-state  analysis 
of  equation  (3-25)  also  yielded  acceptable  agreement  with 
the  maximum  deviation  being  less  than  2%.  In  addition  to 
this  cross-validation  of  codes,  each  calculation  was 
checked  for  the  following  limiting  cases. 
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1.  No  dissociation 

2.  dv=const. 

For  both  of  these  cases,  the  required  Boltzmann  distribu- 
tion was  obtained. 

Having  documented  the  calculations,  let  us  compare 
the  dissociation  coefficient  calculated  under  the  three 
different  models, 

1.  standard  ladder, 

2.  coupled,  and 

3.  uncoupled, 

to  experimental  data  for  mixtures  of  a diatomic  gas  highly 
diluted  in  argon.  From  this  comparison,  we  will  establish 
three  significant  points: 

1.  The  standard  ladder  model  is  inadequate,  yielding 
dissociation  coefficients  that  are  much  too  low. 

2.  Models  including  the  influence  of  rotation  pre- 
dict larger  dissociation  coefficients  than  the  standard 
model  and  these  coefficients  are  in  substantial  agreement 
with  experimental  data. 

3.  Equal  energy,  V-R-T,  transitions  lead  to  an 
increase  in  the  dissociation  rates;  this  increase  is 
especially  marked  in  Br2  and  I2  and  yields  improved  agree- 
ment with  experimental  data. 

For  all  gas  mixtures  considered  here,  the  argon  number 

19-3  . . 

density  was  2.67  x 10  cm  while  the  diatomic  density  was 

2.67  x lO1^  or  lower.  Under  these  conditions,  V-T 
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processes  dominate  the  V-V  processes  for  all  temperatures 
and  diatomic  gases  considered. 


Hydrogen.  Hydrogen  has  been  studied  extensively  both 
experimentally  and  theoretically.  Sufficient  experimental 
data  is  available  to  establish  the  dissociation  coefficient 
to  within  a factor  of  two.  A comparison  of  experimental 
(ref  96)  and  theoretical  values  is  presented  in  Fig.  5-6 
where  the  logarithm  of  the  dissociation  coefficient  is 
plotted  against  the  gas  temperature.  The  theory  is  shown 
in  solid  lines,  experiment  in  dashed.  (This  convention 
will  be  utilized  throughout  the  entire  data  comparison.) 

All  calculated  values  fall  below  the  experimental  data, 
but  exhibit  an  acceptable  temperature  dependence.  The 
ladder  model  is  a factor  of  6 low.  All  rotational  models 

1.  coupled,  equal  energy  transition,  CC; 

2.  coupled,  "over  the  top,"  CA; 

3.  uncoupled,  equal  energy  transition,  UC;  and 

4.  uncoupled,  "over  the  top,"  UA 

yield  essentially  identical  results  and  are  represented  by 
the  curve  marked  "ROT."  Rotation  enhances  the  dissocia- 
tion coefficient  in  hydrogen  by  a factor  of  3.  This  factor 
is  much  less  than  that  anticipated  from  our  equilibrium 
analysis,  and  reflects  the  fact  that  the  vibrational  trans- 
fer rates  are  not  large  enough  to  support  the  full  rota- 
tion enhancement  and  theref o • 3 are  rate  controlling. 
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This  "bottle-necking"  effect  and  consequent  depletion  has 
been  pointed  out  by  Pritchard  (ref  57) . 


Oxygen  and  Nitrogen.  The  dissociation  of  oxygen  and 
nitrogen  has  received  much  attention,  primarily  in  order  to 
elucidate  the  chemical  behavior  of  air  at  high  temperatures. 
The  theoretical  results  for  oxygen  are  compared  with  experi- 
ment in  Figs.  5-7  and  5-8.  Fig.  5-7  parallels  the  format 
and  trends  noted  in  Fig.  5-6  for  hydrogen.  Including 
rotation  increases  the  dissociation  coefficient  once 
again.  For  oxygen,  we  are  now  able  to  observe  slight  dif- 
ferences among  the  rotational  models.  These  differences 
will  become  more  pronounced  in  N2 , CO,  Br2 > I2,  but  the 
general  ordering  will  remain  the  same.  In  Fig.  5-7, 

CC  > CA  >UC  > UA.  We  will  see  that  in  general  CC  > CA  and 
UC  > UA.  This  is  consistent  with  our  equilibrium  analysis 
of  CDS  and  reflects  the  fact  that  equal  energy  transitions 
provide  an  easier  path  for  the  dissociation  process.  The 
coupled  curves,  CC  and  CA,  are  higher  than  the  uncoupled, 

UC  and  UA,  due  to  the  more  uniform  distribution  of  the 
dissociation  rates  with  respect  to  the  vibrational  levels 
in  the  coupled  model.  Stated  simply,  the  vibrational  bias 
in  the  coupled  model  is  less  than  that  in  the  uncoupled 
model.  This  lower  vibrational  bias  reduces  the  vibrational 


depletion  and  leads  to  a larger  dissociation  rate.  Fig. 
5-8  is  an  Arrhenius  plot  over  a wider  temperature  range, 
where  the  logarithm  of  the  dissociation  coefficient  is 


k<j,  cm3  sec 


Temperature,  °K 

Fig  5-7  Dissociation  of  Oxygen  Diluted  in  Argon 
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plotted  as  a function  T An  additional  temperature  scale 
is  provided  at  the  top  to  facilitate  the  conversion.  The 
rotational  model  results  are  bounded  by  the  CC  and  UA 
curves  and  show  good  agreement  with  experiment  over  this 
wide  range  of  temperatures.  Similar  results  are  achieved 
for  nitrogen,  Fig.  5-9. 

Bromine  and  Iodine.  The  dissociation  of  bromine  has 
been  studied  by  Palmer  and  Hornig  (ref  102)  where  they 
employed  a spectrophotometric  measurement  method.  Strong 
(ref  103)  performed  flash  photolysis  experiments  at  300° 
to  430 °K  to  obtain  recombination  coefficient  data  which  can 
be  converted  via  the  equilibrium  constant  to  equivalent  dis- 
sociation rates.  A comparison  of  the  theory  with  these 
experimental  data  is  given  in  Fig.  5-10.  The  ladder  model 
fails  badly  in  bromine.  Note  that  the  ladder  model  results 
as  presented  have  been  multiplied  by  a factor  of  100.  The 
UA  model  is  a factor  of  50  low,  whereas  CC,  CA,  UC  yield 
more  reasonable  agreement.  Equal  energy  transitions  pro- 
vide a more  substantial  rate  enhancement  in  bromine,  than 
that  previously  observed  in  02  and  N2.  Similar  results  are 
obtained  for  iodine.  Fig.  5-11.  The  ladder  model  is  again 
over  a factor  of  100  too  low  with  equal  energy  transitions 
increasing  the  rate  in  the  rotational  models  by  almost  an 
order  of  magnitude.  Enhancement  due  to  eet’s  is  large 
enough  to  result  in  the  UC  curve  exceeding  the  CA  curve. 
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This  reveals  that,  in  iodine,  equal  energy  transitions  and 
vibrational  bias  are  of  equal  importance. 

Carbon  Monoxide.  Carbon  monoxide  is  very  similar  in 
mass  and  vibrational  energy  spacing  to  oxygen  and  nitrogen, 
and  hence  one  would  suspect  that  the  dissociation  coeffi- 
cient comparison  would  parallel  those  in  Figs.  5-8  and  5-9. 
This  is  not  the  case.  Fig.  5-12  shows  an  Arrhenius  plot 
of  the  measured  dissociation  coefficients  of  carbon  monox- 
ide (refs  105,  106,  107)  together  with  the  theoretical  pre- 
dictions. While  the  ordering  of  the  theoretical  predic- 
tions are  very  similar  to  Nj  and  O2/  the  predictions  of 
the  present  models  fall  substantially  below  the  experimental 
data.  The  ladder  model  is  unacceptably  low  again.  The 
\ reason  for  this  discrepancy  may  be  related  to  the  fact  that, 

in  our  calculations,  we  have  considered  only  the  ground 
electronic  state  of  the  molecule.  We  have  chosen  this 
course,  as  mentioned  previously,  due  to  the  scarcity  of 
data  on  e-V  and  e-T  exchanges  and  V-T  and  V-V  exchanges 
for  electronically  excited  states.  However  in  CO,  it  has 
been  pointed  out  by  Fairbairn  (ref  108)  and  by  Appleton, 
Steinberg,  and  Liquornik  (ref  105)  that  excited  electronic 
states  could  play  a significant  role  in  the  dissociation 
reaction.  This  participation  of  electronically  excited 
states  was  inferred  from  the  density  dependence  of  the 

fv 

induction  time.  It  is  difficult  to  accurately  estimate  the 
magnitude  of  the  enhancement  due  to  electronically  excited 
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states  under  the  models  considered  in  this  study.  First 
one  must  know  the  rates  and  spectroscopic  constants  for 
each  state,  with  the  rates  being  the  limiting  factor. 


Second,  and  related  to  the  rates,  we  must  establish  the 
collisional  coupling  distribution,  (e-V,  e-T) , between 
electronic  states.  As  a first  approximation  and  upper 
bound,  one  could  simply  multiply  the  rate  for  the  ground 
electronic  state  by  the  degeneracies  of  all  electronic 
states  and  sum.  This  would  yield  a correction  factor  of 
about  20.  Shui  (ref  109)  has  performed  a more  detailed 
evaluation  of  this  enhancement  factor  within  the  framework 
of  Keck's  statistical  theory  and  obtained  a factor  of  8. 
Certainly  the  enhancement  should  be  temperature  dependent; 
however,  for  lack  of  a better  estimate,  we  could  apply 
Shui ' s correction  factor  directly  to  our  predictions  and 
observe  that  much  better  agreement  with  experiment  is  then 
obtained.  We  would  therefore,  tend  to  agree,  based  on  the 
numerical  experiments,  that  there  is  an  excited  state 
enhancement  in  CO  dissociation  rates. 

» Dissociation  in  Pure  Gases 

. 

The  rotational  models  have  yielded  good  predictions 
for  a wide  variety  of  gases  when  we  considered  the  case 
where  the  diatomic  gas  was  highly  diluted  in  an  argon 
buffer  gas.  The  importance  of  equal  energy  transitions, 
especially  in  B^,  and  1 2,  has  been  pointed  out  as  well  as 
the  role  of  electronically  excited  states  in  CO.  We  now 
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consider  dissociation  in  pure  gases.  We  will  omit  the 
ladder  model  in  future  discussions  due  to  its  failure  in 
the  dilute  case. 

In  general,  the  dissociation  rate  in  pure  gases  is 
larger  than  that  obtained  when  the  diatomic  is  highly 
diluted.  The  extent  of  this  increase  varies  from  gas  to 
gas,  being  about  a factor  1.5  for  B^-B^  and  a factor  of 
10  for  02-02"  Th^s  rate  enhancement  in  pure  gases  has 
been  investigated  within  the  calculational  matrix  of  Table 
V-l,  where  the  influence  of  vibrational-to-vibrational 
exchanges  and  changes  in  collisional  efficiencies  due  to 
mass  variation  have  been  included.  From  this  study,  it 
was  established  that:  (1)  changes  in  the  collisional 
efficiencies  do  not  properly  account  for  either  the 
observed  magnitude  of  the  rate  enhancement  or  its  tempera- 
ture dependence,  and  (2)  V-V  exchanges,  however,  lead  to 
rate  enhancements  similar  in  magnitude  and  temperature 
dependence  to  those  experimentally  observed. 

Hydrogen.  The  Arrhenius  plots  of  the  dissociation 
coefficient  for  H2-H2  and  ^-Ar,  illustrated  in  Fig.  5-13, 
show  a comparison  of  exprimental  and  theoretical  results. 
The  experimental  results  of  Breshears  and  Bird  (ref  110) 
differ  from  those  of  Sutton  (ref  96).  Sutton's  data  show 
the  ratio  (H2-H2)  /k^  (^-Ar)  to  be  almost  constant  and 
equal  to  4.  This  is  in  contrast  to  the  more  recent  data 
of  Breshears , which  reveal  a temperature  dependence  to  the 
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ratio  varying  from  3 at  3500°K  to  10  at  5000°K.  Three 
theoretical  predictions  are  shown.  The  H2~argon  curve 
duplicates  the  rotational  results  obtained  for  the  dilute 
case  which  were  presented  previously  in  Fig.  5-6  in  a dif- 
ferent format.  The  curve  marked  H2  illustrates  the  effect 
of  changes  in  colliding  masses,  (H2~H2  vs  H2~Ar) , and  does 
not  include  the  influence  of  V-V  exchanges.  The  H2~H2 
curve  includes  both  V-V  effects  and  mass  changes,  where  we 
have  utilized  the  theory  of  Rapp  and  Englander-Golden  (ref 
111)  for  the  V-V  rates. 

The  increase  in  the  rate  constant  due  to  mass  changes 
was  less  than  the  observed  experimental  enhancement  and 
insensitive  to  variations  in  temperature.  V-V  exchanges, 
however,  produced  a more  significant  increase  in  the  rate 
constant  and  introduced  a temperature  dependence  very 
similar  to  that  observed  by  Breshears  (ref  110) . The 
observed  rate  enhancement  due  to  V-V  exchanges  is  consistent 
with  the  calculations  of  Bauer  and  Hilden  (ref  87) , who 
used  an  empirical  form  for  the  dissociation  rates.  The 
magnitude  of  the  enhancement  is  fairly  sensitive  to  the 
magnitude  of  the  V-V  rates.  Direct  application  of  the 
Rapp  form  yielded  a factor  of  2.2  enhancement  at  3000 °K 
and  a factor  of  4.4  at  5000 °K;  whereas  if  the  V-V  rates 
were  multiplied  by  ten,  the  enhancement  changed  to  10  and 
30.  Shown  in  the  illustration  are  results  for  the  Rapp 
rate  multiplied  by  3.  Further  theoretical  and  experimental 
work  is  required  to  establish  the  true  magnitude  and 


125 


temperature  dependence  of  the  V-V  exchange  rates.  Summar- 
izing, we  conclude  that  V-V  exchanges  can  account  for  the 
observed  rate  enhancement  in  hydrogen. 


Oxygen . We  now  examine  oxygen  to  see  if  the  results 
obtained  for  H2  find  general  applicability.  A comparison 
of  these  experimental  and  theoretical  results  is  provided 
in  Fig.  5-14.  The  experimental  data  of  Byron  (ref  97) 
shows  that  the  ratio  k^  ((^-C^)  /k^  (02~Ar)  is  about  ten  and 
independent  of  temperature  over  the  temperature  range 
reported.  The  theoretical  predictions  yielded  the  follow- 
ing dissociation  rate  ratios  for  the  values  of  V-V 
exchanges  listed: 


W X f 

kd(02"°2) 

kd(02-A) 

f = 1 

3.0 

f = 3 

3.8 

f = 10 

10 

The  curve  for  f=3  is  presented  in  the  illustration  where 
it  should  be  pointed  out  that  mass  enhancement  alone 
increased  the  dissociation  rate  by  a factor  of  two.  It  is 
gratifying  that,  with  minor  adjustments  in  the  V-V  rates, 
rate  enhancements  in  both  oxygen  and  hydrogen  are  con- 
sistently accounted  for.  This  even  includes  the  tempera- 
ture dependence  of  the  observed  enhancement. 
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Bromine.  Bromine  is  an  interesting,  but  difficult 
case  to  analyze.  It  is  interesting  because  Br2  is  heavier 
than  argon;  and  therefore  from  a collisional  standpoint 
would  be  expected  to  have  a lower  efficiency  for  V-T  trans- 
fer and  result  in  lower  dissociation  coefficient.  However, 
an  enhancement  of  approximately  a factor  of  two  is  observed. 
The  difficulty  with  Bromine  lies  in  the  uncertainty  in  the 
relaxation  rate  data.  Relaxation  data  is  available  only 
for  Br2~Ar  and  these  data  exist  only  for  temperatures  below 
1000 °K.  Extrapolation  to  higher  temperatures  coupled  with 
a mass  scaling  interpolation  to  Br2~Br2  may  introduce  sig- 
nificant uncertainties  into  the  rate  data.  This  fact 
coupled  with  V-V  uncertainties  must  lead  us  to  view  the 
following  results  with  moderated  enthusiasm.  Arrhenius 
plots  of  the  dissociation  coefficient  are  given  in  Fig. 

5-15  where  the  Br2~Ar  curve,  CC,  represents  our  previous 
dilute  results.  The  experimental  data  of  Hornig  (ref  102) 
reveals  a factor  of  two  enhancement.  Mass  effects  lead  to 
a 30%  decrease  in  the  calculated  dissociation  rate;  how- 


ever, when  V-V  exchanges  are  included  the  rate  is  increased, 
yielding  a temperature  dependence  similar  to  that  obtained 
for  H2~H2.  Agreement  with  experiment  is  disappointing  and 
may  reflect  the  previously  cited  rate  inadequacies.  Addi- 
tional experimental  relaxation  and  dissociation  data  are 
needed  to  provide  further  knowledge  about  the  factors  con- 
trolling dissociation  in  Br2- 
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Fig  5-15  Rate  Enhancement  in  Brj 
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Recombination  Coefficients 

Having  established  that  the  rotational  models  have  a 
good  predictive  capability  regarding  dissociation  coeffi- 
cients in  dilute  and  pure  gases,  we  now  examine  the  recom- 
bination coefficients,  kr,  that  would  be  inferred  by 

J 

relating  the  rates  through  the  equilibrium  constant,  K . 

The  equilibrium  constant,  K, , for  the  molecule,  AB,  is  given 
by  (ref  37) , 

I 

with  nA  = number  of  density  of  A atoms 
n„  = number  density  of  B atoms 

D 

n^  = number  density  of  AB  molecules 

fA  = (Z77mA*7fzfA/°A 

fB  ^ 

and  mA,  m^  are  the  masses  of  the  atoms,  gA  and  gg  are  the 
electronic  weights  of  their  normal  states;  pA  and 
are  the  nuclear  spin  weights.  Pertinent  data  is  presented 
in  Table  V-2.  The  molecular  partition  function,  fAB> 

(ref  37)  is  given  by; 

4=  ,-rot 
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with  o-n  the  symmetry  number. 

AS 

ling  the  nuclear  spin  weights  (ref  37) , 

mt 

ft  _ (*V*a  f»e*T)  Z<s~a*  fA 

° (W*'?*  h 3 Q,lirrot  f*m 

Since  Kc  a exp(-$Do)  and  is  approximately  a exp(-BDQ) 
the  recombination  coefficient 


strips  away  the  exponential  factor  in  the  dissociation 
coefficient  and  provides  a detailed  way  of  examining  the 
activation  energy  or  the  pre-exponential  temperature  depen- 
dence. In  looking  at  the  recombination  coefficients,  we 
will  examine  the  UC  and  UA  models  only,  since  these  models 
behaved  similarly  to  the  CC  and  CA  models  differing  mainly 
in  magnitude  and  exhibiting  only  a slightly  stronger, 

(lower  activation  energy) , temperature  dependence.  The 
following  theoretical  predictions  were  made  for  the  dilute 
case,  assuming  that  argon  was  the  collision  partner.  This 
was  the  case  in  the  majority  of  experimental  investigations 
chosen  for  comparison.  In  addition,  we  continue  to  include 
only  the  contribution  of  the  ground  electronic  state, 

(see  Table  V-2) . Our  major  concern  in  all  these  comparisons 
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will  be  the  temperature  dependence,  since  we  have  noted 
variations  in  magnitude  of  the  various  models  previously 
when  we  discussed  the  dissociation  coefficients.  We  will 
observe  throughout  the  recombination  comparisons  that  the 
predicted  temperature  dependence  is  extremely  good, 
especially  when  one  considers  the  temperature  range,  rate 
scaling,  and  variety  of  gases  considered. 

Hydrogen.  The  predicted  kr(H2>,  again  given  by  the 
solid  line,  is  compared  with  the  recombination  coefficients 
inferred  from  the  high  temperature  dissociation  rate  mea- 
surements (refs  112,  113,  114,  115)  and  the  direct  recombina- 
tion rate  measurements  (refs  116,  117)  in  Fig.  5-16.  The 
agreement  between  our  prediction  and  the  high  temperature 
measurements  of  Sutton  (ref  114) , Myerson  (ref  115) , Jacobs 
(ref  112),  and  Hurle  (ref  113)  is  good,  reproducing  the 
temperature  dependence  despite  being  about  a factor  of  three 
low.  At  lower  temperatures,  the  temperature  dependence  is 
too  strong,  and  may  reflect  a breakdown  in  the  assumption 
of  rotational  equilibrium,  an  inadequacy  in  the  scaling  of 
the  exchange  rates  to  low  temperatures,  or  an  increased 
sensitivity  to  the  molecular  potential  and  state  designa- 
tion at  low  temperatures.  The  results  obtained  for  hydrogen 
showed  the  largest  disparity  of  all  gases  considered. 

Oxygen.  The  predicted  rate  coefficients  for  oxygen 
are  presented  in  Fig.  5-17  and  are  compared  to  the  high 
temperature  dissociation  data  of  Wray  (ref  98) , Camac  and 
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Vaughan  (ref  93) , and  Watt  and  Myerson  (ref  118) . Note 
the  disparity  in  the  temperature  dependence  between  the 
experimental  data  of  Watt  and  Myerson  and  Wray  and  Watt. 

Our  calculations  favor  the  latter  dependence  and  exhibit 
good  agreement  with  that  data.  This  agreement,  obtained 
for  both  models  considered,  extends  down  to  the  low  tempera- 
tures where  it  compares  favorably  with  the  direct  measure- 
ments of  Campbell  and  Thrush  (ref  119) . No  significant 
difference  in  temperature  dependence  exists  between  the  UC 
and  UA  curves. 

Nitrogen.  The  agreement  obtained  in  nitrogen  is  good 
and  comparable  to  that  obtained  in  oxygen.  Fig.  5-18 
illustrates  this  comparison  with  the  low  temperature  point 
measurements  of  Campbell  and  Thrush  (ref  120)  and  Clyne 
and  Steadman  (ref  121) , and  the  high  temperature  data  of 
Appleton  (ref  122) , Byron  (ref  123)  and  Cary  (ref  124) . 
Observe  that  at  low  temperatures  the  UC  curve,  incorporating 
equal  energy  transitions,  exhibits  a stronger  temperature 
dependence  than  the  UA  curve  representing  the  "over  the 
top"  model,  and  yields  substantially  better  agreement  with 
experiment. 

Bromine . The  data  comparison  in  bromine.  Fig.  5-19, 
reinforces  the  idea  that  equal  energy  transitions,  repre- 
sented by  UC,  play  a significant  role  in  the  dissociation 
reaction.  The  point  data  at  low  temperatures  are  due  to 
Ip  and  Burns  (refs  125,  126,  127),  while  the  higher 
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temperature  data  represented  as  areas  A and  B are  from  the 
data  summary  given  in  (ref  127)  and  the  emission  measure- 
ments of  Boyd  et  al.  (ref  128)  , respectively.  The  UA 
curve  has  been  multiplied  by  a factor  of  three  to  facili- 
tate the  comparison.  The  UC  curve  again  exhibits  a 
stronger  temperature  dependence  than  UA  and  yields  excel- 
lent agreement  with  the  experimental  data. 

Iodine.  The  predicted  rate  coefficients  for  iodine 
are  presented  in  Fig.  5-20  and  confirm  the  equal  energy 
trend  established  in  nitrogen  and  bromine.  Illustrated  is 
the  point  data  of  Porter  and  Smith  (ref  129)  at  low  tempera- 
tures along  with  the  data  of  Bunker  (ref  130) . The  high 
temperature  data  are  that  of  Britton  et  al.  (ref  131)  and 
Troe  (ref  132) . The  choice  here  is  clear.  The  equal  energy 
transition  model,  not  only  yields  a more  reasonable  magni- 
tude for  the  rate  coefficient  but  also  yields  a temperature 
dependence  consistent  with  the  experimental  data  over  the 
entire  temperature  range. 


Summary 

\ The  thermal  dissociation  process  has  been  parameter- 

ized over  a wide  range  of  temperatures  and  gases.  In  this 
particular  study,  we  have  used  three  simplified  models. 

•I 

1.  Standard  ladder 

2.  Coupled  Model 

3.  Uncoupled  Model 

to  evaluate  the  rotational  contribution  to  dissociation. 


140 


Time-dependent  and  steady-state  solutions  were  obtained  and 
compared  with  experimental  data.  These  comparisons  allowed 
us  to  identify  the  standard  ladder  model  as  an  inadequate 
representation.  The  remaining  rotational  models  yielded 
substantially  more  accurate  predictions  of  the  dissociation 
coefficient. 

The  implications  of  equal  energy  transitions,  (eet) , 
were  examined  within  the  framework  of  the  rotational  models 
and  contrasted  with  an  "over  the  top,"  ott,  mechanism.  We 
established  that  equal  energy  transitions  enhance  the  dis- 
sociation coefficient,  (significantly  in  Bromine  and  Iodine) , 
and  decrease  the  activation  energy  with  respect  to  the  ott 
model,  to  yield  improved  agreement  with  experimental  data. 

An  investigation  of  rate  enhancement  in  pure  gases 
revealed  that  changes  in  collisional  efficiencies  based  on 
mass  variations  can  not  account  for  the  observed  rate 
enhancement.  However,  V-V  exchanges  do  yield  a rate 
enhancement  consistent  with  that  observed  experimentally. 

A comparison  of  recombination  coefficients,  calculated 
from  the  rotational  models,  with  experimental  data  demon- 
strated that  these  models  are  able  to  correlate  both  the 
very  high  and  low  temperature  rate  measurements  for  a 
variety  of  gases.  In  addition,  this  comparison  again 
revealed  the  consistency  of  including  the  eet  mechanism 
in  the  calculations. 
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Chapter  VI.  Dissociation  Under  Highly 
Non-Equi librium  Conditions 


Introduction  to  Non-Equilibrium 
Processes 

In  the  preceding  chapters,  time-dependent  and  steady- 
state  solutions  of  the  Master  Equation  have  been  utilized 
in  an  investigation  of  thermal  dissociation  and  recombina- 
tion. This  analysis  of  elementary  processes  in  excited 
molecules,  including  energy  relaxation  and  vibrational 
exchange,  plays  an  important  role  in  studies  of  combustion 
and  gas  dynamics , but  deals  with  systems  that  are  only 
slightly  perturbed  from  thermal  equilibrium.  In  recent 
years,  however,  there  have  been  significant  advances  and 
interest  in  the  fields  of  chemical  kinetics  and  molecular 
lasers  where  highly  non-equilibrium  situations  are  encoun- 
tered. These  highly  non-equilibrium  conditions  can  be 
achieved  by  electron  impact  excitation  or  radiative  pumping 
of  the  vibrational  manifold.  Given  the  success  of  the  pre- 
sented collisional  analyses  in  addressing  dissociation  under 
essentially  equilibrium  conditions,  it  is  reasonable  to 
attempt  to  extend  these  theories  to  treat  those  highly  non- 
equilibrium conditions  encountered  in  electric  gas  dis- 
charges or  optically  pimped  systems.  In  this  chapter,  we 
will  utilize  the  simple  chemical  reaction  of  dissociation 
to  illustrate  the  effect  of  non-equilibrium  on  chemical 
processes.  We  will  investigate  the  dissociation  process 
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under  highly  non-equilibrium  conditions  and,  for  interest- 
ing plasma  conditions,  compare  the  dissociation  caused  by 
electron  impact  with  the  dissociation  due  to  heavy  particle 
collisions. 

In  this  section  of  the  study,  we  will  consider  only  a 
gas  mixture  of  carbon  monoxide  and  argon.  The  analysis  and 
conclusions,  however,  have  general  application  to  other 
molecular  gases  and  gas  mixtures.*  The  CO/Ar  system  is  an 
attractive  one  to  analyze  since  abundant  data  is  available 
on  electron  impact  cross-sections,  and  the  molecular  trans- 
fer rates  are  reasonably  well  known.  This  is  not  a minor 
point,  since  the  electron  energy  distribution,  vibrational 
energy  distribution,  and  dissociation  rates  are  sensitive 
functions  of  these  input  data.  In  addition,  the  CO/Ar 
system  is  experimentally  versatile;  capable  of  being  either 
optically  pumped  or  excited  in  a gas  discharge.  Reliable 
experimental  data  is  available  for  both  situations,  thus 
permitting  a benchmarking  of  the  calculations.  Here  we 
will  make  comparisons  with  the  experimental  distributions 
obtained  by  Rich  (ref  133)  who  utilized  direct  optical  pump- 
ing of  CO.  Under  conditions  of  either  discharge  or  optical 
pumping,  efficient  excitation  of  CO  is  achieved  and  a highly 
non-equilibrium  state  is  obtained. 

To  address  non-equilibrium  dissociation,  we  utilize 
for  the  vibrational  manifold  of  the  molecules  the  Master 


I i 


♦Selected  calculations  have  also  been  performed  in 
CO/He,  0o/Ar,  N./Ar,  and  0_/He  mixtures  at  various  concentra- 
tions. 2 2 2 
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Equation  formalism  previously  developed  for  the  near 
equilibrium  conditions;  and  for  the  electrons  we  use  a 
steady-state  solution  of  the  collisional  Boltzmann  equa- 
tion. The  presentation  of  these  detailed  calculations  will 
be  prefaced  by  a development  of  the  collisional  Boltzmann 
equation,  and  a discussion  of  the  relevant  electron- 
molecule  collisional  processes.  We  will  examine  the 
molecular  kinetics,  the  characteristics  of  the  vibrational 
distribution  under  strongly  pumped  conditions,  and  review 
an  analytic  approximation  to  the  vibrational  distribution. 
This  will  provide  us  with  a better  insight  when  we  subse- 
quently present  time-dependent  and  steady-state  calcula- 
tions of  the  vibrational  distribution  and  dissociation 
coefficient.  Finally,  we  will  link  the  two  analyses  in  a 
parametric  mapping  of  the  heavy  particle  and  electron  impact 
dissociation  processes  under  non-equilibrium  conditions. 

Electron  Kinetics 

Background . Under  appropriate  conditions,  electron 
impact  excitation  is  an  efficient  method  of  delivering 
energy  to  molecular  gases.  In  this  section,  we  will  pre- 
sent results  of  a numerical  investigation  of  the  electron- 
molecule  energy  transfer  processes  and  the  resulting  elec- 
tron energy  distributions  obtained  in  electrically  excited 
molecular  gases.  Emphasis  is  placed  on  conditions  encoun- 
tered in  low  pressure  discharges,  as  well  as  those  encoun- 
tered at  high  pressures  when  using  an  electron  beam- 
sustainer  configuration.  Both  of  these  discharge  cases 
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are  characterized  by  a low  degree  of  ionization,  - 10  , 

and  electron  mean  energies  in  the  range  of  0.5  to  2.0  ev. 

For  these  conditions,  the  electron  energy  distribution, 
eed,  is  non-Maxwe Ilian  and  exhibits  a structure  that  is 
consistent  with  the  features  of  the  electron-molecule 
excitation  processes.  Since  the  eed  is  generally  non- 
Maxwellian,  a detailed  analysis  of  the  electron  kinetics 
involving  a steady-state  solution  of  the  collisional  Boltz- 
mann equation  is  required  to  evaluate  the  electron  impact 
excitation,  ionization,  and  dissociation  rates.  Relatively 
accurate  and  complete  data  on  electron-impact  cross- 
sections  for  argon  and  CO,  have  been  reported.  These  are 
conveniently  listed  in  the  reference  compilation  of  Kieffer 
(ref  134).  Using  these  data,  the  eed  is  calculated  by 
numerically  integrating  the  Boltzmann  equation  assuming 
that  the  distribution  is  spatially  uniform  and  a steady- 
state  exists. 

Bo ltzmann  Equation  Theory . The  collisional  Boltzmann 
equation  is  the  fundamental  equation  governing  energy  trans- 
fer in  electron-molecule  collisions.  Several  detailed  and 
general  developments  have  previously  been  given  (refs  135, 
136,  137,  138,  139).  Here  we  will  merely  summarize  those 
features  relevant  to  the  present  analysis  and  provide  a con- 
venient review  of  the  development  following  the  previously 
cited  references.  The  eed  is  obtained  by  solving  the 
Boltzmann  equation 
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(6-1) 


where  f = f(v,~r,t)  is  the  electron  velocity  distribution 

/'gfN0 

function,  ’v  the  velocity,  'a  the  acceleration  and  (•g^J 
represents  the  rate  of  change  of  f due  to  electron  colli- 
sions of  all  types.  Note  that  the  terms  on  the  left-hand 
side  of  equation  (6-1)  represent  the  continuum  flow  in 
phase  space  of  the  eed,  whereas  the  right-hand  side 
accounts  for  discrete  transitions  in  the  phase  space  result- 
ing from  elastic  and  inelastic  collisions. 

For  typical  discharge  conditions  we  will  consider;  the 
electric  field,  gas  density,  and  electron  density  are  con- 
stant in  time  and  spatially  uniform.  By  introducing 
spherical  coordinates  in  the  velocity  space,  the  eed  can  be 
expressed  as  an  infinite  series  of  products 


4?,  v,  s)  = 


(6-2) 


where  Pn(0)  is  the  nth  Legendre  polynomial,  e the  polar 
angle  with  respect  to  the  electric  field  direction,  and 
f^  an  expansion  coefficient  which  is  a function  of  the 
electron  speed,  v,  and  time,  t. 

In  this  representation  ths  Boltzmann  equation  takes 
the  following  form  (ref  138) 
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with  y - cos  6.  If  we  next  introduce  the  following  identi- 
ties (ref  140) 


(/-X)J3M  = *&L  (a)  - rnjju)) 

(6-4) 
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equation  (6-3)  yields 
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By  equating  terms  with  the  same  angular  dependence,  an 
infinite  set  of  coupled  differential  equations  is  obtained. 
Under  normal  discharge  conditions,  the  electron  velocity 
distribution  is  nearly  spherically  symmetric,  and  retention 
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of  only  the  leading  terms  in  the  expansion  is  justified 
(refs  135,  136,  137).  Retaining  the  first  three  terms,  we 
obtain  (refs  136,  137,  138) 
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In  order  to  evaluate,  = {'It)  ++0f)  ' let  us  f^rst 

define  the  differential  scattering  cross-section  for  an 
electron  with  speed  v scattered  through  an  angle  a into 
a unit  solid  angle  to  be  q(v,a).  The  collision  geometry 
and  coordinate  definitions  and  relations  are  shown  in  Fig. 
6-1.  Note  that  41  is  the  azimuthal  angle  in  velocity  space 
around  the  vector  "v.  The  total  rate  of  scattering  out  of 
a differential  volume  in  velocity  space  at  v(v,9,<t>),  for 
process  h is  given  by  (refs  136,  138)  0^)  • where 


p Z*1  ft  qQ 

%&>*)  Y 4 ft  *">(«)  cfa  ctf  (6-9) 

Xdt  J o o * AsO 
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and  with  representing  the  gas  number  density  in  the 
initial  state.  The  effective  rate  is  obtained  by  summing 
over  all  processes,  h,  including  elastic  and  inelastic 
collision  channels. 


/■g  f \c+ 

Similarly,  the  positive  term,  , representing 

the  rate  of  scattering  into  the  volume  at  v from  u (u,6'  ,<J>') 


is  (refs  136,  138) 


zn  n ^ 

u j f frCfrO  (6-io) 

o o n-o 


where  J is  the  Jacobian  of  the  transformation  and  repre- 
sents the  ratio  of  that  differential  volume  in  velocity 
space  from  which  the  electron  scatters  to  that  volume  into 
which  the  electron  scatters, 

if  Sin  <x.  c/o<  J^Ju.  to  V^S/noC  Joe  dtc/v 

Consider,  for  illustrative  purposes,  a truly  "elastic" 
collision  in  which  there  is  no  energy  transfer;  but  only  a 
change  in  direction.  Under  these  circumstances 

u-v 

and  the  Jacobian  is  unity.  Then  upon  integrating  over  \p , 
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equations  (6-9)  and  (6-10)  yield  (ref  138) 


(-*¥+(■7$’  • (6‘u) 

' ' ' n*o 

where  the  elastic  cross-section  for  process  h and  order 
n is  given  by 


(6-12) 


l 


with  the  three  lowest  order  elastic  cross-sections  given  by 


Thus  under  the  "elastic"  assumption,  such .collisions  yield 
no  contribution  to  the  zeroth  order,  spherically  symmetric 
equation.  Q-^fv)  is  identified  as  the  experimentally  mea- 
sured momentum  transfer  cross-section,  Q^.,  and  represents 
the  first  non-zero  contribution. 

Now  consider  inelastic  collisions,  such  as  those 
resulting  in  vibrational  or  electronic  excitation.  In 
this  case,  one  must  properly  account  for  the  energy  loss. 
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From  equations  (6-9)  and  (6-10)  we  obtain  (refs  136,  138) 
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(6-13) 


(6-14) 


where 


•S/n  e<  c/<K 


is  the  cross-section  for  the  interaction  in  which  an 

electron  with  speed  v undergoes  process  h upon  colliding 

with  the  molecules  N^,  loses  energy  Ah,  and  terminates  the 

1 2 

collision  with  speed  u in  equation  (6-13)  with  j mv  = 

12  . . 

2 mu  + Equation  (6-14)  represents  the  positive  rate 

which  populates  the  element  at  velocity  "v.  Since  energy 
is  conserved  in  the  collision  process; 

U c/u  = Vc/V 

2 2 

and  the  volume  ratio  given  by  the  Jacobian,  u du / (v  dv) 
reduces  to  u/v.  Thus  for  inelastic  isotropic  scattering, 
the  lowest  order  collision  terms  are  (refs  136,  137,  138) 
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To  obtain  a form  more  suitable  for  calculation,  we 
transform  the  set  of  equations,  (6-6),  (6-7),  (6-8),  from 
functions  of  velocity  to  functions  of  energy,  w,  by  using 
the  following  relations; 
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We  then  obtain  (refs  136,  137,  138), 
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where  w'^+A^f  and  A^  is  again  the  energy  loss  associated 
with  the  inelastic  cross-section  Q^. 

Equations  (6-18,  6-19,  6-20)  represent  the  general 
form  of  the  coupled  differential  equations  used  for 
solving  the  collisional  Boltzmann  equation.  We  now  restrict 
ourselves  to  a two  term  expansion,  and  include,  (refs  135, 
136,  137,  138),  the  effects  of  energy  transfer  in  momentum 
transfer  collisions,  and  superelastic  collisions.  Assuming 
a quasi-steady-state  we  obtain  from  equations  (6-18,  6-19) 


® .6-2!, 

h 

where  is  the  relevant  excited  state  density  for  process  h. 

is  the  sum  of  the  momentum-transfer  cross-sections  for 
each  gas  weighted  by  mole  fraction  and  the  reciprocal  of 
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the  molecular  weight,  whereas  Qr  is  the  sum  of  the  momentum- 
transfer  cross-sections  weighted  simply  by  mole  fraction. 

The  ratio  of  the  electron  to  proton  mass  is  given  by  m/M. 

We  may  write  equations  (6-21,  6-22)  as  a single  equation 
in  fQ  by  integrating  equation  (6-21)  with  respect  to  energy, 
substituting  into  equation  (6-22)  , and  integrating  again 
to  yield  (ref  139) : 

-£<*)  - c £(£%)>■  (tt  £ Qt 


i 
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where  C is  determined  by  the  condition  f f (w)w,5dw=l  and 

o ° 

/*».  NhW'  ^ 

and  Qt=Qr+  2L.  Qh(w)  + ITw  Qh*w'y'  This  inte<?ral  form 

of  the  Boltzmann  equation  (ref  139)  can  be  solved  by 
iteration  for  fQ(w),  substitution  into  equation  (6-22) 
yields  f ^ (w) . 

Electron  Collision  Processes . Having  developed  the 
Boltzmann  formalism,  we  recognize  that  the  solution  requires 
cross-section  data  for  momentum  transfer,  vibrational  and 
electronic  excitation  of  the  gases  constituting  the  spe- 
cific mixture  under  consideration.  In  this  analysis  we 
will  be  concerned  with  mixtures  of  carbon  monoxide  and  argon. 
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Most  of  the  required  data  are  available  from  a JILA  compila- 
tion by  Kieffer  (ref  134).  These  experimental  data  have 
been  obtained  using  high  energy  resolution  electron  beam 
techniques  coupled  with  a detailed  analysis  (ref  141)  of 
electron  transport  data.  The  resulting  data  set  is  quite 
reliable  and  relatively  complete,  providing  an  adequate 
basis  to  quantify  our  results. 

The  momentum  transfer  cross-sections  for  argon  and 
carbon  monoxide  were  obtained  from  Kieffer  (ref  134) . 
Momentum  transfer  collisions  tend  to  reduce  the  anistropy 
in  the  eed  resulting  from  the  applied  field  and  channel 
some  of  the  discharge  power  into  "elastic  heating"  of  the 
gas . 

Vibrational  excitation  is  tne  dominant  excitation 
process  for  the  conditions  considered  in  this  study  because 
the  cross-sections  are  very  large  at  the  mean  energies 
characteristic  of  laser  discharges.  Vibrational  cross- 
sections  and  excitation  rates  are  required  for  an  analysis 
of  the  molecular  kinetics.  In  carbon  monoxide,  vibrational 
excitation  proceeds  predominately  via  "shape  resonances." 

The  cross-section  characterizing  this  excitation  displays  a 
structure  related  to  the  lifetime  of  the  compound  state, 

AB~;  and  the  characteristic  vibrational  frequency.  When  the 
lifetime  is  short,  as  it  is  in  H2 , the  cross-section 
exhibits  a broad  maximum.  In  CO,  the  lifetime  of  the  com- 
pound temporary  hegative  ion  is  comparable  to  the  vibra- 
tional period  and  the  cross-section  exhibits  broadened 
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spikes  corresponding  to  the  energy  resonances  of  the  com- 
pound state.  Erhardt  (ref  142)  and  Schulz  (ref  143)  have 
reported  electron  cross-section  data,  a , for  vibrational 

O /V 

excitation  from  the  ground  state,  v=0,  to  vibrational  level 

v for  v=l  to  8.  This  data  is  presented  in  Figs.  6-2  and 

6-3.  Referring  to  Fig.  6-2,  note  that  near  threshold  the 

effective  cross-section*  has  a value  of  approximately 
-17  2 

2 x 10  cm  which  rapidly  increases  to  a value  approaching 
-15  2 

10  cm  at  1.7ev  before  decreasing  by  orders  of  magnitude 
in  the  4ev  range.  Fig.  6-3  presents  the  individual  cross- 
section  data,  a , of  Schulz  for  excitation  from  the 

O f V 

ground  state  to  vibrational  state  v.  The  cross-sections 
all  exhibit  the  resonant  spikes,  and  decrease  in  magnitude 
with  increasing  v.  The  magnitude  and  energy  dependence  of 
these  cross-sections  will  be  reflected  in  the  calculated 
distribution  functions,  and  fractional  power  transfer, 
which  will  be  examined  later. 

Since  we  anticipate  highly  non-equilibrium  vibrational 
distributions  exhibiting  a substantial  excited  state  popula 
tion,  we  have  not  restricted  the  electron-molecule  vibra- 
tional excitation  to  excitation  only  from  the  ground  vibra- 
tional state.  Instead,  we  assume  a = o , for  m=l  to  8 

o,m  v,v+m 

Electron  excitation  cross-sections  from  excited  states  are 
not  available;  however,  guided  by  the  theoretical 

*The  effective  cross  section  is  defined  as  the  summa- 
tion of  all  vibrational  c.  oss-sections  from  the  ground 
vibrational  state. 
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calculations  of  Chen  (ref  144)  in  N2  and  based  on  the 
isoelectronic  relationship  of  N2  to  CO,  this  assumption 
is  reasonable.  Detailed  calculations  of  these  excited 
state  cross-sections  (ref  145)  and  experimental  measure- 
ments (refs  146,  147)  are  the  topic  of  current  research. 

Electronic  excitation  cross-sections  for  the  excita- 
tion of  a^II,  ALn,  b^E,  C^T,  E^JI,  and  a 13.5  ev.  loss 
(ref  148)  were  included,  as  well  as  the  dissociative 
attachment,  ionization  and  dissociative  ionization  cross- 
sections.  Under  conditions  typical  of  self-sustained  or 
e-beam  sustained  operation,  only  a small  fraction  of  the 
i energy  is  channeled  into  these  electronic  states.  However, 

it  is  necessary  to  include  these  states  in  order  to 
properly  evaluate  the  rate  of  dissociation  due  to  electron 
impact.  Consistent  with  our  reasoning  concerning  the 
highly  non-equilibrium  distribution  of  vibrational  states, 
we  have  investigated  the  implication  of  assuming  that  the 
dissociative  ionization,  DI,  cross-section  from  excited 
vibrational  states  is  equal  to  that  from  the  ground  state 
with  the  threshold  displaced  by  the  vibrational  energy  rela- 
tive to  v=0.  We  will  expand  on  this  point  later  when  we 
discuss  and  compare  dissociation  rates.  The  dissociative 
ionization  cross-section  measured  by  Rapp  (ref  149)  is  given 
in  Fig.  6-4.  A threshold  energy  of  16.5  ev  should  be  noted 
as  well  as  the  small  magnitude  in  comparison  with  the  vibra- 
tional cross-sections  given  previously.  In  addition  to 
dissociative  ionization,  dissociative  attachment,  DA,  occurs 
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and  will  be  included  in  this  analysis.  Fig.  6-5  displays 

the  cross-section  for  the  dissociative  attachment  process 

as  measured  by  Chantry  {ref  150) . This  cross-section,  in 

contrast  to  that  of  discociative  ionization,  is  resonant 

-19  2 

in  character  exhibiting  a pronounced  maximum  of  2 x 10  cm 
at  lOev.  Also  note  that  the  threshold  for  this  excitation 
is  approximately  7 ev.  below  that  of  dissociative  ioniza- 
tion. This  fact  coupled  with  the  relative  magnitudes  of 
the  two  cross-sections  results  in  dissociative  attachment 
being  the  dominant  electron  impact  dissociation  process. 

Calculated  Electron  Energy  Distributions . Using  the 
cross-section  data  and  the  integral  form  of  the  Boltzmann 
equation  previously  described,  the  electron  energy  dis- 
tribution was  calculated  for  a 10/90  mix  of  CO /Argon, 

(mole  fraction) , for  a wide  range  of  values  of  the  ratio 
of  the  electric  field  strength  to  the  total  neutral  gas 
number  density,  E/N.  The  selected  values  of  E/N  yield 
mean  electron  energies  in  the  range  of  .5-3.0ev.  These 
mean  energies  are  typically  encountered  in  self-sustained 
and  e-beam  sustained  discharges.  Fig.  6-6  is  a semi-log 

plot  of  the  calculated  energy  distribution,  with  the  dis- 

% 

tribution  density  at  energy  w given  by  fQ(w)w  dw.  In  this 
representation,  a Maxwellian  distribution  would  appear  as 
straight  line  with  a negative  slope  of  (l/kTe) , with  Tg 
being  the  "electron  temperature."  The  calculated  distribu- 
tions depart  significantly  from  a Maxwellian.  The 
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distributions  reflect  the  dominant  electron-molecule  col- 


lision  processes,  by  exhibiting  a rapid  decrease  in  the 
energy  region,  l-2ev,  of  significant  vibrational  excitation. 
The  large  resonant  vibrational  cross-sections  act  as  a bar- 
rier to  the  electrons , and  prevent  many  of  them  from 
obtaining  energies  in  excess  of  2ev.  As  the  field  strength 
is  increased,  more  electrons  are  able  to  penetrate  this 
barrier  and  the  distribution  develops  a high  energy  tail. 

The  dominance  of  the  vibrational  excitation  is  con- 
firmed by  examination  of  the  fractional  electron  power 
transfer  curve  presented  in  Fig.  6-7  for  the  conditions 
of  Fig.  6-6.  For  values  of  mean  electron  energy  below 
1.5ev,  E/N=20  Townsend,  vibrational  excitation  is  the 
dominant  energy  transfer  channel.  Electronic  excitation 
becomes  significant  only  at  higher  values  of  E/N,  and 
elastic  collision  losses  are  always  relatively  small.  The 
energy  loss  associated  with  dissociation  and  ionization  are 
completely  insignificant  in  comparison  with  the  other 
losses  in  the  E/N  range  considered. 

When  the  number  of  excited  molecules  increases,  an 
additional  energy  exchange  process,  a superelastic  colli- 
sion or  a collision  of  the  second  kind,  becomes  important. 
The  influence  of  superelastic  collisions  on  the  electron 


energy  distribution  can  be  understood  by  referring  to 


equation  (6-23) . The  bracketed  term  in  the  summation  is 
proportional  to  net  inelastic  excitation  rate.  Consider 


from  the  ground  state  and  the  reverse  process.  In  this 


I 

[ 


r 
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case  Nh  represents  the  number  density  in  the  state  v=0,  and 
N^,  the  number  density  in  v=l.  A significant  excited  state 
population  will  decrease  the  value  of  the  integrand,  reduce 
the  net  excitation  rate,  and  result  in  a lowering  of  the 
effective  vibrational  barrier.  As  a result,  the  tail  of 
the  electron  energy  distribution  will  be  raised  above  that 
value  obtained  when  superelastic  processes  are  omitted. 

This  effect  is  observed  in  Fig.  6-8  where  we  compare  the 
eed  with  and  without  superelastic  collisions.  The  vibra- 
tional distribution  assumed  in  this  calculation  was  of  the 
analytic  form  given  by  Gordiets  (ref  151) , which  will  be 
discussed  subsequently  in  the  vibrational  kinetics  section 
of  this  chapter.  The  0-1  vibrational  temperature  in  this 
calculation  was  2000°K  and  the  translational  temperature 
300 °K.  Note  that  the  degree  of  enhancement  in  the  eed 
tail  population  increases  with  decreasing  E/N,  and  that 
there  is  a smaller,  compensating  decrease  in  the  distribu- 
tion at  low  energies.  It  is  also  interesting,  yet  antici- 
pated, that  the  distributions  are  proportional,  or  self- 
similar, in  that  energy  region  beyond  the  vibrational 
exchanges.  Thus  vibrational  superelastic  processes  lead 
to  changes  in  the  electron  energy  distribution  function  and 
enhance  electronic  excitation  and  dissociation  rates. 

Summary.  We  have  examined  the  electron  kinetics  of  a 
CO-Ar  gas  discharge  within  the  framework  of  the  collisional 
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Boltzmann  equation.  For  this  purpose  we  developed  the 
integral  form  of  the  Boltzmann  equation,  specifically 
including  superelastic  collisions.  It  was  established  that 
the  electron  kinetics  are  dominated  by  the  large  resonant 
vibrational  cross-sections,  which  permit  efficient  molecular 
excitation  to  be  achieved.  The  dominance  of  the  vibra- 
tional cross-sections  was  reflected  in  the  calculated  eed's 
which  were  highly  non-Maxwellian.  Due  to  the  non-Maxwellian 
nature  of  the  distribution,  the  possibility  of  using  the 
Maxwellian  function  as  an  approximation  for  calculating 
transfer  rates  is  precluded.  Superelastic  collisions  were 
considered  and  shown  to  have  a significant  effect  on  the 
eed  leading  to  enhanced  electronic  and  dissociative  excita- 
tion rates.  The  individual  electron-molecule  excitation  and 
dissociation  rates  were  computed  for  a wide  range  of  E/N  values 
and  gas  temperatures.  These  rates  are  used  in  a following 
section  where  we  parameterize  the  dissociation  process 
under  non-equilibrium  conditions. 

Molecular  Kinetics 

Background.  Problems  of  molecular  and  vibrational 
kinetics  play  a central  role  in  the  fields  of  chemical 
kinetics  and  laser  physics.  Currently,  interest  is  centered 
around  an  analysis  of  energy  exchange  processes  occurring 
when  significant  amounts  of  energy  are  introduced  into  a 
gas.  All  interactions  be  they  thermal,  chemical,  optical, 
or  electrical  result  in  a redistribution  of  state 
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populations.  The  ultimate  distribution  then  determines 
the  nature  of  chemical  and  physical  processes  occurring  in 
the  resulting  non-equilibrium  state.  Here  we  will  be  pri- 
marily concerned  with  the  investigation  of  the  non- 
equilibrium vibrational  distributions  of  a molecular  gas 
excited  by  an  electric  discharge  or  obtained  through 
intense  optical  pumping  of  the  vibrational  states.  For 
the  moment,  we  will  conveniently  assume  the  existence  of 
some  "pump"  term  delivering  energy  to  the  vibrational  mani- 
fold and  proceed  to  characterize  the  molecular  kinetics 
and  the  vibrational  distribution.  Within  this  section  we 
first  review  the  fundamental  kinetic  processes  involved, 
their  relative  strengths  and  characteristic  time  scales. 
Based  on  a time  scale  argument  and  guided  by  the  results 
of  detailed  numerical  calculations,  we  will  establish  that 
there  are  3 characteristic  regions  of  the  vibrational 
energy  distribution,  VED.  A simplified,  analytic  approach 
to  the  calculation  of  vibrational  distributions  (refs  151, 
152)  under  highly  non-equilibrium  conditions  will  then  be 
reviewed.  The  dependence  of  the  VED  on  input  energy,  V-V 
and  V-T  rate  magnitude,  or  translational  temperature  will 
be  examined  using  this  approach.  Finally,  for  completeness 
and  for  the  additional  insight  they  reveal,  selected 
detailed  numerical  solutions  of  the  VED,  both  time  depen- 
dent and  steady-state,  will  be  presented.  These  distribu- 
tions will  be  compared  with  experiment  to  establish  the 
degrees  of  non-equilibrium  that  can  be  achieved  and  lead 


us  into  our  parameterization  of  non-equilibrium  dissocia- 
tion. 


Characteristic  Regions  of  the  VED.  Let  us  examine 
the  factors  that  determine  the  populations,  Nv,  of  the 
vibrational  energy  levels.  In  general  these  factors  can 
be  classified  as  follows. 

1.  Net  rate  of  external  excitation  or  de-excitation 
for  each  level  (pumping),  Nv 

pump 

2.  Rates  of  collisional  V-V  energy  transfer,  N 


3.  Rates  of  collisional  V-T  energy  transfer,  N 


We  will  omit  consideration  of  radiative  processes  as  loss 
mechanisms,  while  recognizing  the  limitation  they  impose 
when  considering  optical  pumping.  This  is  appropriate  for 
the  high  pressures  considered  in  this  study.  The  rate 
equation  for  each  vibrational  level  v,  from  equation  (3-5)  , 
can  now  be  written  in  the  simple  form: 


pump 


(6-24) 


By  examining  the  relative  strengths  and  time  scales 
of  the  three  terms  in  equation  (6-10) , we  will  establish 
that,  for  highly  non-equilibrium  conditions,  the  VED 
exhibits  3 characteristic  regions. 
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The  mechanism  for  vibrational  excitation  in  electric 
discharge  systems  is  electron  impact  collisions,  which  pro- 
ceed via  the  formation  of  temporary  compound  excited  states. 
As  noted  in  the  previous  section,  the  rate  of  excitation 
depends  not  only  on  the  eed,  but  also  on  the  population  of 
the  levels.  The  pump  term,  when  considering  optical 
excitation,  is  also  a function  of  the  level  population  of 
the  pump  transition  as  well  as  the  intensity  of  the  radia- 
tion. 

Subjecting  a diatomic  gas,  initially  at  a temperature 
T,  to  external  excitations,  results  in  a departure  of  the 
vibrational  distribution  from  its  initial  Boltzmann  dis- 
tribution given  by: 


(6-25) 


The  resulting  non-equilibrium  distribution  and  the  relaxa- 
tion process  will  be  governed  by  the  Master  Equation,  and 
depend  upon  the  various  energy  exchange  rates,  the  molecu- 
lar structure,  and  the  strength  of  external  excitation. 

We  are  interested  in  the  quasi-equilibrium  distribution  of 
Nv  under  conditions  of  intense  pumping. 

Representative  probabilities  of  the  various  exchange 
processes,  given  in  Fig.  6-9,  show  that  V-V  collisions 
dominate  the  relaxation  process  for  low  vibrational 
quantum  numbers.  The  V-V  exchanges,  which  involve  colli- 
sions between  excited  molecules  and  are  a near  resonant 
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process,  result  in  a rapid  redistribution  of  the  externally 
pumped  vibrational  quanta  prior  to  their  ultimate  deactiva- 
tion by  the  slower  V-T  process.  Guided  by  the  data  in 
Fig.  6-9  at  300°K,  let  us  assume  that,  for  sufficiently 
low  vibrational  quantum  number,  we  can  neglect  VT  exchange 
and  assume  rotational-translational  equilibrium  at  tempera- 
ture, T.  Under  these  conditions,  the  vibrational  distribu- 
tion is  given  by  the  Treanor  distribution  (ref  153) , 


K ‘H*  ev> 


(T)  exp  (tv) 

£1 


(6-26) 


In  order  to  gain  insight  as  to  the  meaning  of  <J>,  we 
define  the  population  of  the  v-th  level  by  the  relation: 


/i  = M0exp(- 


For  a Boltzmann  distribution. 


~7^-  constant  VmO)/i>  • -jV 


(6-27) 


Using  this  definition,  we  can  rewrite  4>  in  terms  of  , 
the  effective  1-0  temperature. 


j (CrE. o)  (T,-  •£>) 

r-  AT  AT, 


(6-28) 


I 


The  Treanor  distribution  becomes  upon  substitution  for  <f> , 


//K->V0  exp 


v ( 5-E0) 

*T, 


4- 


A r 


(6-29) 


Thus  we  can  interpret  <p  as  representing  the  influence  of 
an  external  pumping  of  the  0 -*■  1 transition,  resulting  in 
the  vibrational  temperature  T^. 

For  a simple  harmonic  oscillator,  Ev-Eq=v (E^-Eq) , and 
consequently  the  Treanor  distribution  degenerates  into  a 
simple  Boltzmann  characterized  by  a temperature  T^  For 
an  anharmonic  oscillator,  however,  v(E^-Eq)  > (Ev-Eq) . 

For  T^/T  > 1,  the  second  term  of  the  exponential  in  equa- 
tion (6-29) , now  leads  to  a distribution  with  an  increased 
population  relative  to  the  Boltzmann  at  T-^.  The  dependence 
of  the  Treanor  distribution  on  T and  T^  is  presented  in 
Figs.  6-10  and  6-11.  The  distributions  pass  through  a 
minimum  and  exhibit  a population  inversion  for  values  of 
T^/T  > 1.  Such  significant  inversions  are  not  obtained  when 
the  Master  Equation  is  solved  in  detail  with  V-T  losses 
being  included.  The  V-T  losses  prevent  these  large  total 
inversions.  Instead,  the  dominance  of  these  V-T  processes 
in  the  high  vibrational  levels  results  in  the  establishment 
of  a Boltzmann  "tail,"  at  a temperature  T,  for  the  popula- 
tion of  high  vibrational  levels.  This  behavior  can  be 
confirmed  by  noting  in  Equation  (3-13)  that  if  V-T  terms 
were  completely  dominant  the  entire  distribution  would  be 
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Fig  6-10  Treanor  Distribution  Parameterized  in 
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described  by  a Boltzmann  distribution  at  the  gas  tempera- 


ture. 

Numerical  solutions  of  the  complete  Master  Equation 
are  presented  in  Figs.  6-12  and  6-13,  where  effective  pump- 
ing rates,  , are  shown.  Three  regions  of  the  vibrational 
distribution  can  be  identified.  At  low  vibrational  levels, 
where  the  V-T  probabilities  are  insignificant  relative  to 
the  V-V  probabilities,  the  distribution  is  Treanor.  Beyond 
the  minimum  of  the  Treanor  distribution,  the  "resonance” 
or  "plateaxx"  region  is  established.  In  this  region  the  dis- 
tribution is  well  approximated  by  Nvv=constant . The  third 
region,  existing  at  high  vibrational  quantum  numbers,  is 
a Boltzmann  at  the  gas  temperature,  reflecting  the 
dominance  of  the  V-T  process.  The  existence  of  the  Treanor 
and  the  Boltzmann  regions  is  consistent  with  our  earlier 
discussion.  We  next  discuss  the  "resonance"  or  "plateau" 
region . 

The  vibrational  quanta,  being  pumped  into  the  gas  pre- 
dominately at  low  vibrational  levels,  are  deactivated  by 
local  V-T  collisions  when  the  pumping  is  weak.  However, 
when  strong  pumping  is  imposed,  the  local  V-T  deactivations 
are  insufficient  to  balance  the  source  and  the  distribution 
reacts  to  establish  a quanta  flow  balance.  This  balance  is 
achieved  through  the  action  of  the  rapid  V-V  exchanges, 
which  due  to  the  anharmonic  nature  of  the  energy  levels , 
preferentially  direct  the  excess  quanta  to  the  higher 
vibrational  levels  where  they  can  be  de-activated  by  the 
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Fig  6-12  Numerical  Solutions  of  VED  Parameterized  in 
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faster  V-T  collisions.  A "resonance"  or  "plateau"  region 
of  the  VED  results.  The  "resonance"  region  consequently 
exists  only  at  high  levels  of  pumping  and  low  translational 
temperatures,  and  serves  merely  as  a "connecting  link" 
between  the  pump  and  the  V-T  losses.  The  characteristic 
shape  of  the  distribution  can  be  explained  (refs  151,  152, 
157)  by  noting  that  for  the  stated  conditions,  resonant  V-V 
collisions  will  dominate  the  transfer  of  vibrational  quanta 
in  this  region.  Thus  the  constant  quanta  flux,  F,  supplied 
by  the  pump,  is  delivered  to  the  high  levels  by  resonant 
V-V  transfer,  with  the  flux  equation  given  by 


'v « L K 
/■' 


i'i 


Nv 

Y+t  v 


Nj’  VvP‘jrtjNy+,]sConsUnt  (6-30) 


An  approximate  evaluation  of  the  sum  (refs  151,  152,  157) 
yields,  within  the  SSH  formalism: 


/V^  V—  constant 


(6-31) 


> which  is  consistent  with  our  numerical  results. 

Analytic  Approach  to  VED1 s.  Having  reviewed  the  funda- 
mental kinetic  processes,  and  established  the  existence  of 
the  three  characteristics  regions  of  the  VED,  we  now 
highlight  the  development  of  an  analytic  expression  for 
the  vibrational  distribution  under  strongly-pumped 
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conditions.  The  details  of  this  development  are  given  in 
a series  of  papers  by  Gordiets  (refs  151,  154,  155).  An 
independent  and  more  general  approach  by  Lam  (ref  152) 
yields  a similar  form.  Here  we  follow  Gordiets  (ref  151) , 
and  start  out  with  the  Master  Equation  for  the  vibrational 
populations,  equation  (3-5);  with  the  notation  given  by 


- /v-m/ 


) 


(*-v) 


VT^  = /?o  Cm)eJ'"v 
E0v 0- <S0/+O) 

aE~  <fE0 


Gordiets  used  the  diffusion  approximation  to  obtain  an  ana- 
lytic form  for  the  populations  of  the  vibrational  levels. 
Omitting  here  the  intermediate  calculations,  we  reproduce 
his  final  result  for  Nv: 

(6-33a) 

(6-33b) 


_ (WaEJJ  V<£ 
- K*S”V , 

v t '*!  Q^KJvtCvh) 
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with  the  following  additional  notation  introduced 


c = Gt+i)  N,e»f>  f/ e*p(<r„  yj 


U - r,  a.  f 77  Y2-  Oto  ~1'36  &.E 
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(6-35) 


The  quantity  is  again  the  vibrational  "temperature"  of 
the  first  level;  v and  v are  respectively  determined  from 
the  conditions  for  the  minimum  of  function,  equation  (6-33a), 
and  for  the  vanishing  of  the  function,  equation  (6-33b) . 
where. 


I/  - 5T- 
2AE  X 


ex/°(JirT?)=  C<TvtQ,0K/F?q 


Let  us  now  examine  the  form  of  this  analytic  solution 
and  correlate  its  predictions  of  the  vibrational  population 
with  the  three  characteristic  regions  of  the  distribution 
previously  identified.  We  will  then  examine  the  influence 
of  the  V-V  and  V-T  rate  magnitude,  gas  temperature,  and 
pump  rate  on  the  calculated  distribution. 
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Referring  back  to  equation  (6-29)  and  taking  Eq  as 
the  zero  energy  reference  observe  that 

* v£Q  cT(V-/)  - VAE(Y-t)  (6-36) 

Inserting  this  identity  into  equation  (2-29)  yields 
equation  (6-33a) . Thus  the  analytic  distribution  for  v<v 
is  merely  the  Treanor  representation  of  the  distribution 
characterized  by  T.^  and  T.  This  is  reasonable  in  light  of 
our  previous  discussion.  Beyond  the  Treanor  minimum, 
v>v,  we  note  that 

A fy(V+0  ~ constant  (6-37) 

reflecting  the  characteristic  "resonance"  region  of  the 
distribution.  Thus  the  Gordiets  solution  explicitly  con- 
tains two  of  three  characteristic  VED  regions.  The  Boltz- 
mann "tail"  is  not  a limiting  form  of  the  analytic  solu- 
tion; but  the  onset  of  the  tail  correlates  well  with  the 
value  of  v.  Therefore,  by  limiting  the  falloff  of  the  dis- 
tribution to  that  of  the  Boltzmann  at  temperature  T,  reason- 
able agreement  with  detailed  numerical  solutions  are 
achieved.  Fig.  6-14. 

When  we  consider  the  dissociation  process,  we  will  be 
specifically  interested  in  the  exact  population  of  the 
Boltzmann  "tail"  region  and  can  not  rely  on  this  simple 
analytic  model.  However,  the  formalism  of  equation  (6-33) 
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is  very  useful  in  showing  explicitly  the  effects  of  the 
strength  of  pumping,  the  molecular  anharmonicity , and  the 
V-V  and  V-T  rates.  In  addition,  it  offers  a method  for 
quantitative  determination  of  kinetic  rate  information  for 
highly  excited  vibrational  states  (refs  133,  152). 

Referring  to  equation  (6-33) , let  us  examine  the 
explicit  effects  of  variations  in  different  parameters. 
First,  let  us  consider  variations  in  the  strength  of  pump- 
ing at  a fixed  temperature,  T.  An  increase  in  is  equiva- 
lent to  an  increase  in  the  net  pumping.  The  strength  of 
pumping  determines  the  location  of  v,  the  Treanor  minimum, 
and  subsequently  the  value  of  v required  to  balance  this 
source  term  with  a V-T  loss.  Note  that  as  T^  increases, 
v decreases,  being  proportional  to  T^  1.  This  results  in 
a significant  increase  in  Nv»  Thus  as  T^  increases  the 
location  of  the  Treanor  minimum,  v,  moves  to  a lower  vibra- 
tional quantum  number  and  achieves  a larger  population. 

Since  this  population  is  given  by  equation  (6-33b)  as 
C/ (v+1) , C also  increases.  Noting  that  the  second  term  in 
equation  (6-33b)  is  independent  of  T^,  we  see  that  v must 
increase  as  T^  increases.  This  behavior  is  confirmed  in 
Fig.  6-12  for  the  strongly  pumped  conditions  of  T^=3000, 
2000,  1500 °K  at  T=200°K. 


The  dependence  of  the  VED  on  translational  temperature, 
T,  is  slightly  more  complex,  since  <5^,  6W,  Aw,  Q^Q, 

Q1Q,  and  are  all  functions  of  this  parameter.  Fig. 

6-13  presents  the  detailed  numerical  parameterization  of 
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the  VED  with  respect  to  T.  The  gross  behavior  observed  can 
be  understood,  once  again,  in  terms  of  the  variation  of  v. 
In  this  case,  note  that  v is  proportional  to  T.  Therefore, 
v decreases  as  T decreases,  resulting  again  in  an  increase 
in  Nv*  As  before,  since  Nv  increases,  v will  also  tend  to 
increase  as  T decreases.  To  establish  the  variation  of 
v with  T,  we  must  also  consider  the  temperature  dependence 
of  the  second  term  in  equation  (6-33b) . For  the  cases  con- 
sidered. 


Q,0«  * T 


p * exp  (A/T 1/3 ) T>3oo 

Tic  7&  T<Z3°°  ’* 


(See  Fig.  6-15  for  examples  of  this  temperature  dependence 
of  the  V-T  rates.)  The  temperature  dependence  of  these 
terms  would  therefore  reinforce  the  increase  of  v as  the 
temperature  decreases.  This  increase  in  v,  however,  is 
tempered  by  the  T dependence  of  6^  given  by  the  SSH 
theory.  Thus  the  increase  in  v tends  to  saturate  at  the 

lower  temperatures. 

* 

The  mixture  dependence  of  the  VED  can  be  appreciated 
by  noting  in  this  case  that  P1Q  and  Q10  of  equation  (6-33b) 


now  represent  the  averaged  V-T  and  V-V  strengths  using 
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the  appropriate  mole  fractions  of  diluent  and  diatomic 


gases.  Depending  on  the  variations  of  the  ratio 


P10e 


«VTV 


°10K 


the  knee  will  move  out  to  higher  energy  or  in  to  lower 
energies  as  the  mixture  is  varied.  This  sensitivity  of  the 
location  of  the  knee  to  the  rate  magnitude,  suggests  an 
experimental  method  for  quantitative  determination  of 
kinetic  rate  information  for  highly  excited  vibrational 
states  (refs  133,  152). 


V-V  and  V-T  Rates  at  Low  Temperatures.  Before  proceed- 
ing with  the  parametric  mapping  of  dissociation  under  non- 
equilibrium conditions,  we  present  the  details  of  the 
kinetic  rates  used  in  this  section  of  the  analysis,  and 
later  compare  the  calculated  VED's  with  those  inferred  by 
Rich  (ref  133)  in  his  optical  pimping  studies. 

In  this  section,  we  will  continue  to  utilize  the  theo- 
retical form  for  the  heavy  particle  rates  developed  within 
the  Schwartz-Slawsky-Herzfeld  (ref  20)  theory  as  modified  by 
Keck  and  Carrier  (ref  28) . The  rates  have  been  cast  in  a 
form  that  incorporates  the  theoretical  quantum  number  depen- 
dence and  temperature  dependence,  and  is  fitted  empirically 
to  available  experimental  data.  At  low  temperatures,  the 

V-T  rates  depart  from  the  Landau-Teller  form  of  In  (Px) 

-1/3 

proportional  T , as  observed  in  Fig.  6-15.  In  addition, 
long  range  force  contributions  to  the  V-V  rates  become 
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significant  and  are  included  utilizing  fits  to  the  calcu- 
lations of  Sharma  and  Brau  (ref  81) . 

The  V-T  rates  used  are  based  on  the  semi-empirical 
modification  of  the  SSH  theory  proposed  by  Keck  and  Carrier. 
This  modification  provides  for  a smooth  variation  of  rates 
from  states  with  low  vibrational  quantum  number,  experi- 
encing near-adiabatic  collisions,  to  high  vibrational 
states,  where  collisions  approach  the  impulsive  limit.  This 
modification  attempts  to  extend  the  first  order  perturba- 
tion theory  results  of  the  SSH  theory,  which  are  invalid  for 
transitions  occurring  with  high  probability.  The  expression 
used  for  the  V-T  transition  probabilities  was  given  pre- 
viously in  equation  (3-3). 

The  V-V  rate  model  incorporates  a sum  of  the  contribu- 
tions from  a short  range  force  (SSH)  and  a long  range  force 
contribution  which  is  most  effective  for  near  resonant  col- 
lisions. The  present  calculations  incorporate  a weighting 
factor  for  each  contribution,  which  permits  a matching  to 
experimental  data  over  the  temperature  range  of  interest. 
Following  Rockwood  (ref  89)  and  Rich  (ref  99) , the  expres- 
sions used  for  the  V-V  transition  probabilities  for  CO  are 
given  by  equation  (3-4)  for  v+l+v  and  w->w+l: 
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where  S corresponds  to  the  short  range  and  L the  long 
range  contribution.  The  short  range  contribution  is  func- 
tionally similar  to  that  encountered  in  the  V-T  model  and 
is  given  by 


and  ::or  CO 


c/-  a ooS9Q 
<9  * 4.45  X/o6 

SXt)-*/.2Sx/6'sT',>‘’ 


The  long  range  contribution  for  CO  self-collisicns , 
which  arises  for  molecules  with  permanent  dipole  moments, 
is  included  following  the  theory  of  Sharma-Brau  (ref  81) . 
Their  rate  expressions  are  algebraically  complex,  therefore, 
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we  have  retained  the  functional  form  of  the  rate  being  pro- 
portional to  the  dipole  matrix  elements  and  specified  the 
resonant  character  by  a Gaussian  function  of  the  energy 
defect  (ref  89) . This  yields 

l°(t)g  &-w)  q <v*  o o (6-39) 


with 


l°(r)  - J.sx/o"  T~& 


The  V-V  rates  are  explicitly  formulated  to  satisfy  detailed 
balance , so  that 


vv  • 


V+I)V 


V.V+/ 

VVw+>jW  exp 


(6-40) 


The  resulting  fit  of  the  combined  long  and  short  range  con- 
tribution is  depicted  in  Fig.  6-16  where  it  is  compared 
with  the  experimental  data  of  Hancock  and  Smith  (ref  33) , 
Stephenson  (ref  160) , Sackett  (ref  159,  and  Powell  (ref  156) . 
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Note  that  this  form  provides  a reasonable  fit  to  the  experi- 
mental data,  being  a factor  of  2.0  higher  than  the  data  of 
refs  (33,  160,  159)  and  a factor  of  3 lower  than  the  recent 
data  of  Fushiki  (ref  158) . The  selection  of  this  particular 
fit  will  be  discussed  when  we  next  compare  the  predicted 
VED's  with  Rich's  (ref  133)  inferred  populations  from 
optical  pumping  experiments. 

Comparison  with  Optically  Pimped  VED's.  From  a 
kinetic  standpoint,  optical  pumping  offers  a simple  and 
contaminant  free  method  of  analyzing  molecular  energy 
exchange  processes.  Here  we  utilize  the  experimental  data 
of  Rich  to  benchmark  our  non-equilibrium  analysis.  Rich 
(ref  133)  used  a flowing  gas  absorption  cell  containing 
CO/Ar/He,  in  which  vibrational  excitation  was  obtained  by 
direct  optical  pumping  from  an  electrically  excited,  super- 
sonic flow,  CW,  carbon  monoxide  laser.  Utilizing  this 
technique,  mean  vibrational  energies  of  0.3  to  0.4  eV  per 
CO  molecule  were  obtained  at  a translational  temperature 
of  300 °K.  Vibrational  populations  were  inferred  from 
infrared  spectra  from  the  first  overtone, (Av=2) , emission 
bands  of  CO.  These  populations  are  presented  in  Fig.  6-17. 
Note  that  the  "resonance"  region,  in  the  zero  helium  mix, 
extends  out  to  vibrational  quantum  number,  v=40.  Addition 


of  helium  causes  the  knee  to  move  down  to  lower  quantum 
number,  reflecting  the  dominance  of  the  helium  V-T  rate 
over  that  of  argon  or  of  CO.  This  behavior  is  consistent 
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A comparison  of  measured  and  calculated 
vibrational  distributions  in  CO  mixes 
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with  our  previous  discussion  of  Gordiets'  analytic  form, 
6vtv 

since  P1Qe  ^®10K^  increases  with  helium  addition.  Also 

observe  that  the  remaining  lower  state  populations  are 
essentially  unchanged. 

Extracting  the  vibrational  temperature,  T^,  from 
Rich's  data,  a steady-state  calculation  of  the  VED  was  made 
for  T1=2700°K  and  T=300°K,  with  the  results  given  by  the 
dashed  line  in  Fig.  6-16.  Agreement  is  excellent,  the 
theoretical  curve  being  displaced  only  to  permit  convenient 
comparison.  Initial  calculations  using  V-V  rates  a factor 
of  2.25  lower  which  would  have  yielded  a better  fit  to  the 
experimental  data  of  ref  (33)  in  Fig.  6-16,  failed  to 
accurately  reproduce  the  data  and  resulted  in  the  distribu- 
tion for  the  0.6/94.4/5  mix  essentially  overlapping  the 
experimental  distribution  for  the  0.6/87.4/12  mix.  Based 
on  numerical  experiments,  it  was  established  that  agreement 
could  be  obtained  by  either  decreasing  the  magnitude  of 
CO-He  V-T  rate  or  increasing  the  magnitude  of  the  V-V  rate, 
while  retaining  the  same  vibrational  quantum  number  depen- 
dence. (Note  that  this  adjustment  confirms  the  dependence 

^VTV 

of  the  knee  location  on  the  ratio  P1Qe  /(Q1QK).  We  have 

selected  the  latter  adjustment  placing  more  confidence  in 
the  V-T  rate  measurements  and  noting  the  scatter  in  the 
measured  V-V  rates.  Fixing  the  value  of  the  V-V  rates  at 
this  value  in  all  subsequent  calculations,  the  zero  helium 
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mix  was  investigated.  Low  temperature  data  on  CO-Ar  V-T 
probabilities,  in  contrast  to  CO-He  and  C0-H2,  do  not 
exist;  we  therefore  extrapolated  the  high  temperature  data 
using  the  empirical  Pt  fit  of  Millikan  -and  White  (ref  5) . 
Using  this  form,  curve  1'  in  Pig.  6-l(/  was  obtained.  By 
increasing  the  magnitude  of  the  CO-Ar  V-T  rate  by  a factor 
of  5,  which  is  consistent  with  the  low  temperature  behavior 
of  CO-He  V-T  exchanges , good  agreement  with  experiment  could 
be  obtained.  The  result  of  this  adjustment  is  given  by 
curve  1 . 

Parameterization  of  Dissociation 

Having  discussed  the  electron  and  vibrational  kinetics 
and  benchmarked  the  calculational  techniques  and  kinetic 
rates,  we  are  now  prepared  to  investigate  the  dissociation 
process  under  highly  non-equilibrium  conditions.  In  this 
parameterization  we  will  compare  the  rate  of  dissociation 
resulting  from  heavy  particle  collisions  to  the  rate  of 
electron  impact  dissociation.  From  this  analysis,  we  will 
establish  two  characteristics  of  non-equilibrium  dissocia- 
tion. First,  non-equilibrium  conditions  lead  to  a signifi- 
cant rate  enhancement  of  both  types  of  dissociation. 

Second,  under  typical  conditions  encountered  in  self- 
sustained  and  e-beam  sustained  discharges,  electron  impact 
dissociation  exceeds,  by  orders  of  magnitude,  the  heavy 
particle  dissociation  rate. 
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A detailed  parametric  analysis  of  the  dissociation  pro- 
cess under  non-equilibrium  conditions  would  involve  a time- 
dependent  solution  of  the  coupled  Boltzmann  and  Master 
equations  for  electron  impact  dissociation,  or  similarly 
an  absorption  calculation  coupled  to  the  Master  Equation 
in  order  to  handle  the  case  of  optical  pumping.  Inherent 
in  the  analysis  would  be  an  extensive  parameterization 
with  respect  to  E/N,  electron  density,  gas  temperature  and 
radiation  intensity.  The  codes  developed  in  this  study  are 
well  suited  for  such  an  approach.  However,  rather  than 
embarking  on  such  an  ambitious  course  for  an  initial  sur- 
vey, we  have  chosen  a simplified  approach  which  adequately 
treats  the  dominant  processes  of  non-equilibrium  dissocia- 
tion and  permits  an  evaluation  of  magnitude  and  functional 
dependencies  of  the  heavy  particle  and  electron  impact 
dissociation  rates.  In  this  approach,  we  decouple  the 
electron  and  molecular  kinetics.  As  in  the  equilibrium 
case  previously  addressed,  we  are  primarily  interested  in 
the  behavior  of  a quasi-steady-state.  On  the  time  scales 
of  interest,  both  the  eed  and  the  VED  achieve  such  a steady 
state.  The  vibrational  distribution  is  relatively  insensi- 
tive to  the  form  of  eed  depending  mainly  on  the  total  net 
vibrational  excitation  and  gas  temperature.  The  net  excita- 
tion rate  of  the  vibrational  manifold  is  insensitive  to 
the  form  of  the  vibrational  distribution  for  conditions 
considered  in  this  study.  For  this  reason  we  have  param- 
eterized the  heavy  particle  dissociation  rate  in  either  T^, 
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or  pump  rate,  and  gas  temperature.  The  eed,  however, 
does  reflect  the  state  and  distribution  of  vibrational 
excitation,  especially  when  electronic  and  dissociative 
collisions  are  examined.  To  account  for  this  interaction, 
while  retaining  an  uncoupled  approach,  we  have  bounded  the 
vibrational  temperature,  T^  for  typical  discharge  condi- 
tions and  incorporated  Gordiets'  (ref  151)  analytic  form  of 
the  vibrational  distribution  into  the  electron  kinetics 
solution. 

We  will  consider  a gas  mixture  of  CO/Ar,  10/90,  as  an 
illustrative  and  representative  mix.  The  electron  impact 
dissociation  process  will  be  parameterized  in  E/N,  T^ , and 
T.  We  will  consider  E/N's  ranging  from  3 to  40  Townsend  and 
gas  temperatures  ranging  from  100  to  3000 °K  for  ^=2  000 
and  3000 °K.  This  space  adequately  covers  the  normal  opera- 
ting parameters  of  gas  discharges.  In  parallel,  the  heavy 
particle  dissociation  process  will  be  examined  over  a 
similar  parameter  space.  In  addition,  we  will  examine  the 
heavy  particle  dissociation  process  for  conditions  encoun- 
tered in  optical  pumping  studies.  Here  we  will  consider 

constant  pump  rates  which  span  the  optical  pumping  regime, 

3 7 -1  -1 

10  to  10  quanta  molecule  sec  , for  gas  temperatures 

ranging  from  100  to  2000 °K. 


Electron  Impact  Dissociation, 
ation  occurs  in  gas  discharges  and 


Electron  impact  dissoci- 
initiates  a change  in 
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the  discharge  chemistry,  gas  composition,  and  electron  and 
molecular  kinetics.  Often  these  changes  lead  to.  an  eventual 
degradation  of  laser  discharge  efficiency  or  a transition  of 
the  glow  discharge  to  an  arc  (ref  1)  . 

In  this  section,  we  will  present  data  on  the  rate  of 
dissociation  due  to  electron  impact  collisions.  These  data 
were  obtained  from  a numerical  solution  of  the  Boltzmann 
equation.  We  will  examine  the  influence  of  dissociation 
from  excited  states  and  the  influence  of  superelastic  col- 

• P • 

lisions.  The  dissociation  coefficient  k^  is  related  to  the 
rate  of  electron  impact  dissociation,  R^,  by 

^ = kd  HUt  <«-<!> 

-3 

with  ng  = electron  number  density  (cm  ) 

-3 

Ntot  = total  neutral  gas  number  density  (cm  ) 

6 3 

k^  = dissociation  coefficient  (cm  /sec) 

Here  we  will  consider  two  different  electron  impact 
dissociation  reaction  channels: 

a.  dissociative  ionization 

b.  dissociative  attachment 

We  will  first  examine  the  rate  of  dissociative  ionization 
and  establish  that  this  reaction  rate  can  be  significantly 
enhanced  under  non-equilibrium  conditions.  The  rate  of  dis- 
sociative attachment  will  then  be  presented  and  compared 
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with  the  rate  of  dissociative  ionization.  The  influence  of 
superelastic  collisions  on  both  reaction  rates  will  be  dis- 
cussed . 

The  dissociative  ionization  cross-section,  presented 
in  Fig.  6-4  was  utilized  for  the  calculation  of  k^  (E/N, 
TjT^).  The  calculated  coefficient  of  electron  impact  dis- 
sociative ionization,  k1|,  is  presented  in  Fig.  6-18.  In 
order  to  evaluate  the  influence  of  excited  states  on  the 
reaction  rate,  we  have  assumed  that  the  cross-section  from 
excited  vibrational  states  is  equal  to  that  from  the  ground 
state  with  the  threshold  being  displaced  by  the  vibrational 
energy  relative  to  the  v=0  state.  The  excited  state  vibra- 
tional distributions  were  calculated  using  the  Gordiets' 

(ref  151)  formalism.  When  dissociation  is  assumed  to  occur 
only  from  the  ground  vibrational  state,  the  calculated  rate 
is  low  in  magnitude  and  exhibits  a strong  dependence  on  E/N. 

Dissociative  ionization  from  excited  states  results  in 
a rate  enhancement  of  many  orders  of  magnitude.  This 
enhancement  is  due  to  the  highly  non-equilibrium  distribu- 
tion of  both  the  VED  and  the  eed.  The  excited  state  dissoci- 
ation rate  is  dominated  by  the  contributions  from  the  vibra- 
tional states  in  the  high  energy  end  of  the  VED  plateau. 

This  dominance  can  be  understood  in  terms  of  the  population 
and  the  threshold  for  dissociative  ionization  of  the  high 
lying  vibrational  states.  These  states  obtain  a signifi- 
cant population  for  large  T^  and  low  T as  a result  of 
anharmonic  pumping  and  the  formation  of  the  "resonant 


plateau."  Recalling  that  Nvv~  constant  in  this  region,  we 
note  that  there  is  only  a slight  decrease  in  state  popula- 
tion from  the  Treanor  minimum,  v,  to  the  knee  location,  v. 
Since  the  plateau  can  be  maintained  over  a wide  range  of 
the  vibrational  manifold,  (over  5ev) , a significant  reduc- 
tion in  the  threshold  for  dissociative  ionization  is 
achieved  with  only  a minor  reduction  in  population.  The 
magnitude  of  the  excited  state  dissociation  rate  reflects 
the  sensitivity  of  the  process  to  excitation  threshold.  To 
ellucidate  the  dependence  of  the  dissociation  rate  on  the 
location  of  the  knee  in  the  VED,  we  have  calculated  the 
excited  state  dissociation  coefficient  for  values  of  v 
ranging  from  60  to  30.  The  values  obtained  are  given  in 
Pig.  6-18. 

The  fact  that  the  eed  is  highly  non-Maxwellian  magni- 
fies the  contribution  of  dissociation  from  excited  states. 
Observe  in  Fig.  6-6,  that  the  eed  exhibits  two  plateaus. 

These  plateaus  exist  in  energy  regions  where  there  are  no 
significant  inelastic  losses,  0.0  to  l.Oev.  and  2.5  to  6.0ev. 
The  plateaus  are  terminated  in  a region  of  significant 
inelastic  excitation  and  exhibit  a Maxwellian  behavior  in 
these  regions.  The  characteristic  temperature,  Tc,  in  these 
regions,  where 

/(it.)  = ex/o(-u/(JT))  (6-42) 
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is  much  less  than  T , where 

avg 

^ AT  “ mean  Gner£ y (6-43) 

This  results  in  an  enhanced  sensitivity  of  the  rates  to  the 
threshold  energy  of  the  excitation  and  points  out  the  inade- 
quacy of  any  Maxwellian  model  of  the  eed. 

Superelastic  collisions  also  enhance  the  dissociation 
rate,  with  this  effect  being  strongest  at  low  E/N.  The 
influence  of  superelastic  collisions  on  the  dissociative 
attachment,  DA,  and  dissociative  ionization,  DI,  rates  is 
presented  in  Fig.  19.  The  DI  curve  includes  dissociation 
from  excited  states.  The  source  of  rate  enhancement  due 
to  superelastic  collisions  and  the  self-similar  form  of  the 
eed  were  discussed  previously  in  reference  to  Fig  6-8. 

Dissociation  Due  to  Heavy  Particle  Collisions:  Electric 
Discharge  Conditions.*  We  now  examine  the  heavy  particle 
dissociation  rate  over  the  discharge  parameter  space  and 
compare  this  rate  with  the  electron  impact  rate.  We  will 
find  that  electron  impact  dissociation  substantially  exceeds 
the  heavy  particle  rate  for  typical  discharge  operating  con- 
ditions . 

Here  we  will  utilize  the  standard  ladder  model  for  the 
heavy  particle  analysis.  This  is  a convenient  and  reasonable 

*It  is  assumed  throughout  that  we  are  dealing  with  a 
C.W.  discharge  and  not  a pulse  discharge  or  afterglow. 
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approach,  since  at  the  low  temperatures  considered  here, 
the  various  heavy  particle  dissociation  models  approach 
the  same  limit  and  no  rate  variation  of  significance  in 
this  parameterization  results.  The  heavy  particle  coeffi- 
cients, k.  , are  parameterized  in  T and  T.  . With 

Q„  J. 


(6-44) 


The  results  of  this  analysis  are  presented  in  Figs.  6-20 

and  6-21.  Fig.  6-20  compares  the  electron,  k^,  and  heavy 

particle,  k.  , coefficients  for  conditions  of  vibrational 
aH 

equilibrium,  T^=T.  The  heavy  particle  rate  coefficient 

retains  its  Arrhenius  form,  and  is  very  low  in  comparison 

to  the  electron  impact  coefficient  at  the  typical  discharge 

operating  temperatures  below  700 °K.  The  influence  on  k. 

aH 

of  discharge  induced  departures  of  VED  from  equilibrium  is 

established  in  Fig.  6-21.  Here  k.  is  presented  as  a func- 

aH 

tion  of  T for  T1=2000°K  and  3000 °K.  The  departure  of  the 
VED  from  equilibrium  enhances  the  population  in  the  high 
energy  states  and  substantially  increases  the  dissociation 
rate  above  its  equilibrium  value.  Despite  this  increase, 
the  heavy  particle  dissociation  rate  remains  many  orders 
of  magnitude  below  k®  and  never  obtains  a physically  sig- 
nificant magnitude.  The  magnitude  and  temperature  depen- 
dence of  the  heavy  particle  rate  can  be  understood  in  terms 
of  the  simplified  analysis  of  the  VED  previously  presented. 
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E/N  (TDWN5END) 

Fig  6-20  Equilibrium  Comparison  of  Electron  and  Heavy  Particle  Dissociation 


Fig  6-21  Non-Equilibrium  Heavy  Particle  Dissociation  Constant 


Within  the  framework  of  the  ladder  model,  the  dissoci- 
ation rate  is  proportional  to  the  occupation  of  the  last 
bound  state.  The  occupation  of  the  last  bound  state  is 
substantially  increased  due  to  the  strong  influence  of  V-V 
exchanges  and  the  subsequent  anharmonic  pumping.  This  state 
population  may  be  approximated  as 

= Nt  **p[-nrrd]  ,6-45> 

where  we  have  neglected  the  slight  decrease  in  the  popula- 
tion experienced  in  the  plateau  region.  Nv  represents  the 
state  occupation  at  the  Treanor  minimum.  This  population, 
according  to  equation  (6-33)  increases  with  decreasing  gas 
temperature.  In  addition,  as  the  gas  temperature  is 
lowered  E-  increases  as  the  knee  moves  out  to  higher 
energies.  Thus  the  last  state  population  increases  as  the 
temperature  is  lowered,  but  this  behavior  does  not  persist. 
At  temperatures  below  350 °K,  these  population  enhancement 
effects  begin  to  saturate  and  the  inverse  temperature  depen- 
dence in  the  exponential  of  equation  (6-45)  dominates  the 
expression.  The  population  of  the  last  state  then  falls  off 
rapidly. 

Dissociation  Due  to  Heavy  Particle  Collisions:  Optical 
Pumping.  In  order  to  complete  the  parameterization  of  the 
heavy  particle  rate,  we  now  consider  the  case  in  which 
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excitation  is  achieved  through  intense  optical  pumping. 
Optical  pumping  has  advantages  over  discharge  excitation : 


1.  A greater  degree  of  vibrational  excitation  can  be 
achieved  in  comparison  to  electric  discharge  excitation. 

2.  Gas  mixture,  pressure,  and  temperature  need  not  be 
compatible  with  discharge  operation. 

These  advantages  force  one  to  consider  the  possibility  of 
achieving  a significant  heavy  particle  dissociation  rate 
under  conditions  of  intense  optical  pumping.  In  this 
section  we  parameterize  the  dissociation  reaction  under 
conditions  of  intense  optical  pumping  with  respect  to  gas 
temperature  and  effective,  v=0-*-l,  pump  rate.  The  tempera- 
tures considered  ranged  from  100  to  1000°K,  while  pump  rates 

extended  from  those  characteristic  of  discharge  excitation 

3 7 

to  gas  breakdown  threshold,  10  to  10  quanta/ (molecule- sec . ) 

Paralleling  our  treatment  of  thermal  dissociation,  we 
will  first  present  the  results  of  a time-dependent  solution 
of  the  Master  Equation,  Fig.  6-22.  This  solution  estab- 
lishes the  existence  of  a quasi-steady-state  and  reveals  the 
temporal  evolution  of  the  Treanor,  Plateau,  and  Boltzmann 
regions  of  the  VED.  We  here  consider  a 10/90,  CO/Ar  gas 

mixture  with  pump  rate  P=103  and  T=300  at  a total  number 

19  -3 

density  of  2.67  x 10  cm  . The  initial  distribution  is 
assumed  to  be  the  equilibrium  Maxwellian  distribution  at 
the  gas  temperature  T.  The  quasi-steady-state  is  achieved 
in  90  ysec  and  exhibits  the  three  characteristic  regions 
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discussed  previously.  In  the  early  stages  of  the  transient 
evolution,  the  relaxation  is  dominated  by  the  0+1  pumping 
and  off  rc-sonant  V-V  exchanges  with  the  first  excited  vibra- 
tional state.  These  exchanges  drive  the  low  vibrational 
states  toward  a Treanor  distribution  and  lead  to  a signifi- 
cant increase  in  the  population  of  the  higher  vibrational 
levels.  This  evolution  proceeds  until  a sufficient  popula- 
tion is  achieved  in  these  high  energy  states.  Then  vibra- 
tional exchanges  become  dominated  by  near  resonant  exchanges 
and  the  plateau  begins  to  be  established.  The  input  quanta 
are  conveyed  along  the  plateau  towards  the  Boltzmann  deacti- 
vation zone.  The  zone  moves  toward  higher  energies  to 
finally  achieve  a balance  with  the  energy  input  rate  and 
the  quasi-steady  state  is  established. 

A similar  VED  evolution  and  attainment  of  a quasi- 
steady-state was  obtained  for  all  pump  rates  and  gas  tempera 
tures  considered.  In  no  instance  was  it  possible  to  extend 
the  plateau  region  to  the  dissociation  limit  and  obtain  a 
VED  without  a Boltzmann  "tail."  As  a result,  the  heavy 
particle  dissociation  rates  obtained  for  conditions  of 
intense  optical  pumping  were  extremely  slow  and  not  com- 
parable to  those  obtained  by  electron  impact  excitation. 

The  heavy  particle  dissociation  coefficients  obtained 

3 6 —1 

for  pump  rates,  P,  of  10  and  10  are  plotted  against  T 

in  Fig.  6-23.  The  dissociation  coefficient  decreases  mono- 

tonically  with  temperatures,  falling  off  more  rapidly  at 

temperatures  below  400°K.  This  fall  off  is  due  to  the 
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Fig  6-23  Optical  Pumping,  Heavy  Particle  Dissociation  Constant 


exponential  temperature  dependence  of  the  final  state  popu- 
lation discussed  previously  with  regard  to  the  discharge 
parameterization  and  equation  (6-45) . The  gross  tempera- 
ture dependence  is  different  from  that  obtained  in  the  T^ 
parameterization.  This  can  be  understood  by  noting  that 
although  T^  is  proportional  to  a pump  rate,  and  an  increase 
in  T^  at  a fixed  T corresponds  to  an  increase  in  pumping, 

T^  being  equal  to  a constant  does  not  imply  a constant  pump 
rate  if  the  gas  temperature  is  varied.  Note  that  by 
increasing  the  pump  rate  by  three  orders  of  magnitude,  the 
dissociation  rate  has  increased  over  10  orders  of  magnitude. 
This  enhancement  is  caused  by  the  increase  in  plateau  popu- 
lation and  the  extension  of  the  "knee"  to  higher  energies 
as  the  pumping  is  increased.  This  behavior  saturates,  how- 
ever, as  shown  in  Fig.  24,  where  the  dissociation  constant 
is  plotted  as  a function  of  pump  rate  at  a fixed  temperature 
of  300°K. 

Evaluation  and  Comparison  of  Heavy  Particle  Dissocia- 
tion Rates.  Based  upon  the  discharge  and  optical  pumping 
analyses,  we  conclude  that,  within  the  framework  of  the 
ladder  model,  it  is  difficult  to  achieve  substantial  heavy 
particle  induced  collisional  dissociation  under  non- 
equilibrium conditions.  Quantitative  experimental  data  on 
non-equilibrium  dissociation  rates  are  not  available,  how- 
ever, qualitative  information  has  recently  been  obtained  by 
Rich  (ref  133)  in  his  optical  pumping  studies.  These 
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studies  were  motivated  by  interest  in  high  energy  CO  lasers 
and  infrared  avionics  programs,  and  were  directed  towards 
measurement  of  the  V-T  deactivation  rates  for  high  vibra- 
tional quantum  levels  of  CO  in  collisions  with  helium  and 
argon.  During  the  investigation,  however,  significant  non- 
equilibrium dissociation  of  CO  was  observed  to  occur,  the 
inferred  dissociation  rate  being  similar  in  magnitude  to 
that  obtained  under  electric  discharge  conditions.  This 
observation  is  inconsistent  with  the  predictions  of  our  non 
equilibrium  calculations  and  requires  a re-evaluation  of 
the  kinetic  model  used  in  the  analysis. 

First,  let  us  review  Rich's  observations  (ref  133). 

1.  Carbon  is  deposited  on  cell  walls  during  excita- 
tion. 

2.  At  gas  pressures  below  6 atmospheres  in  CO/Ar,  a 
visible  blue  emission,  (Swan  Bands  of  C 2) , is  observed.  Th 
intensity  increases  as  the  total  pressure  decreases. 

3.  A comparison  of  the  relative  intensities  of  the 

13  12  12  12 

C C v=l-+-0  and  the  C C v=l-*0  band  heads  suggest  an 
isotopic  enrichment  of  "several"  hundred  percent. 

4.  Observed  isotopic  enrichments  appear  to  be  essen- 
tially independent  of  gas  temperature. 

5.  Emission  from  the  CO  electronc  bands  is  not 
observed. 

Now  let  us  re-examine  our  model  predictions  of  non- 
equilibrium dissociation  rates  in  light  of  these  observa- 
tions and  the  calculated  VED's.  The  obvious  deficiency 


lies  in  the  magnitude  of  the  dissociation  constant.  For 

-43  -43  -6  -3  -1 

k.=10  , a dissociation  rate  of  10  N N =10  cm  sec 

d co  tot 

is  obtained.  This  rate  would  have  to  be  increased  about  20 
orders  of  magnitude  to  yield  a rate  consistent  with  experi- 
ment. 

A disparity  of  this  order  of  magnitude  might  result 
from  inaccuracies  in  the  kinetic  rates  used  in  the  calcula- 
tion. To  provide  a measure  of  the  sensitivity  of  the  cal- 
culation to  kinetic  rates  and  molecular  potential  and  obtain 
a perspective  as  to  what  variation  in  the  VED  would  be 
required  to  produce  such  a significant  variation  in  the  dis- 
sociation rate,  consider  the  following  two  cases.  First, 
assume  that  the  V-V  and  V-T  models  did  not  accurately 
identify  the  "knee"  location.  If  the  knee  actually  was 
displaced  by  1 ev  to  a higher  energy,  the  final  state  popu- 
lation and  rate  would  be  increased  by  e^®~  10^  at  300 °K,  or 
a factor  of  (55)^  for  every  tenth  of  an  ev.  displacement,  j. 
As  a second  situation,  and  really  a limiting  form  of  the 
first  case  mentioned,  consider  the  possibility  of  obtaining 
a VED  in  which  the  plateau  extends  to  the  dissociation 
limit.  Under  these  circumstances , neglecting  vibrational 
depletion  due  to  dissociation,  the  rate  would  approximately 
be  increased  by  a factor  g,  where 


It  is  therefore  critical,  in  a non-equilibrium  ladder 
climbing  analysis  of  dissociation,  to  know  the  magnitude 
and  level  scaling  of  the  V-V  and  V-T  rates.  This  data  is 
currently  not  available.  The  current  work  of  Rich  (ref 
133) , Lam  (ref  152)  , and  Fisher  (ref  161)  on  V-V  and  V-T 
rate  measurements  for  high  vibrational  states  will  be 
extremely  beneficial  in  this  respect. 

Before  setting  aside  the  delicate  question  of  rates, 
let  us  return  to  the  reality  of  experiment  and  examine 
Rich's  VED  data  as  previously  presented  in  Fig.  6-16. 

In  view  of  the  reasonable  agreement  obtained  between  theory 
and  experiment,  we  conclude  that  while  the  heavy  particle 
dissociation  rate  is  extremely  sensitive  to  the  magnitude 
and  level  scaling  of  the  rates,  the  rates  used  in  the 
calculation  are  adequate.  We,  therefore,  do  not  attribute 
the  disparity  in  the  dissociation  rate  to  inaccuracies  in 
the  V-V  and  V-T  rates. 

Continuing  now  with  an  examination  of  the  experimental 
data,  we  establish  that  the  "ladder"  model  itself  is  incon- 
sistent with  the  experimentally  observed  dissociation 
reaction  in  CO/Ar.  The  Boltzmann  "tail"  of  the  VED  of  experi- 
mental curve  onsets  at  less  than  3/4  of  the  11  eV.  dissoci- 
ation energy  of  CO.  The  approximate  population  of  the  last 
bound  state  would  be  10  ^ of  the  total  number  density 
yielding  a k.  of  approximately  10  6^.  This  magnitude, 
although  consistent  with  numerical  calculations,  is  not 
consistent  with  the  experimentally  observed  rate.  We 
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therefore  conclude  that  under  these  conditions  dissocia- 
tion does  not  only  proceed  from  vibrational  states  near  the 
dissociation  limit.  A more  efficient  reaction  channel  must 
exist. 


Alternative  Processes  and  Reaction  Mechanisms.  Let's 
examine  various  alternative  processes  and  reaction  mecha- 
nisms. First,  we  assume  from  the  description  of  the  experi- 
ment, ref  (133),  that  under  the  present  conditions  of 
optical  pumping  it  is  reasonable  to  neglect  any  electron 
impact  excitation  process.  We  can  also  neglect  any  uncer- 
tainties in  knee  location  since  the  experimental  distribu- 
tions were  inferred  from  the  radiative  properties  of  the 
molecule.  We  consider  next  the  possibility  of  a multi- 
photon, photo-dissociation  of  vibrationally  excited  CO. 


CO(v)  + n(h> ') 


C +o 


(6-46) 


The  pump  intensities  used  in  the  experiment  do  not  favor 
this  process,  especially  for  large  n~10,  (ref  171).  For 
n~3,  states  accessed  would  lie  in  the  Boltzmann  tail  and 
yield  a rate  much  lower  than  observed.  We  therefore  con- 
clude that  the  observed  dissociation  results  from  a heavy 
particle  collision  process.  For  the  dissociation  process 
to  be  efficient,  it  must  access  states  in  either  the  plateau 
or  Treanor  regions  of  the  VED.  We  suggest  two  possible 
reaction  mechanisms: 
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a.  Multi-quantum  V-V  exchange 

b.  Vibrational- to-Electronic  exchange 

For  mechanism  (a)  the  dissociation  reaction  proceeds  via, 


COCv)  + CO(v')  C + 0+  CO(v")  (6-47) 


Whereas  for  mechanism  (b) , the  dissociation  reaction  is 
related  to  the  following  reaction  sequence. 


co(v)  + CO(V')  — CO*  + COCv") 


(6-48A) 


CO*  + coCv)  -*•  CQ  + C 


(6-48b) 


where  CO*  is  an  electronically  excited  state  of  CO.  Both 
reactions  require  the  exchange  of  significant  amounts  of 
vibrational  energy,  with  |v-v'|~30.  Experimental  data  on 
multiple  quantum  exchanges  of  this  order  does  not  exist. 
However,  the  recent  data  of  Brechignac  (ref  162)  for  two 
quantum  exchange,  |v-v' | =2,  suggest  a large  cross-section, 
up  to  20%  of  that  of  the  single-quantum  V-V  transition. 
Theoretical  analyses  of  multiple  quantum  exchanges  are 
also  very  limited  and  complex.  These  analyses  do  not 
address  quanta  exchanges  of  the  order,  (Av=30) , considered 
here.  However,  the  theory  of  Dillion  and  Stephenson,  (ref 
163) , for  av=2  appears  to  qualitatively  agree  with  the  data 
of  Brechignac.  It  is  therefore  difficult  to  directly 
estimate  the  transition  rates  for  the  reactions  given  in 
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equations  (6-47)  and  (6-48).  We  can,  however,  estimate 
the  magnitude  of  the  transition  rate  that  would  be  required 
to  be  consistent  with  the  observed  dissociation  rate. 
Assuming  that  reaction  equation  (6-48a)  is  rate  control- 
ling, both  equations  (6-47)  and  (6-48a)  may  be  written  as 


«co  ^ 


w 


*u  H,„ 


(6-49) 


with 


g = 


diss 


N = 


the  specific  collision  frequency  (cm  / sec) 

the  exchange  probability 

the  vibrational  state  population 


Rich  has  stated  (ref  164)  that  the  observed  dissociation 
resulted  in  a variation  of  less  than  1%  in  the  CO  density  in 
a time  of  60  sec.  We  can,  therefore,  place  an  upper  bound 


-21  3 

on  the  dissociation  coefficient  k,  , at  10  cm  /sec. 

dH 


This  would  require  a maximum  value  of  for  both 


_ 3 _ c. 

mechanisms,  to  be  of  the  order  of  10  -10  . In  view  of 


Brechignac's  results  (ref  162),  we  consider  this  probabil- 
ity to  be  obtainable  and  continue  with  analysis  of  dissoci- 
ation via  mechanisms  (a)  and  (b) . 

In  discussing  non-equilibrium  dissociation  via  a multi- 
quantum V-V  exchange  or  a V-e  transition,  we  will  attempt 
to  characterize  dependence  of  each  mechanism  on  gas  tempera- 


ture and  T^.  In  addition,  we  will  examine  the  consistency 
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of  the  proposed  mechanisms,  a and  b,  with  experimental 
observations.  Finally,  simple  experiments  and  parameteriza- 
tions  will  be  suggested  to  establish  and  evaluate  the  rela- 
tive importance  of  each  mechanism. 

Both  mechanisms  involve  vibrational  states  in  either 
the  Treanor  or  Plateau  regions  of  the  VED.  Consequently, 
they  will  exhibit  a similar  dependence  on  and  T.  The 
reaction  rate,  if  primarily  proportional  to  the  level  occupa- 
tion, will  increase  with  increasing  pump  rate  and  increase 
with  decreasing  gas  temperature  at  a fixed  pump  rate.  This 
behavior  is  consistent  with  Rich's  observation  of  increased 
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the  isotope  in  this  binary  mixture  of  diatomics  (refs  18, 
153,  133) . Using  the  theory  of  Gordiets  (ref  18) , Rich 
(ref  133)  has  calculated  the  ratio  8 of  the  specific  dis- 
sociation rate  for  to  the  corresponding  rate  for 

12  16 

C 0 . These  results  (ref  133)  are  presented  in  Fig. 

6-25  where  the  separation  coefficient,  6,  is  plotted  as  a 

function  of  temperature  T,  for  various  values  of  T^.  The 

essence  of  the  solution  can  be  understood  by  noting  that 

the  calculation  assumes  that  the  VED  of  the  isotope  is 

12  16 

dominated  by  resonant  exchange  processes  with  CO,  the 

normal  isotopic  abundance  being  only  1.1%.  Since  the 

vibrational  energy  levels  of  the  isotope  are  more  closely 

12  16 

spaced  than  that  of  C 0 , this  resonant  exchange  occurs 

between  different  vibrational  levels  of  the  isotope  and 
12  16 

C 0 . The  resonance  shift  is  approximately  two  vibra- 

tional levels.  Thus  the  isotope  VED  essentially  reproduces 
that  of  the  with  a two  level  shift;  i.e. 


H,(cmo*)  _ />U  (c“dl) 

"m(c"oV~  "tot  (cacr) 


This  yields  a T^  for  the  isotope  greater  than  that  of 
12  16 

C 0 , and  results  in  a higher  relative  excited  state  popu- 

lation in  the  isotope  VED.  Note  that  higher  values  of  n 
in  Fig.  6-24,  corresponding  to  lower  T^,  give  a larger 
separation  coefficient,  but  yield  a lower  rate  magnitude. 
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TRANSLATIONAL  TEMPERATURE,T,°K 

ISOTOPE  SEPARATION  COEFFICIENT  FOR  LASER 
INDUCED  CO  DISSOCIATION 

Fig  6-25  Isotope  Separation  Coefficient  for  Laser 
Induced  CO  Dissociation  (ref  133) 


In  addition,  this  analysis  reveals  that  the  separation 
coefficient  is  fairly  insensitive  to  gas  temperature,  in 
quantitative  agreement  with  Rich's  observation  (ref  133). 

Based  on  the  previous  discussion,  both  a V-V  or  V-e 
exchange  mechanism  are  consistent  with  experimental  observa- 
tions. Both  mechanisms  exhibit  a similar  dependence  of 
the  dissociation  rate  and  isotopic  separation  coefficient 
on  and  T.  In  order  to  evaluate  the  relative  role  played 
by  each  mechanism,  the  following  experiments  are  suggested. 

Examination  of  reactions  (6-47)  and  (6—4 8b)  reveals 
a significant  difference  in  the  primary  dissociation  pro- 
ducts. The  V-V  mechanism  results  in  the  production  of  both 
atomic  carbon  and  oxygen,  while  the  V-e  exchange  mechanism 
produces  C02  and  atomic  carbon.  A mass  spectrometric  analy- 
sis of  the  gas  establishing  the  relative  of  amounts  of  0, 

02,  and  C02  would  provide  great  insight.  The  presence  of 
significant  amounts  of  either  0 or  02  would  favor  the  V-V 
mechanism,  whereas  their  absence  and  the  presence  of  C02 
would  tend  to  indicate  the  dominance  of  the  V-e  exchange 
mechanism. 

In  addition  to  the  difference  in  the  dissociation  pro- 
ducts of  mechanisms  (a)  and  (b) , we  note  that  the  activation 
energies  for  the  V-V  and  V-e  mechanism  also  differ  substan- 
tially. For  the  V-V  mechanism,  it  is  required  that 
Ev+Ev, >Dq,  11  eV.;  while  for  the  V-e  mechanism  Ev+Ev,>E*, 

E*  denoting  the  energy  of  the  electronic  state,  6.2  eV.  for 
3 

CO  (a  If)  . Capitalizing  on  this  difference,  a sequence  of 
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experiments  in  which  controlled  amounts  of  helium  are  intro- 
duced into  the  gas  mix,  while  the  Swan  band  intensity  is 
simultaneously  being  monitored,  would  establish  the  depen- 
dence of  the  dissociation  rate  on  the  location  of  the  knee 
in  the  VED.  This  information  could  then  be  used  in  conjunc- 
tion with  the  results  of  the  gas  mass  analysis  to  determine 
the  relationship  of  v and  v'  for  each  mechanism.  If  the 
V-V  process,  a,  proceeds  through  states  v and  v' , with 
v'»v,  then  a rapid  fall  off  in  Swan  band  intensity  would 
occur  only  when  the  knee  drops  below  5.5  ev.  If  the  process 
proceeds  dominantly  with  v'>v,  then  a more  continuous  reduc- 
tion in  intensity  will  be  observed  as  the  knee  moves  to 
lower  energies. 

It  was  previously  believed  that  the  V-e  transfer  process 
proceeded  dominantly  through  those  vibrational  levels  nearly 
resonant  with  the  electronic  state  (ref  168) . More  recent 
data  has  revealed  that  this  situation  does  not  always  hold 
(ref  169) . The  helium  addition  experiments  would  shed  light 
on  this  aspect  of  the  V-e  transfer.  If  a cutoff  of  the  dis- 
sociation reaction  is  observed  when  the  knee  drops  below 
6.2  ev. , a resonant  exchange  would  be  suggested.  A "softer" 
fall  off  of  the  dissociation  rate,  on  the  other  hand,  would 
favor  an  off  resonant  exchange. 

Further  work,  both  theoretical  and  experimental,  is 
required  in  the  area  of  V-e  and  V-V  exchanges.  However, 
based  on  Rich’s  experimental  data  and  our  analysis,  both 
mechanisms  offer  the  possibility  of  a more  efficient  channel 
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for  the  non-equilibrium  dissociation  of  diatomic  molecules 
and  provide  a simple  method  of  achieving  significant  isotope 
separation  factors  using  optical  pumping. 
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Chapter  VII . A Review  and  Some 
Considerations  for  Future  Study 


General  Review 


In  the  present  study,  we  have  performed  a detailed  and 
comprehensive  analysis  of  the  collisionally  induced  dissoci- 
ation of  diatomic  molecules  under  equilibrium  and  non- 
equilibrium conditions.  In  this  chapter  we  will  highlight 
the  interesting  and  important  results.  By  referring  again 
to  Table  II,  one  can  establish  the  scope  of  previous  works 
and  note  the  contradictions  that  exist  in  these  studies. 

This  study  has  resolved  many  of  these  inconsistencies.  We 
have  utilized  both  time-dependent  and  steady-state  analyses 
of  the  Master  Equation  while  incorporating  an  experimentally 
consistent  set  of  exchange  rates.  Within  this  computational 
framework,  we  have  critically  evaluated  three  different  dis- 
sociation models  and  established  the  significant  influence 
of  rotation,  vibrational  depletion,  equal  energy  transi- 
tions, V-V  exchanges,  and  electronic  states  on  the  dissoci- 
ation reaction.  The  scope  of  this  analysis  has  been  broad, 
including  calculations  for  Hj , 0_ , Nj,  CO,  Br^,  and  and 
examination  of  both  dissociation  and  recombination  coeffi- 
cients. In  our  parameterization  of  dissociation  under  non- 
equilibrium conditions,  we  have  considered  excitation  due 
to  electron  impact  and  optical  pumping  and  evaluated  the 
relative  importance  of  electron  impact  and  thermal  dissoci- 
ation. 


w 


n 

New  Aspects  of  Dissociation  for 
Near- Equilibrium  Conditions 

The  emphasis  of  the  near  equilibrium  aspect  of  this 

study  has  been  directed  toward  the  critical  evaluation  of 

collisional  theories  of  the  thermal  dissociation  of  diatomic 

molecules.  Semi-classical  exchange  probabilities  based  on 

a distorted  wave  treatment  were  utilized  in  time-dependent 

and  steady-state  solutions  of  the  kinetic  equations.  The 

dissociation  process,  in  pure  diatomic  gases  and  diatomic 

gases  highly  diluted  in  an  inert  gas,  was  parameterized 

over  a wide  range  of  temperatures  and  gases.  In  this 

parameterization,  to  evaluate  the  rotational  contribution 

to  dissociation,  three  kinetic  models  of  dissociation  were 

considered : 

1.  Standard  Ladder 

2.  Coupled  Model 

3.  Uncoupled  Model 

In  addition,  the  implications  of  equal  energy  transitions 
were  examined  within  the  framework  of  the  rotational  models, 

(2)  and  (3) , and  these  results  contrasted  with  the  more 
standard  "over-the-top"  mechanism. 

A time-dependent  solution  of  the  Master  Equation 
established  that  the  VED  rapidly  attains  a pseudo-steady- 
state,  Figs.  5-2  and  5-3.  This  steady-state  distribution 
deviates  substantially  from  the  equilibrium  distribution 
because  of  vibrational-dissociative  coupling.  The  resulting 
depletion  of  population  substantially  reduces  the 

1 
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dissociation  rate  from  that  which  would  be  obtained  if  the 
coupling  were  not  considered.  Induction  times  were  calcu- 
lated for  the  various  models.  No  significant  variation 
in  induction  times  was  observed  among  the  dissociation 
models  considered.  We  concluded  that  this  indicated  that 
the  induction  time  is  primarily  controlled  by  the  relaxa- 
tion or  exchange  rates  and  is  relatively  insensitive  to  the 
vibrational  biases  of  the  dissociation  process  considered 
here.  A comparison  of  predicted  and  experimental  induction 
times.  Fig.  5-5,  showed  good  agreement  with  experiment 
except  at  very  high  temperatures.  At  temperatures  in  excess 
of  10,000°K,  the  calculated  incubation  times  were  too  long. 
We  concluded  that  this  deviation  at  high  temperatures  can 
be  attributed  to  a more  rapid  relaxation  due  to  the 
increased  importance  of  multiple  quantum  V-T  transitions 
in  this  temperature  range.  It  is  only  in  this  temperature 
range,  however,  that  collisions  are  sufficiently  impulsive 
that  multiple  quantum  transitions  need  be  considered.  When 
sufficient  transition  rate  data  becomes  available,  a realis- 
tic attempt  at  substantiating  this  hypothesis  can  be  made. 

The  quasi-steady-state  dissociation  coefficient  was 
calculated  for  the  three  models  and  compared  to  experimental 
data  for  mixtures  of  a diatomic  gas  highly  diluted  in  argon. 
Gases  considered  included  , 02,  N2 , CO,  Br2,  and  I2*  From 
this  comparison,  we  establish  three  significant  points. 
First,  the  standard  ladder  model  is  inadequate  and  yields 
dissociation  coefficients  that  are  much  too  low.  Second, 
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models  including  the  influence  of  rotation  predicted  larger 
dissociation  coefficients  than  those  obtained  from  the 


standard  ladder  model.  These  coefficients  exhibited  sub- 
stantial agreement  with  experiment  for  all  gases  considered. 
Within  the  rotational  models,  equal  energy  transitions  were 
found  to  yield  an  additional  increase  in  the  dissociation 
rates;  (the  increase  was  especially  marked  in  Br2  and  I2) , 
and  decrease  the  activation  energy  in  comparison  to  the 
over-the-top  model.  Agreement  with  experiment  was  further 
improved  by  their  inclusion.  This  latter  observation  is 
consistent  with  the  recent  results  of  trajectory  calcula- 
tions for  dissociation  in  (ref  92) . 

Dissociation  in  pure  diatomic  gases  was  then  considered 
in  an  effort  to  establish  the  cause  of  the  observed  rate 
enhancement.  This  analysis  revealed  that  the  observed  rate 
enhancement  can  not  be  attributed  to  changes  in  collisional 
efficiencies  based  upon  mass  variation.  However,  V-V 
exchanges  lead  to  rate  enhancements  similar  in  magnitude 
and  temperature  dependence  to  those  experimentally  observed. 

Having  established  the  predictive  capabilities  of  the 
rotational  models  regarding  the  dissociation  coefficients 
in  dilute  and  pure  gases,  a comparison  of  the  calculated 
recombination  coefficients  with  experiment  was  presented. 
This  comparison  expands  the  temperature  range  over  which 
the  models  can  be  evaluated  and  provides  a detailed  examina- 
tion of  the  pre-exponential  factor  in  the  dissociation 
coefficient.  Reasonable  agreement  with  experiment  was 


f I 


achieved  for  all  gases  considered,  demonstrating  an  ability 
to  quantitatively  correlate  very  high  and  low  temperature 
rate  measurements  for  a wide  variety  of  gases.  The  equal 
energy  transition  model  was  found  to  yield  not  only  a more 
reasonable  magnitude  for  the  recombination  coefficient,  but 
also  a temperature  dependence  consistent  with  experimental 
data  over  the  entire  temperature  range. 

As  a consequence  of  this  aspect  of  the  study,  we  reach 
the  broad  conclusion  that  collisionally  induced  thermal 
dissociation  rates  of  diatomic  molecules  can  be  accurately 
calculated  using  a simple  molecular  model  which  includes 
both  rotational  structure  and  equal  energy  transitions. 

Non-egui libr ium  Dissociation — Heavy 
Particle  Vs . Electron  Impact 

An  examination  of  dissociation  under  highly  non- 
equilibrium conditions  was  utilized  to  illustrate  the 
effect  of  non-equilibrium  on  a simple  chemical  process.  A 
parametric  mapping  of  heavy  particle  and  electron  impact 
induced  dissociation  was  performed.  For  the  vibrational 
kinetic  analysis,  we  utilized  the  Master  Equation  formalism 
previously  developed  for  the  near  equilibrium  case.  The 
electron  impact  analysis  of  dissociation  was  based  on  a 
steady-state  solution  of  the  collisional  Boltzmann  equation. 

The  electron  impact  analysis  showed  that  significant 
dissociation  occurs  under  normal  discharge  operating  condi- 
tions. The  magnitude  of  the  dissociation  rate  exhibits  a 
strong  dependence  on  E/N,  and  is  substantially  enhanced  due 
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to  the  influence  of  superelastic  vibrational  collisions. 
Within  this  analysis,  the  influence  of  dissociation  from 
excited  vibrational  states  was  considered.  A considerable 
increase  in  the  rate  of  dissociation.  Fig  6-18,  was  noted. 
The  magnitude  of  this  enhancement  reflected  the  highly  non- 
equilibrium nature  of  both  the  VED  and  eed  under  normal 
discharge  operating  conditions. 

The  dissociation  rate  due  to  heavy  particle  collisions 
was  examined  over  the  discharge  parameter  space.  A time- 
dependent  solution  of  the  Master  Equation  established  the 
existence  of  a quasi-steady-state  for  the  VED  under  the 
non-equilibrium  conditions  considered  and  revealed  the 
temporal  development  of  the  three  characteristic  regions 
of  the  VED.  The  calculation  was  benchmarked  against  the 
experimental  data  of  Rich  (ref  133)  and  showed  good  agree- 
ment with  experimentally  inferred  distributions.  The  pre- 
dicted heavy  particle  dissociation  rate  was  compared  with 
the  electron  impact  rate.  Electron  impact  dissociation 
was  determined  to  substantially  exceed  the  heavy  particle 
rate  for  typical  discharge  operating  conditions.  In  order 
to  complete  the  parameterization  of  the  heavy  particle  rate, 
excitation  achieved  by  means  of  intense  optical  pumping  was 
considered.  Based  upon  this  analysis,  within  the  framework 
of  the  ladder  model,  we  concluded  that  it  is  difficult  to 
achieve  substantial  heavy  particle  collision  induced  disso- 
ciation under  non-equilibrium  conditions.  This  conclusion, 
however,  was  inconsistent  with  Rich's  experimental 
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observation  of  substantial  dissociation  in  CO.  We,  there- 
fore, decided  that  a more  efficient  reaction  channel  for 
CO  dissociation  must  exist.  A multi-quantum  V-V  exchange 
mechanism  and  a V-e  exchange  mechanism  were  proposed  and 
were  shown  to  be  consistent  with  available  experimental 
data.  Experiments  were  suggested  to  evaluate  the  effec- 
tiveness of  each  of  the  proposed  exchanges. 

Considerations  for  Future  Study 

In  the  course  of  this  study  we  have  examined  the  dis- 
sociation process  under  equilibrium  and  non-equilibrium  con- 
ditions. As  a result  of  these  analyses,  the  following 
broad  areas  can  be  identified  as  warranting  further  theo- 
retical and  experimental  investigation: 

a.  Heavy  particle  energy  transfer  rates 

b.  Excited  state  reactions 

These  areas  are  not  mutually  exclusive,  but  are  generically 
different  in  the  way  they  affect  the  dissociation  process. 
The  heavy  particle  transfer  rates  form  the  kinetic  basis  for 
a Master  Equation  approach.  Whereas,  excited  state 
reactions  dictate  the  critical  reaction  channels  to  be  con- 
sidered. 

Significant  strides  have  been  made  in  the  theoretical 
field  of  heavy  particle  kinetics  using  semi-classical  or 
exact  quantum  mechanical  calculations  of  inelastic  colli- 
sion processes.  Much  of  this  work  is  based  on  the  avail- 
ability of  improved  interaction  potentials.  The 
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experimental  and  theoretical  determination  of  these  poten- 
tials is  critical  if  accurate  predictions  of  energy  transfer 
rates  are  to  be  achieved. 

New  experimental  techniques  (refs  133,  161)  are  comple- 
menting and  capitalizing  on  the  recent  theoretical  predic- 
tions, however,  the  following  rate  deficiencies  still  exist: 

1.  There  is  a lack  of  sufficient  experimental  data 
relating  to  the  dependence  of  V-T,  V-V,  V-R-T,  V-R-V-R-T, 
vibrational  rate  coefficients  on  the  initial  state  quantum 
number . 

2 . Multiple  quantum  V-V  and  V-T  exchanges  have 
received  little  critical  consideration.  In  view  of 
Brechniacs'  recent  results  and  their  implications,  further 
work  should  be  devoted  to  this  topic. 

3.  Finally,  vibrational  to  electronic,  V-e,  or  elec- 
tronic to  vibrational,  e-V,  transfers  should  be  given  atten- 
tion commensurate  with  the  unique  possibilities  they  offer. 
The  possibilities  of  new  lasers  based  on  electronic  transi- 
tions, isotope  separation,  and  new  energy  storage  concepts 
using  metastable  states  should  spur  new  research  in  this 
area. 

The  experimental  observation  (refs  105,  108)  and  our 
numerical  confirmation  of  the  participation  of  excited 
electronic  states  in  CO  dissociation  process  suggest  that 
further  efforts  and  new  techniques  be  applied  to  establish 
the  role  of  excited  states  in  the  dissociation  process. 

Under  non-equilibrium  conditions,  we  have  established  that 
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excited  vibrational  and  electronic  states  can  significantly 
alter  reaction  rates,  discharge  chemistry,  and  laser  effi- 
ciency. Further  work  is  required  in  this  area  to  establish 
excited  state,  vibrational  and  electronic,  electron  impact 
excitation  and  ionization  cross-sections.  Detailed  knowl- 
edge of  these  cross-sections  offers  the  opportunity  to  con- 
trol electric  discharge  chemistry,  optimize  discharge  effi- 
ciency, and  temporally  tailor  discharge  excitation. 
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Appendix  A.  The  Influence  of  Anharmonicit^ 
on  Vibrational  Energy  Transfer 


In  this  appendix,  we  briefly  review  the  semi-classical 
theory  of  vibrational  energy  transfer,  V-T,  and  use  this 
theory  as  a framework  to  discuss  the  influence  of  anharmon- 
icity  on  the  transition  rates.  This  discussion  is  based  on 
the  recent  work  of  McKenzie  (ref  32)  and  the  earlier  papers 
of  Mies  (ref  27) . From  this  discussion  and  subsequent 
examination  of  calculated  matrix  elements,  we  will  conclude 
that  when  anharmonic  effects  are  considered,  two  factors 
must  be  incorporated  into  the  transition  element  analysis: 

a.  Morse  matrix  elements 

b.  Assymetric  correction  factor 

When  both  factors  are  properly  evaluated,  the  scaling  of 
transition  elements  is  found  to  deviate  only  slightly  from 
that  of  a simple  harmonic  oscillator. 

Let  us  now  consider  the  simple  problem  of  determining 
the  probable  state  of  an  oscillator,  initially  in  a speci- 
fied state,  after  it  has  been  perturbed  by  undergoing  a 
collision  with  a structureless  atom  (ref  32) . Here  we  will 
utilize  the  simple  collision  geometry  given  in  Fig.  A-l 
where  a structureless  particle  of  mass,  in  , impacts  a 

a 

diatomic  molecule  with  nuclear  masses,  m^,  and  mc»  The 
impacted  oscillator  nucleus,  m^,  is  located  a distance  yr 
from  the  molecular  mass  center,  where  y=mc/ (n^+m^,)  . We 
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Coilinear  Collision  Geometry 

Fig  A-l  Coilinear  Collision  Geometry 


adopt  a three  body  center  of  mass  reference  frame  in  which 
the  relative  collision  speed  is  given  by  u,  and  the  inter- 
action potential  is  specified  by 


V(^  =A 


-*/L 


(A-l) 


where  L and  A are  constants.  Defining  the  potential  in 
terms  of  the  oscillator  coordinates  and  the  center  of  mass, 
we  obtain 


. - % *r/L 


The  system  Hamiltonian  is  given  by 


z -z 


(A- 2) 


with 


//<>=  'n6mc/(/nb+r»c) 


(A- 4) 


(a-5) 


The  symbols  with  subscript  0 denote  oscillator  quantities 
while  the  other  symbols  refer  to  the  incident  particle 
variables.  Let  us  write  the  total  wave  function  as 


(A-6) 


and  treat  the  incident  particle  classically.  Then  the  equa- 
tions of  motion  for  the  incident  particle  yield  (ref  32) 


-£***<•'*> 


(A- 7 ) 


with 


yr/L>=7^  tUyr/L^(r,t) 


(A- 8) 


Note  that  the  potential  constant.  A,  only  influences  the 
distance  of  closest  approach,  and  may  therefore  be  removed 
by  the  transformation 


. £ - *Jl 


(A- 9) 


where  E=*syu  is  the  semi-classical  relative  collision 
energy  before  the  interaction,  and 


fa  WCr)eYr/L%(r)  Jr 


(A-10) 


is  the  time- independent  diagonal  matrix  element,  and  (r) 
the  initial  oscillator  stationary-state  eigenfunction. 

Defining  a new  interaction  coordinate  as  Z=x-xo,  and 
substituting  into  equation  (A-7)  yields  (ref  32) 


£ " 


y*  ^ = -r  e RCt) 


(A-ll) 


where  R (t)E<eYr/^L> /V^.  If  the  oscillator  is  assumed  to 
be  fixed  in  its  initial  eigenstate,  then  R(t)~l  for  all  time 
(ref  32) . This  assumption  reduces  the  classical  equation  of 
motion  to  the  equation  for  a constant  energy  trajectory: 


/jl 


(A-12 ) 


This  equation  can  be  integrated  analytically  (ref  170)  to 
obtain  the  interaction  potential 


V(t,r)  = K e-Zft,Vr/" 


(A-13) 


in  terms  of  time,  by  using  the  result 


e Z™/l 


(A-14) 
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The  oscillator  motion  is  treated  by  expanding  its  wave 
function  4i(r,t)  in  terms  of  the  stationary-state  oscillator 
eigenfunctions  \J>n(r),  thus 


= ZL  cnCt)  ) 


(A-15) 


with  a)  =E  /R  and  E the  eigenenergy  of  state  n.  The  wave 
n n n 

function,  \p(r,t)  describes  the  oscillator  response  during 
a collision  and  is  the  solution  of  the  time  dependent 
Schrodinger  equation. 


it  Jjfpt)  = j7jL-£-+V(r)+  (A.16) 


where  VQ(r)  is  the  molecular  potential.  The  solution  of 
equation  (A- 16)  is  reduced  to  a set  of  coupled,  linear  dif- 
ferential equations  for  the  expansion  coefficients,  Cn« 

By  defining 


vn: 


(A- 17) 


and  introducing  equation  (A-14) , the  coefficients  are  deter' 
mined  to  vary  in  time  as 


i*  e‘ * Ct)/L  Tc/fi)e<H'‘i')V . (A-18) 

fi  v 

The  probability  that  an  oscillator,  initially  in  state  k, 
will  exist  in  state  n after  the  interaction  is 


(A-19) 


If  only  the  initial  and  final  states  are  coupled,  equation 
(A-18)  can  be  written  in  the  integral  form  (ref  32) : 


where 


(«^)<Wkk)Efak  SzchXT r) 


(A-21) 


Mies  (ref  27)  has  integrated  equation  (A-20)  and  obtained 


P (£) - /tt-  2, * 

A+n  / VtA  smhfrf)  P > 7 (a-: 


where 


L C^'WA)  2-  fVnn~V**  ) (A- 2 3 ) 

t * > * ( W 

and  M(l+ig,  2,  i2X)  is  the  confluent  hypergeometric  series 
with  complex  arguments. 

Following  Mies  (ref  27),  we  redefine  this  probability 
as 


lA’24) 


where 

Vfij* 

Then  A(E,X)  can  be  identified  as  an  anharmonic  correction 

factor  to  the  "symmetric"  transition  probability,  Pk+n' 

where  the  symmetric  transition  probability  is  obtained  by 

assuming  V,.=V  . Values  of  A(E,X)  for  0-*-l  transition 

kk  nn 

have  been  tabulated  by  Mies  (ref  27)  for  02,  N2,  Br2»  and 

-1  -2 

CO.  Values  of  A(E,X)  range  from  10  to  10  . Thus  a 

general  effect  of  removing  the  symmetric  approximation  and 
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including  oscillator  anharmonicities  is  to  decrease  the 
transition  probability. 

Mies  (ref  27)  also  calculated  the  anharmonic  factor, 

A (E , X)  for  a limited  range  v=0^4  of  initial  and  final 

oscillator  states.  He  established  that  for  single  quantum 

transitions,  A(E,X)  decreases  as  the  vibrational  quantum 

number  of  the  initial  state  increases.  To  investigate  this 

scaling  in  more  detail,  we  have  expanded  the  range  of  states 

investigated  to  v=0-*40  for  the  CO-Ar  system.  The  results 

of  this  analysis  are  presented  in  Fig.  A-2,  where  we  plot 

2 

the  normalized  transitioned  elements  <VIF>  > f°r  F=I+1,  where 


for  two  cases : 

a.  Morse  Oscillator 

b.  Simple  Harmonic  Oscillator 

For  the  SHO  case,  we  have  retained  only  the  first  two 
terms  in  the  expansion  of  equation  (A-2) . This  yields  the 
frequently  assumed  form  for  the  interaction  potential  being 
linear  in  r.  The  Morse  elements  were  evaluated  according 
to  equation  (A-17)  using  the  formalism  and  numerical  tech- 
niques developed  by  McKenzie  (ref  32) . 

Note  that  the  Morse  elements  scale  much  more  rapidly 

vibrational  quantum  number  than  the  linear  behavior 
••  • « for  the  SHO  case.  However,  when  the  Mies 


Value  of  Transition  Element 


Fig  A- 2 Transition  Probabilities,  V-T 
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assymraetrical  correction  is  applied,  dashed  curve,  we 
observe  that  the  effective  transition  elements  scale  less 
rapidly  than  those  of  the  SHO. 

We,  therefore,  conclude  that  to  properly  identify  the 
magnitude  and  level  scaling  of  the  transition  elements  we 
must  consider  both  the  Morse  matrix  elements  and  the  Mies 
assymmetric  correction  factor. 
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Appendix  B.  Number  of  Bound 
Vibrational  States 


In  any  collisional  theory  of  energy  transfer  and  dis- 
sociation, it  is  necessary  to  enumerate  all  the  discrete 
vibrational  or  vibrational-rotational  states  to  be  con- 
sidered. There  exists  some  inconsistency  in  the  specifica- 
tion of  these  states.  Various  approaches  to  their  specifi- 
cation will  be  discussed  here. 

The  simplest  molecular  model  is  that  of  a simple 
harmonic  oscillator  with  energy  level  spacing  given  by 

hv  • The  number  of  bound  vibrational  states  for  this  model 
e 

is  simply  (DQ/(hve) 

In  a more  refined  approach,  one  would  attempt  to 
select  a model  which  at  least  accounts  for  the  molecular 
anharmonicity . The  Morse  potential,  Fig.  4-1,  is  a widely 
used  potential  curve  for  an  anharmonic  oscillator.  Morse 
proposed  that 

V- (/-e  *(r  ) (B-i) 

with  De  and  a being  constants,  (see  Chapter  IV),  related 
to  the  dissociation  energy  and  vibrational  frequency,  vg, 
of  the  ground  state.  Within  this  representation,  the 
vibrational  energy  levels  are  given  by  (ref  78) 
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(B-2) 


Ev-(v*my,  - (JsJ/C(Dt)(v<‘/zf 


Solving  for  the  last  bound  vibrational  state  v*,  by  setting 
Ev#=De  we  obtain  with  y=(v+*s) 


. . , I At-  4/ni^k 

zWi/(4vr) 


2^/CH) 


(B-3 ) 


Thus  in  a Morse  formalism,  the  number  of  bound  vibrational 
states  is  twice  that  obtained  using  the  truncated  SHO 
formalism  (ref  78) . We  have  adopted  the  Morse  formalism 
in  this  analysis  because  it  accurately  reproduces  the 
observed  spectroscopic  dependence  of  the  vibrational  energy 
structure  and  is  amenable  to  a direct  incorporation  of 
rotational  effects  (ref  77) . 

In  general,  anharmonicity  may  be  accounted  for  by 
expressing  the  vibrational  energy  levels  in  terms  of  a 
power  series  in  y: 


(B-4) 
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the  value  of  the  coefficients  of  the  expansion  are  given 
by  Herzberg  (ref  78).  If  this  formalism  is  adopted,  rather 
than  the  Morse,  and  the  expansion  truncated  at  the  quad- 
ratic term,  y*  is  given  by  the  smaller  of  the  two  roots  of: 

Thus 


Atj  - UVS-  4^'X'D' 


and 


Lj*x^  = 1/2-  )'/i  ” (B"5) 

Thus  in  this  formalism  y*<Y*Morse*  As  an  example,  consider 
oxygen.  Using  the  data  of  Herzberg  we  obtain  y*Morse  = 

53  and  y*=35. 

Preliminary  investigations  of  the  sensitivity  of  the 
dissociation  rate  to  state  specification  have  yielded  the 
following  general  conclusions  for  the  case  of  thermal  dis- 
sociation. 

a.  The  temperature  dependence  of  the  dissociation 


rate  is  relatively  insensitive  to  the  formalism  selected. 


b.  The  magnitude  of  the  dissociation  rate  is  smaller 
when  additional  vibrational  states  are  assumed  bound. 

This  variation,  however,  was  not  significant  enough  to 
alter  the  general  agreement  and  conclusions  obtained  with 
the  Morse  formalism. 

In  non-equilibrium  studies,  the  dissociation  results 
are  much  more  sensitive  to  the  assumed  number  of  bound 
states.  This  is  just  another  aspect  of  the  rate  sensitiv- 
ity analysis  that  was  discussed  in  Chapter  VI.  By  using 
the  power  series  expansion  with  a sixth  order  fit,  no  sig- 
nificant variation  from  the  Morse  value  in  calculated  bound 
states  was  obtained  in  CO.  However,  if  oxygen  or  nitrogen 
are  considered,  the  variations  are  significant.  For  example, 
using  the  truncated  series  approach  one  would  conclude  that 
in  oxygen  significant  heavy  particle  dissociation,  via  the 
ladder  model,  can  occur.  This  question  of  last  bound  state 

r j | 

designation  and  rate  scaling  requires  additional  experimental 
and  theoretical  attention. 
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b.  The  magnitude  of  the  dissociation  rate  is  smaller 
when  additional  vibrational  states  are  assumed  bound. 

This  variation,  however,  was  not  significant  enough  to 
alter  the  general  agreement  and  conclusions  obtained  with 
the  Morse  formalism. 

In  non-equilibrium  studies,  the  dissociation  results 
are  much  more  sensitive  to  the  assumed  number  of  bound 
states.  This  is  just  another  aspect  of  the  rate  sensitiv- 
ity analysis  that  was  discussed  in  Chapter  VI.  By  using 
the  power  series  expansion  with  a sixth  order  fit,  no  sig- 
nificant variation  from  the  Morse  value  in  calculated  bound 

states  was  obtained  in  CO.  However,  if  oxygen  or  nitrogen 

l 

are  considered,  the  variations  are  significant.  For  example, 
using  the  truncated  series  approach  one  would  conclude  that 
in  oxygen  significant  heavy  particle  dissociation,  via  the 
ladder  model,  can  occur.  This  question  of  last  bound  state 
designation  and  rate  scaling  requires  additional  experimental 
and  theoretical  attention. 
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